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Stresses in a Semi-Infinite Plate Containing 
a U-type Notch Under Uniform Tension 


: e, By M. Seika 


| 1. Introduction. The effect of a semi-circular notch on a semi-infinite plate under a uniform 
_ tension was first studied by Maunsell! through the use of polar coordinates. The study was later 
- extended by Ishibashi?, Weinel®, Ling*, and Udoguchi® with the aid of bipolar coordinates to cover 
_ a single notch of the shape of a circular arc. The case of an elliptic notch or a rectangular one 
_ with rounded corners in a semi-infinite plate under a uniform tension was discussed by Higuchi 
_ and Suzuki® and by Suzuki’, respectively, by means of the complex-variable method and a method 
__ of successive approximation. 

__ In the present paper, the problem of determining the effect of a single U-type notch on a semi- 
_ infinite plate under a uniform tension applied parallel to the straight boundary will be treated by 
_ the complex-variable method, associated with the name of Muskhelishvili8. The parametric coeffi- 
cients included in the solution will be determined by a method of perturbation. The shape of the 

notch was introduced by Greenspan®. By varying the parameters in the mapping function, some 
_ other shapes of practical interest may also be obtained, which include the semi-circle, the semi- 
_ ellipse, and the rectangle with rounded corners; the solutions to be obtained therefore can also be 
_ applied to all these shapes. The results when properly reduced will be seen to agree with those 

given by previous authors. 


2. The method of analysis. Let us consider an isotropic homogeneous semi-infinite plate with 


a single U-type notch subjected to a uniform tension T acting parallel to its straight boundary, 
Fig. 1. 
Now, we shall introduce the mapping relation due to Greenspan® 


s=a0) =RC+E +3} (R> 0, 1> |p| =0, 1> |q| = 9), (1) 


PR ‘aie, | 


where R, p and q are real constants. Denoting z = x + iy and ¢ =o e'®, we have 


+= R(o +f) cos O +z 00530 ; y= Rl(o—F)sin@ — Fsin3 0 5; (2) 


_ When 0 = 1, by varying the values of the constants R, p and q, equation (1) maps various shapes 
__ of the opening in the z-plane onto the unit circle in the ¢-plane. The center of the opening is map- 
__ ped into the center of the unit circle. Some shapes of practical interest are as follows: (a) if p =q =0, 
 acircle of radius R is defined; (b) if g = 0, an ellipse is defined; (c) if p = 0 and q is positive, a 
- square with rounded corners of which diagonals are on the x- and y-axes can be approximated; 

(d) if p = 0 and q is negative, the preceding “square” is rotated 45°; (e) if p == 0 and q + 0, equa- 

tion (1), by proper choice of these constants, can be made to approximate an ovaloid opening, 
which is made up of a square and two semi-circles. 


1 FG. Maunsell, Phil. Mag. 21 (1936) p. 765. 

2 T. Ishibashi, Mem. Faculty of Engrg., Kyushu Imp. Univ. 9 (1940) p. 131. 

3 E. Weinel, Z. angew. Math. Mech. 21 (1941) S. 228. 

4 C. B. Ling, J. Math. Phys. 26 (1947) p. 284. 

5 H. Udoguchi, Trans. Japan Soc. Mech. Engrs. 16, No. 55 (1950) p. 44 (in Japanese). 

6 S, Higuchi and M. Suzuki, Technol. Rep. Tohoku Univ. 14 (1949) p.95. 

7 M. Suzuki, Reports of the Laboratory for Stress and Elasticity, Tohoku Uniy., No. 2 (1949) p. 61 (in 
Japanese). 

8 N. I. Muskhelishvili, Z. angew. Math. Mech. 13 (1933) S. 264. 

9 M. Greenspan, Quart. Appl. Math. 2 (1944) p. 60. 
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Now, we apply the mapping function (1) in the notched problem. The edge of the notch is” 
taken as unity (9 = 1), and the middle surface of the plate treated here occupies the right half 


of the z-plane as shown in Fig. 1. : 
Let the plate be in a state of generalized plane stress. Then the stresses, 


when averaged across the thickness, are derivable from a single stress func- 
tion F’, satisfying 


Ahn (3) 


The solution of this biharmonic equation can be expressed in terms of two 
analytic functions y,(z) and y,(z) of the complex-variable 2: 


F=Re EQ) +n@), ne =Sy@)&. (4) 
Overbars are used to indicate conjugate complex quantities and Re signifies 
“the real part of”. 
The stresses, in terms of the two functions, are 
0, + oy = 2[D,(2) + D,@)] . 


oi ek 6) 
where @,(z) = gilz), Py(z) = yilz) - 


Fig. 1. 
Notched plate under tension. 


By use of the mapping function z = w(¢), the stresses referred to the | 
curvilinear coordinates (0, 0) will be expressed as | 


0, +09 =2[G(0) + OO], 


Co — 0 ti2 te = om LE) OO) +o OVO, | 
where 9(C) = a2) =% {co(C)} > 0) =W%) =~{old)}. 
ee. _ ¥@) = _v) 
D(C) — P, (2) ee w’(£) ? YC) aE (2) =e w’(E) ° J 


The boundary conditions for the present problem are given as follows. 
[1]. The plate subjected to a uniform tension T in the direction of y-axis: 


109; oy = T, Care (7) | 
[2]. On the straight boundary, free from loading: | 


O=455 +r 9O+HO =0. (8) 
[3]. At the rim of the stress-free notch, denoting e'? =o: | 
ya heey | 

g=1; 90) sn? e +9 (3) =i) (9) 


Now the problem is reduced to the construction of two analytic functions y() and y(¢) which > 
satisfy the boundary conditions (7), (8) and (9). 


3. Determination of the parametric coefficients. We shall consider the functions 
92) = 90) + 9%), —-wZ) = yO} + OOD), (10) 


ee g(C) and y)(¢) satisfy the conditions (7) and (8), and can be written in terms of 2(x, y) coor: 
inates as 


1 1] 
PQ) ={7Tz, ys) =— Ts. , (11) | 


Through the use of the mapping function given by equation (1), we transform the functions in equa- 
tions (11) into 


I 
P= ETRE EE a, YC) =a TERETE +a). (12)) 
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Equations (12) are the solutions in the case where the notch is absent, and they cannot satisfy the 
condition (9). In order to satisfy this condition, p™(¢) and y(¢) are introduced. 
Now we may write ~)(C) and w()(¢) as follows: 


(1) = — . peseine —2 (n + 1) TS —(2n 
PUG ae BOS] reatae ss coe) 


ye) = TRIED Dd) [on—Or49 + 45-2640) (13) 


Oe 
= ee ee) pee dn C—@n+) 
+ >, |e) In+l1 


Because of the symmetry of the stress distribution in the plate, the coefficients a,, b,, Cy, and d,, 
included in the foregoing functions are real. Evidently, the stresses derived from these functions 
tend to zero when €—> 00. 

Substituting equations (10), (12) and (13) in the condition (8), we have the relations 


Gerster, [i Obs (14) 


| 
| 


Substituting equations (10), (12) and (13) in the condition (9), using the relations (14), and then 
separating into real and imaginary parts, we obtain two expressions. Since these expressions 
contain cosines and sines of odd and even multiples of @, it is necessary to expand these cosines 
and sines in Fourier series of cosines and sines respectively, of odd multiples of O valid from 


0 = —2/2 toO =72/2. The required formulas are 
4 wy(—1)"*# 2s+1) 
cos 2n (0) a a a @ s+1?—4n2 - COS (2 Ss + 1) (0) P) (1s) 
see eas i Oe) On 


% del (28 + 1)?—4n? 

Taking into account these formulas, we finally have the following equations for a, and 6,: 

1 

le P29 3 pg — 2p’ —3 9°) o,. + + p—3q—5 pg) or, + (1— 39) ¢ 02,1 

2 BY 
FUL +4P +39) 8.5 Q) dua] bo +[(5 +P +24)dn0+(1+ FP) one] 
4 3 

+ [a600+(3 +P + $9) 810|bs +9510 b — (1 —d,) (1-31) 


i 24s 1) ¢ 
’ eo. Dee? 


2 3 2 1)+(2s+3)p+3q 
oan ).— ae + Bs Oenay ab.+2| 


oe) Sta, {2m + 3) [4 (wv + 2)? — (2s + 1)? — 8] 0. 0." 0," 
4 (2m + 1) pastas; —3 (2 —1) [4 (nw? +2) —(2 5 + 1)] gas {% ai} = 0 (8 = 0), (16) 
= 1 +p+2q+3pq—3@)o,.+(+p+349+3pq)o,s+ (1 +349) 4 02,6] 
+ (1 +34) 0,5 +61.) b+ |(g +P —24)d02 + +P) b1,4]6 
+ |0d..+(F +P —$4)O.0| be + 961.06 + —dr,.) 0d.) 
BGA) a sor 3) P34 

tt = 


1+2(s+1)q ¢ 
x | oF 46,1 +(1 +p) b+ ra 1+ qbe42 


> (— 1)" + * a, (2 n + 3) [4n(n 4+ 1) (n + 2) —(n—1) (28 + 1)] 0 0. a, 


at 
$n (2n +1) pai of — 3 (2 —1) 


x [4n (n— 1) (n+ 2)—(n +1) (284 1?) qari Bets ot 0 (sv 0), (17) 
20* 
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where 6; ; means Kronecker delta and 


) 1 
Wik= aE Qstok ln Cs =o heen Czy te) 


Eliminating b,, between equations (16) and (17), we obtain the infinite set of equations con- 
taining the unknowns a, only. This set of equations may be solved for a, by a method of suc- 
cessive approximation, for instance, when the values of p and q are given. Once a, have been 
determined, the unknowns 6, are to be found. 


4, Application of a perturbation method. It is easily seen by equations (16) and (17) that the 
unknown coefficients a, and b, are the functions of the constants p and q. Further, the constants 
p and q are the parameters which govern the shape of the notch. Hence, to solve the infinite set 
of equations (16) and (17), it will be convenient to apply a perturbation method in which p and q | 
are the perturbation parameter. 


Now we shall seek a solution to equations (16) and (17) in the form 


ee Ma q’ pt : b, = ma 2a 6H q’ p' : (19) | 


ine 


where “a() and “b are independent of p and q. Substituting equations (19) in equations (16) | 
and (17), and equating to zero the coefficients of the various powers of p and q, we obtain the | 


following expressions: 


orem 0handutea0)s 


Sy fe AY eo (0) (0) eS 
Be NS oa 
~ pee \ re ma) 0). 0) 2s 
BoA ee AGRE Oa ae prt (20) 
an ee pen 0 (s = 2), | 
n=0 
= n 5 
3 (—Wrabh Pa? = —F (811 + Gone + yon) 
n=0 


ia yy (— 1)" aon ay {(2 m — 1) [(2 m — 3) git—9) _ (2m + 9) Mgt-D] 
m=0 
+ (4m? — 16 m — 1) “at-?} , 


= Sd NC 9 (n) (n) (0) Oa ry , } ; | 
oh Sen SO Ce g(t 4 f=1] @a 


PEs es Ly [(2 m + 1) yao? BY (0) alm » + 4 (m — 1) ) af? Bo (0) Goa ine 
m=0 
Ss (— HN (en + 3) os) BS ae en Ss (—1) ™* 0%") [(2m +1 2a (f-—2) 
a0 m=0 
ci eee 
== [ny (271 — 2s 1) (Ge) ye an = al (siz 2) 


by A, Pe. 
SY = Xr = = » ( Lys (2 m— 1) (2 m ule 9) af ay” Caer 


m=0 


Be — fae ta 7) (— 1)" af eer on (22) | 


pO __ 8s ie 


Sy) (— Lye (2 m + 3) (2m—2s+ 3) as" ys" eae. (S22); 
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OHO = (3 61,4 + 62,1) 


4 co 
=F DS (HM oft oft {(2 m — 1) [2 m— 3) Pah? — 4 a) 


CO 


+ (4 m* — 16 m— 1) afi3 3, 


3 
) OTP et or 62, t ed Coie ae 9 Ob) 


—s 5 S (— 1)" alt? {a8 (2 m +1)? ait-? + 4m (2m + 3) at-?] [eer] 


10 a) 
Cpe Cees: pre 0) _ 8@s+)) os yak 3 


nA s 


m=0 


xX {(2m— 2s + 3)[(2 m + 1) at ee 24 (2m + 3) a oO] 
+ [(2 m + 3) (2m—2s+41)4+ (2m+4+1)(2m—2s + 5)] ai-?} 


($2). 
For r= 1 and t=0; 
~ n 1M) (),0) 3% TH r—wpo) = m v(m) (r—1),(0) 
Do (— 1) XO,1 an =Saeay Vie bY ar 3 oe (— 1) O11 Am > 
n=0 m=0 
> C1 On 4) as a, Pg BE MG a DRO 1 G-vpo 4. 2 r-ppo 
ar Coe EVENS 5 ae 28 7 
+5 SF rm +1) [Cm + 5) afk —@ m + V) af af Pa 
>) (— |i ek (2 n - 3) a™ pee 1) gO — pce 6 wel og om EES AO 
—. Ee io Ya aa ae) 2s+3 
352 oy 3 ~ mt+s m, m, 
4(2s +5) Sea aD) (2 m — 1) [22 m—2s + 3) Yn0 Yu —8 
—(2m—2s+4+7)y (m) ie Tats Sep (ge 72) 


co 3 2 2 
ee Dyan os a 62,7 (61,2 + 62, +) —= eget ifr = AD eae. 
n= 
3 
JF DAE-D 1 9 Dp t-D : O-DEE-D 4 Po) 


oe Feit 6 (4m? — 1) af? af? D git —2) 


710) 
=26 [o$"? pee: (4 m2 — 1) age "3 1—-Dgt—D Rie? Oi aD ge 
— af? af? {(2 m — 1) [(2 m — 3) Pak ® — (2 m + 9) Pat?) 
+ (4m? — 16m—1) Pam? }]> 
3 


1% 1 


Ze fo Li (2 n = 1) ay} as ae a =" & ay Bo Ale ae ee 
1 G—-Dpo 1 (r— Dp) 3 (r—-DpCt) 

es bo ae by seer bs 
a S (—1)" + aft (2 m — 1) {[(2 m + 3) af? — (2m — l) a ri Naame ait-? li : 

m= 0- = 1 
+ [(2 m + 5) of? — (2m + 1) a§2] 7-? ai?) 

) 

— (2 m + 1) aft BYP Pan” — 4 (m — I) ah Bo” Pam" | , 
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SAI Rn $3) als BP Pa = — 195d 07 31,5 
m [r-D_E=D 3 F -DHE=D (r—Dpt=D (F- DH 
8 14 Qs=) © 405409), | 2s 2eo) me era) 
aU De 1 3 REL. 
eres = 
a 3 m m 
~ ee macsare (ee 1) 7% yrs | 
—(2m—2s ae ©) yar Yor oe Dq Me 1) {2 m—2 st 3) yo) yor_s 
—(2 m—2s ai 7) my eo es D a?) a (2 m =e 1) a Iai Wer Mise 2) 
+ {[m Eee We ane | 
— |[m (2m — 2s + 1)—2 (s 1) p82} «4 ean" | (s 22), | 
| 
OHO — ae di ig DAO 4 Zr-DAO 4 Yo 9) 
a - > (1°43 [Bag 2 nt — 1) os of] Dg 
m=0 
=" (2 ni — 1) (2 Te a et ay Pq@r - 
3 
(DH = = on 2G i 140) ig = 2-9) 
cs (26) 
ee = Pa ye [3 (2 me 1) (2 m ais 5) age ag”? rn Sh a ay? at . 
7 2s+1 4(2s +1) 
PLO (r — 1) p00) (=) We 
ee ee ae ee AP a rece bs 1 
8 (2 iby 2 
SEO S Hunt BCm—Y m—28 + YY 26 Pay 
— (2m + 3) (2m—=—2 5 + 3) al} yo Ma] (32), J 
Cape) ae Sey) (3 Ort ae do, 1) ft (ones Je IK = rps” ] 
we 3 oe “ ae Naar ae Sey at 3 TaN bY? gine) 
sie =e i pe {6 (4 m? — 1) af yet as 1) aie 2) FG 
x es — (4 m2 — 1) (a8 + af) of] -Pat-P 
— 3 [3 af? — 2 (4 m2 — 1) af? af] 7 Pa? 
— ada aft [(2 m — 1) (2m — 3) Pak-® 4+ 2 (4m? — 12m 4 1) Pat-? 
Ba aa ea ice tes ONC 
Pee = 3 b1.-(81 = zoni)— 3(- DB a: 1) + Be eS DOS iD) + >, Gas eA 
Bee 1 co | 
SReeae ee cet! ee a yy’ (— I" {3 als @-VgC-D [t = 1] (27)) 
+ 3 (2 1m —1) (2 m + 5) o67 oft 7 Pa! — of aoe (2m +. 1)2 Pat-” 
-|- 4m (2 m aL 3) Oe es >| — ay” Oa aoe 
r 3 (2 1 
eS | oan ae isa Amgue ar a = 2 il) ear es = ae Sips De 
+ OD ames a 9) oe) 


hee 2) 


—+ 2G ae 3 yi Nee de m— 1) 


x i =) S13) a eens a ape ee 1) Sea) 
1 ¥s0 {(2m—2s + 3) [(2m +1) Pat? 4 (2m + 3) Pa] 
+ [@m4 1) (2m—2s +5) + (2m +3) (2m—2s + 1)] Pas-?}] 


(s =). | 


XXVII. Band 1960 M. Seika: Stresses in a Semi-Infinite Plate Containing a U-type... 291 


The symbols used in these expressions are the following: 0; ; means Kronecker delta and 


(mM __ i! 

PUES mn ea. -ekh + lm Os—2h- 1° 

(m) __ 1 

Bs ~ 2m+@s+)] m+ @s+5)° (s = 0) (28) 
Cn) aye 1 
Sse 1d) [am se kb | 


These equations are to be solved for the coefficients a‘? and ‘°b in the series (19) suc- 
cessively. To solve the infinite set of equations (20), (21), (24) and (25) for af’, the laborious 
method of elimination of unknowns has to be used in which the set is first reduced to finite by 
taking only the first ny equations involving the first ny unknowns. This will be performed finally 
for fifteen equations as treated by Maunsell‘, and the values of the first eleven ‘"a‘? will be used 
in the subsequent calculations. 

Tt may be noted that the foregoing method is valid in so far as the series (19) are convergent. 
The proof of convergence seems rather difficult to achieve. However, it may be seen later in the 


numerical examples that this method gives practically satisfactory results for certain values of p and q. 


5. Numerical examples. The numerical calculations have been carried out for the cases of 
elliptic notches and U-type notches. The values of og along the periphery of the notch and those 
of o, across the symmetrical axis of the plate have been given for two cases. 

a) Elliptic notches. If we choose 9 = 1 and q = 0 in equation (1), an elliptic notch can be 
obtained for a certain value of p. Fig. 2 shows the shapes of the elliptic notches treated here. The 
values of p, a/b and r/a corresponding to these shapes are tabulated in Table 1. 


Table 1. Values of p, a/b, and r/a 


Pp —0.20 —0.13 0 0.13 0.20 0.27 0.30 
a/b 0.66 0.77 1 13 1.5 IL) 1,86 
r/a 2.25 | 1.69 1 0.69 0.45 0.33 0.29 


In this case, the coefficients “a? and 6‘ are calculated by equations (20) ~ (23). Espe- 
cially, for the case of p = 0, i. e., a semicircular notch, the coefficients a’ and “b? are ob- 
tained by equations (20) and (22), and the stresses derived from these coefficients coincide 
with those given by Maunsell? and the others®. Previously, Higuchi and Suzuki* have dealt with 
the case of a/b = 1.3 (p = 0.13), and we can also obtain an agreement withtheir results. 

Fig. 3 shows the distributions of o, along the rim of the notch for several values of a/b. The 
distributions of o, across the symmetrical axis of the plate in relation to the ratio a/b are shownin 
Fig. 4. 

b) U-type notches. When the breadth of the notch 2 6 is constant in the preceding case, 
the radius of curvature at the bottom of the notch changes with the depth of the notch. Here, 
to clarify the effect of the change of the depth of notch only, we shall take the radius of curva- 
ture at the bottom of the notch equal to one half of the breadth of the notch. Fig. 5 shows the shapes 
of the U-type notches treated here. The values of p, qg, a/b and r/a corresponding to these shapes 
are tabulated in Table 2. The coefficients “a? and “b? required in this case are calculated by 


Table 2. Values of p, q, a/b, and r/a 


0.195 | 0.264 0.290 


p 0 | 0.069 | 0.128 | 

q 0 | —0.008 | —0.016 | —0.024 | —0.032 | —0.035 
ae P15 ies aero 1.75 1.86 
rla 1 0.87 0.77 | 0.67 0.57 | 0.54 


equations (20) ~ (27). Fig. 6 shows the distributions of og along the rim of the notch for several 
values of a/b. The distributions of 7, across the symmetrical axis of the plate in relation to the 
ratio a/b are shown in Fig. 7. 

The distributions of og and o, for this case are almost similar as those for the preceding case. 
But, the maximum value of og or o, corresponding to each value of a/b is smaller than the value 


F.G. Maunsell, 1. c. 

F. G. Maunsell, 1. c. 

T. Ishibashi, etc., l. c. 

S. Higuchi and M. Suzuki, 1. c. 


a  ~ 
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Fig. 2. Elliptic notches. 
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Fig. 6. og at the rim of the notches. Fig. 7. oy across the line of symmetry. 
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for the case a). This may be chiefly due to the fact that for the same value of a/b the radius of 

curvature at the bottom of the U-type notch is larger than that at the bottom of the elliptic notch. 
6. The stress concentration factor K. It is of practical interest to calculate the stress concen- 

tration factor K, or the ratio of the maximum stress at the bottom of the notch to the mean stress T. 
Now, the factor K for this case is readily computed by the expression 


4 
K=1 ieee ede (29) 


In the case of the elliptic notch, equation (29) is reduced to 
K=1+ i (0.5165 + 0.5435 p + 0.0421 p? + 0.0326 p? — 1.333 p* + 4.167 p5 


— 6.448 p® + 8.186 p? — 12.46 p® + 14.01 p® — 23.55 p!) , (30) 
where —0.20 S p < 0.30. Equation (30) may be approximated by a simple formulae: 
K = 0.941 + 2.125 — . (31) 


where 0.66 < a/b < 1.86. The factor K of the corresponding problem for the case of an infinite 
plate with an elliptic hole was given as! 


K=1+25 or K=1+42)/%. (32) 
The results obtained from equations (30) and (32) are tabulated in Table 3. 


Table 3. Stress concentration factor 


a/b Wrote ream PARE Sk ist ei ol 86 
Notch 235 | 258 | 3.06 3.70 | 4.13 | 4.66 | 4.90 
Hole 2.32 | 2.54 4 3 $6024 4 4.5 4,72 


In the case of the U-type notch, the coefficients “A in equation (29) are shown in Table 4, 
where 0.290 = p = 0 and —0.035 <q <0. The factor K for the case of an infinite plate having 
an oval hole is given as 


Table 4. Values of “A 


ay 0 1 | 2 | ey | 5 | 6 Tet fore 9 | 10 
| 
0 | 0.5165 | 0.5435) 0.0421 | 0.0326 1.333 4.167 | —6.448 8.186 | 12.46 | 14.01 [23.55 
1 | 2.675 | 2.828 |1.660 |—50.92 | 95.00; —140.9| 251.9] —441.0| 364.1 
2 | 6.614 [18.54 | 13.79 117.1 | 127.0| 403.2) 7244 | 
3 |15.64 |—85.08 |—13.02 |—2551 | 3582 | | | 
4 165.76 | 3482 | 2252 
2 
yee | p+6q aie (33) 


1—p—3q' 0—) —p—30) 
The results obtained from equation (29) by using Table 4 and from equation (33) are tabulated 
in Table 5. 


Table 5. Stress concentration factor 


ab ere kien 3 ae Bie | 1.86 
| 


Notch 


3.06 | 3.22 | 3.36 
Hole 


3 


These results obtained above are shown in Fig. 8, where I, the case for the elliptic notches; 
II, the case for the elliptic holes; III, the case for the U-type notches; and IV, the case for the 
oval holes. 


at CB. Inglis; Engineering 95 (1913) p. 415. 
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Fig. 9 shows the stress concentration factor K of the notch treated here in relation to the ratio 
r/a. The curve corresponds to the elliptic notches and © marks indicate the values for the U-type 
notches. It may be said that as far as the ratio r/a is taken from 0.54 to 1 the factor of stress con- 


centration is quite similar in these two cases. 


k 
Teal oni ane a eae 
6) 
I. 
7, et 
ME I 
IV Ca 
é} 
3 
‘Hh 
; a/b 
L_{|___|_ 
066 =—08 10 12 14 16 18 186 2G 05 10 145 20 Gar 
Fig. 8. Stress concentration factor. — I, Factor for the elliptic notches; Fig. 9. Stress concentration factor. Curve, Factor for the elliptic 
II, Factor for the elliptic holes; III, Factor for the U-type notches; notches; o marks, Factor for the U-type notches. 


IV, Factor for the oval holes. 


In conclusion, the author wishes to express his hearty thanks to Prof. Emeritus Dr. Seiichi 
Higuchi of Tohoku University for his valuable advice throughout the progress of the present in- 
vestigation. The author also wishes to express his gratitude to Prof. Dr. Osamu Tamate of Tohoku 
University for his kind discussions during the course of the work. 


(Eingegangen am 20. November 1958.) 


Anschrift der Verfasser: Assistant Masaichiro Seika, Sendai (Japan), Tohoku University, Department of 
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Zur Theorie der nahe der Wasseroberflache fahrenden Tragflichen 
Von W.-H. Isay 


1. Einleitung. Unterwassertragfliigel sind in den letzten Jahren mehrfach behandelt worden und 
zwar sowohl als raumliches Strémungsproblem!7 als auch in zweidimensionaler Betrachtungs- 
mrcige®? 4) 6, 7, 

Betrachtet man wie in der hier vorliegenden Arbeit die Strémung um parallel zur Wasserober- 
flache fahrende Tragfliigel als ebenes Problem, so liegt es nahe, hierfiir die in der gewéhnlichen 
Tragfligeltheorie bewahrte, auf Birnbaum zuriickgehende Methode der tragenden Wirbelflache zu 
verwenden. In den bisherigen Arbeiten ist dieser Weg nur bis zu einem bestimmten Grade be- 
schritten worden. Insbesondere wurde die Berechnung der Zirkulation des Unterwasserfliigels ent- 
weder gar nicht® oder unter Verwendung von naherungsweise giiltigen technischen. Auftriebsbei- 
wertrelationen durchgefiihrt*, und Ansatze zur Lésung der Integralgleichung beschranken sich auf 
bestimmte Sonderfalle® ©. 

In der vorliegenden Arbeit soll nun eine vollstandige Theorie angegeben werden, die auch den 
Fall von zwei hintereinander fahrenden Tragflachen umfaft (Tragfliigelboot). Es wird die zur 
Bestimmung der Auftriebsverteilung (Zirkulationsverteilung) notwendige Auflésung der Integral- 
gleichungen durchgefiihrt, und zwar in einer Form, die sich vollkommen an die aus der Schaufel- 
gittertheorie bekannten Methoden anlehnt. AuSerdem wird die Form der Wasseroberflache in der 
Nahe der Tragfliigel berechnet. Fiir eine Theorie der Unterwassertragfliige] mit Erfiillung der 
Randbedingung am Fliigelprofil unter Verwendung der Methode der konformen Abbildung sei 
auf Nishiyama’ verwiesen. 


2. Das Geschwindigkeitsfeld einer Tragflache unter Beriicksichtigung der freien Wasseroberflache. 
Wir betrachten eine Tragflache in der Tiefe h unter der bei y © 0 liegenden Wasseroberflache 
(Abb. 1) in der homogenen Anstrémung uy, wobei u und v die x- bzw. y-Komponenten der Ge- 
schwindigkeit bedeuten. An der Wasseroberfliche ist der 


Druck iiberall konstant. Aus der Bernouillischen Gleichung y 
erhalt man dann unter der Voraussetzung, dafi das Qua- bee e 
drat aller durch den Fligel auftretenden Stérgeschwindig- cies 
keiten u,, v, klein ist gegen up, die Oberflachenbedingung : Es ! 

| 


in der Form 


uyu, + g Y = konst. 
Dabei gibt Y(x) die genaue Form der Oberflache an. Da 


jedenfalls weit vor dem Fliigel fiir x ~ — © sowohl u,, als 
auch Y(x) genau Null sein miissen, ist auch die Konstante gleich Null, und man erhalt fiir die 


Form der Wasseroberflache 
Y(x) uz Uy l 
4 = ok (1) 

§ Uo 
mit up/V/g h als Froudesche Tiefenzahl. Aus (1) erhalt man weiter mit 0Y/0x = v,/u) die bekannte 
kinematische Bedingung der freien Oberflaiche 


s is ae 9 
ib Ot + 9, =0 (fiir y = 0) (2) 

Die betrachtete Tragflache wird mit einer Wirbeldichte y(&) (— a <€ < a) belegt; das wohl- 
bekannte von dieser Wirbeldichte induzierte Geschwindigkeitsfeld geniigt der Gleichung (2) natiir- 
lich nicht, sondern mu8 durch zusatzliche Geschwindigkeitsanteile erginzt werden. Im Hinblick 
auf die zu erfiillende Randbedingung (2) der freien Oberflache hat es sich als zweckmaBig erwiesen®, 


1 K. Krienes, Die tragende Flache in einer Stromung mit freier Flissigkeitsoberflache, unveréffentl. Be- 
richt Berlin 1951. 

2 Y. fT. Wu, J. Math. Phys. 33 (1954) S. 207. 

3 K. Krienes, Die von einem tragenden Wirbel an der Fiiissigkeitsoberflache hervorgerufene Wellen- 
bewegung, Bericht 8 (1946) des Konstruktionsbiiro RofBlau (unveroffentlicht). 

4 C, Possio, Atti Accademia Scienze Torino. 76 (1941) S. 365. 

5 H. B. Squire, Proc. Royal Soc. A 243 (1957) S.48. 

8 FE, Cumberbatch, J. Fluid Mechanics 4 (1958) S. 466. 

? T. Nishiyama, Society of Naval Architects of Japan; 60 th Anniversary Series, Vol. 2. (1957); S596 u.S. 116. 
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das von der Zirkulationsverteilung y(¢) induzierte Geschwindigkeitsfeld in der Form zu schreiben: | 


a 


tT peony deen ges So 
Uy ate (y+ hy [7@ fe Reins ey 2 ie | 
is) eo 6 i y>—h) Bal 


<5, | Gar eG ye wae fre if cH +0 sin Ge —8) da | 
a 


Fir die zusatzlichen Geschwindigkeitsanteile u,v? liegt daher der Ansatz nahe 


—a 


w(t) = Sh [re fF) ) ety cos (x —é) dd dé, 
very — b [ v(é) | F(A) ety sind (x —&) dd dé. 
—a 0 
Dabei ist die Funktion F(A) so zu bestimmen, da8 u, = ul) + ul) und v, = of) + vf) fir y = 0 
die Bedingung (2) erfiillen. Dies ergibt nach einer elementaren Rechnung 
2 e—Ah 


i eee 


& 


HQ) ee 


> 


und somit wird 


a 


SF — if, / Bie (k—y) 

= v(é) (h — y) dé & 1 e 0 & 

(I = RLS PRON See eA teen te ieee © (y—&) dud 
OD (h ae Ape ay ua ee 1 cos LL ua (x &) Lb & > 


| | jetoss 

1 v(&) (x — &) dé Bee Cae . g 

(I) — ete oS £8 (6, — 

Oy Bat (hy) Ee ee nat we) ae ee 
2a —a 0 


2 
Es ist dabei w = aA gesetzt. 


] 
| 
+ (A) 
| 
J 


Fiir die weitere Behandlung von (4) kénnen die zuerst von Krienes! angegebenen Integral- 
formeln 


co 


a oiniel "Ge? dd | 
[ice X du =e C,(X) cos X + S(X in an (PE al 
6 
| E 
eu H : = Xe dd 
[F2psn a Xu = a #]129 om x — C00 sa ee (ae | 
5 0 
| 


(+ fir X > 0; — fir X < 0) 


verwendet werden, deren Beweis wir in Ziff. 8 nachtragen. S; und C; sind Integralsinus und Inte- | 
gralcosinus. Fiir X = 0 ist im unteren Integral (5) der Wert Null zu nehmen. 

Nun darf jedenfalls weit vor dem Fliigel fiir « + —oo nur die ungestérte Anstrémung vorhan- | 
den sein, d.h. es muB gelten | 


Tima wateee=ritees lim’ v = 0. 6) 


x—>— 0 x—>— © 


TGicholbuBnote sven Seieeoos 
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Diese Bedingung wird durch unsere bisherige Lisung noch nicht erfiillt, denn wegen lim C(x) = 0 
und lim S,(x) = —z/2 ergibt sich ae 


os OO) 


Es muf somit noch ein weiterer Geschwindigkeitsanteil ul), fl) iiberlagert werden; dieser ist 


so zu bestimmen, daB er einmal fiir sich selbst der Bedingung (2) geniigt und auBerdem das ganze 
Geschwindigkeitsfeld 


u =u, + ul) + ul)+ ull), vy = off) + off) + oft) 


die Relation (6) erfiillt. Damit ist das Geschwindigkeitsfeld einer in der Tiefe h unter der Wasser- 
oberflache fahrenden Wirbelflache endgiiltig bestimmt; es lautet 


22 1 [ v(y+hjdé ai v(&) (y — h) dé 
Pato t a5 | @—) + +he! ia | («— £2 + (y—hp 


er —a 


C7) 
y dd 
+ Oe® 5 dé—u+u, (+x*x>& —x<6), 
oe 4+ 8 (x—ép 
0 0 
: : (7) 


1 y(&) (x —&) dé 1 | y(E) (« —&) dé 


e— 3x.) @—S'+ +h) 2a | @—F +O —h 


— —a 


oan g Rai g 
Sr 0 [r@ cos £ @ —2) Ge oa “0 [ro \+ Ae dg 
es a?) E g re 
Serko v5) } Sits FOE ae) ae Ean 8) 
Pe) 
& o dé £ 
a <3 (x —é)e = 2 dé =v, (taSé, -—x~< 6). 
3 Panera Sarmac) 
0 


3. Berechnung der Zirkulationsverteilung der Tragfliche. Genau wie bei einem gewohnlichen 
Tragfliigel im unbegrenzten Medium ist die Zirkulationsverteilung y(¢) aus der Randbedingung am 
Fliigel zu bestimmen. Diese lautet mit f(x) als Kontur der Profilskelettlinie 


u(x, —h)-f’(x) =v(x%,—h) (—aSx Sa). (8) 


Dies ist eine Integralgleichung fiir y(é), wenn man u und v gemaf (7) einsetzt; ihr Kern soll nun 
genauer untersucht werden. Dazu ist zunachst zu bedenken, daf fiir die meisten praktisch wichtigen 
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Falle 2h g/u2 < 0,7 sein wird; dann gilt aber in sehr guter Naherung 


H H 
He? dd b 1 1 he Gear X*\, H?+ xX? 
[aren fara! bo + peg ota ea — 
0 0 
xX? H x? HH? HB 
—X(1—"5 Jaret tS ee aa? (9) 
Xda f X 1 1 x X*\, H2+X? 
Sexe ‘ tas 9 wees Hake 
escent ct Bena” +5 8)d0 =F (1 =) ta X? 
0 0 
x? H > H 


Dabei ist zur Abkiirzung 


2 ee ard £ (x —é) =X 


gesetzt. In dem Integral (10) ist fiir X = 0 auch auf der rechten Seite der Wert Null zu nehmen, 
und zwar ist das Integral an der Stelle X = 0 unstetig, und der Wert Null ist das arithmetische | 
Mittel der beiden Grenziibergange X > + 0 und X > —0. 

Ferner benutzen wir die bekannte! Darstellung fiir den Integralcosinus 


CX) = 0,577 + + InXt— XP. (11) 


Damit nimmt unsere Integralgleichung (8) die Form an 


4 


2 uy f(z) == [re bee +K@ola (—asf<a), (12) 


und K(x, &) ist eine im ganzen Bereich in beiden Variablen stetige Funktion 
2hg 
(x —6)—2hf'(x aa . BN Ane 
K(x 6) = GSO IMO) 4 obo ob f(a) fwsin Ee —2) 4 singe —8) + 


+ ot red (a —&) + Si(o*7q sim ge 8) 


up 3 
ab hes 2 i 4h? + (x —é)? g 2 2h 
+3[ Ha s))In cp ag 9) Fae 24 ae ee 
2hg 2 h g\2 1 /2hg\3 
See ee) 
eae ae eh ea 
§ u2 = 
tan G fr cos E(w —8) £9 cosas (2 8) +5; (g ua Joos (x —€) 
=H F\ = & a 4h? —&) 
~ G(s 7 join (x Deere eet) eet ricer qf ee 
eae ees h gh eg h 
: 2s) arc te gp OL + Salt > 8 —2 <8). (13) 
Speziell fir € = x hat man 
1 ony Qhe | 2hg , (hg, 1/2hg\ = 
te, ena of Sy pe Py, 8g 8 8g ue 
K@,2)=—apf) +256 ples +m Ake 4 ABs (Re) 4 PA) og ee a. 


(14) 


Die Integralgleichung (12) mit dem Kernanteil (13) ist genau vom gleichen Typ wie die in der 
Schaufelgitterstrémung auftretenden Integralgleichungen, fiir die der Verfasser in friiheren Arbeiten 
eine vollstandige Auflésungstheorie vorgelegt hat® 3. 


1 Vel. z. B. Lésch-Schoblik, Die Fakultat, S. 122. Leipzig 1951. 
> W. H. Isay, Z. angew. Math. Mech. 33 (1953) S. 404. 
° W. H. Isay, Z, angew. Math. Mech. 37 (1957) S. 322. 
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Wir werden uns daher an dieser Stelle mit einer kurzen Darstellung der Lésungsmethode be- 
_ gniigen und verweisen fiir genauere Einzelheiten auf die genannten Arbeiten® 3, 


Wir setzen 


x=—acost, €=—acost (0s i Sa (15) 
T 
und ferner 
, 1 1 
— 2 f'@) == Ct), y@) == oft). (16) 


AuBerdem approximieren wir den stetigen Kernanteil K im Bereich 0 < i <a durch ein sechs- 
T 


gliedriges Fourier-Polynom in der Form 


5G 
asntK(t,t)= 2 2 b,,sinwicosyt. (17) 


w= 17=0 


Dabei sind die b,,, mit Hilfe einer harmonischen Analyse aus folgender Formel zu berechnen: 


5 5 7 ~ 
by = > sings f > sin tj K(Qj, 4) 008 » t+ SN OE 1K, 
a = (18) 
(w=, 55 eat = he), 
Fir y = 0 und y = 6 ist die Halfte des sich aus Formel (18) ergebenden Wertes zu nehmen. 
Nunmehr nimmt unsere Integralgleichung (12) die Form an 
1 ; sin t 2 < : 
G(t) = = Joo oe ee se ias sin st cos vt) dt (19) 
0 


Die bekannte Funktion G(t) kann unschwer (eventuell mit Hilfe einer harmonischen Analyse) in 
die Form 


60) = ne Sy sinew (20) 


= 1) 


gebracht werden; fiir die gesuchte Lésung w(t) machen wir den Ansatz 


w 2s 
one) = 28 + YZ 2 pcos ir. (21) 
Hiermit geht die Integralgleichung (19) in das lineare Gleichungssystem fiir die Koeffizienten w, 


1 Ie 
a = Cu — FPao Mo FO Pap MB (Gree Vie) (22) 


iiber, welches durch Iteration sehr leicht auflésbar ist, denn, wie sich an Beispielen gezeigt hat, 
sind die b,, im allgemeinen klein gegen eins. Die in den w (8 =1,...,5) noch enthaltene freie 
Konstante w, wird aus der AbfluBbedingung an der Profilhinterkante 


y(a) =0, bzw. w(x) = 0 


ermittelt. Damit ist die gesuchte Zirkulationsverteilung y(é) des Unterwassertragfliigels in der 
Form (16), (21) endgiiltig bestimmt. Die Gesamtzirkulation ist 


I =f y(é) dé = yaw. (23) 


4. Zwei hintereinander fahrende Unterwassertragflichen. Wir sind nunmehr ohne weiteres in der 
Lage, auch den fiir Tragfliigelboote wichtigen Fall von zwei im Abstand L hintereinander fahrenden 
Unterwassertragflichen zu behandeln (Abb. 2). Dabei setzen wir voraus, daB die EKintauchtiefe 
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der beiden Fliigel wahrend der Fahrt konstant bleibt, so daB es sich um ein stationéres Str6mungs- 
Fir unsere Untersuchung kénnen wir die beiden Tragflachen dann als ruhend 

in der homogenen Anstrémung u, betrachten. Die 
y Tragflache 1 liege bei —a, <& Sa, in der Tiefe hy 


problem handelt. 


unter der Wasseroberflache y ~ 0; die Tragflache 2 | 


liege bei —a,+L<&<a,+L in der Tiefe h,. 


“0, | &, 
att schwindigkeitsfeld u — ug = u, v = v, in der Form (7) 


iibernommen werden, und das Gesamtfeld beider 


Tragflachen wird durch Addition erhalten. Die | 


Abb. 2. Randbedingung an den Tragflachen 


[ug + Uy, (% 1, — My) + Uy, (%, — hy)] filx) = Vy, (x1, — hy) + vy, (41, — hy). 
[Uy + ty, (%2, — hy) + Uy, (2%, — he)] f2(%2) = vy, (%22 — hg) + vy, (2, — he) 


fihrt auf ein Integralgleichungssystem fiir die Zirkulationsverteilungen y,(é,) und y2(,); dieses 


lautet 
Le 1 ar 
— 2 Ug fil%y) = me | yi(é1) ee + Ky (%, 2) dé, + . | Ve(Eo) Kyo(%1, 2) 4. , | 
ay —a,+L | 
ay, a,+L 

, Ik af Sic y 
—2 uy file) = | rlG) Kalen bas +2 | rudd [ede + Kulensa| dias | 9) 

—a, —a,+L is | 

(-asSa, —_+Ls?satl), J 


Alle vier K-Funktionen sind in beiden Variablen stetig. K,, und K,, haben genau die Form (13). 
K,, und K,, kénnen aus (7) leicht entnommen werden; es ist nur 2h durch h, + h, zu ersetzen; 
auBerdem kénnen wegen (g L/uj)? > 0? an Stelle der Integralformeln (9), (10) nunmehr die ein- 
facheren Relationen 


Be? dd eeoyie a 8 
we § ¢ dd jel = || 
Neane af) an (a! cham irene ae en) z ~~ * (25) 
uj 2 Y) ua (%, a) Ceo = ua (x, ey uz (x,— &) 
0 


0 


verwendet werden. Damit wird 


K. ’ ao ae Sees (hy + he) fo(%e) X29 61 + (hy — hy) fo(%2) 
21(Xp» §1) (x, —&)? + (hy + h,)? si (x,—&)? + (hy aa hy)? 


Xo a \_- 1 (h 2 = he 
-} s,(e a sin § (x. —&,) 4 mth p24 C el ) Sie aaa é) 
0 


2 (x, — ey 


uz (an al &,) 


Analog ergibt sich K. . 4 § 
7/2 a ritice Vgl. ee dabei ist 2u beachten, daB wegen x, < é an Stelle von 32/2 nur 


Die Auflésung des Integralgleichungssystems (24) 


Gleichung id aatohei sate eae ge ee olgt genau in gleicher Weise wie die der 


arstellung des Verfassers verwiesen werden 


Ly wa 
aera j Fiir jede der beiden Tragflachen kann das Ge- | 
Se 
ay 
20, 


a 
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. 


es 


kann’. Durch eine zu (15), (16), (17) analoge Transformation geht (24) ber in 


1 j 5 6 
G,(t,) = — | ox(t)( oes + & 2 60:4 sin wt, cos y 2] dt, 


; cos T; — cos ty pee ee 


if is 5 6 
= al (T) 2 2 61?) sin wt, cos y T, dr,(0 = Sn), 
4=1ly=0 


F L Ty 
a 
I ip . ; 
Ca(te) = — a(t) & 2 be) sin wt, co8 vt, dt, si) 
F B=1 v=0 
era 
Lee sint ee 5 
ear | (o(T) be 7 Sor t & 2 be.) sin py t, cos » n dt, 
é 2 2 B=) y=0 
0<%<zx), 
yg 


und an Stelle von (22) haben wir jetzt die Gleichungssysteme 


1 
iN) : : 2 ; ; 
lee CC y (b0) co) a3 b(1.2) wo )— (60 1) @(1) + 0,2) w)), 


1 5 
) ie 
yen 
(2) (2) 1 (2,1 1 2,2 2 1 . 2,1 1 2,2 2 ce) 
co?) = C0) — (V2) al!) + 1G) oD) — FS (HED cl + GP wH), 


y2 
zur Bestimmung der Koeffizienten w®) bzw. w®) (a = 1,...5) der beiden in der Form (21) an- 


gesetzten Lésungsfunktionen @,(t,) und @,(t,). Die Konstanten @{) und w?) werden durch die 
beiden AbfluBbedingungen w,(z) = 0 und w,(z) = 0 festgelegt. 


5. Berechnung der Fliigelkrafte. Die resultierenden Fliigelkrafte K, und K, in x- und y- Richtung 
berechnen wir nach dem Kutta- Joukowskischen Satz: Die resultierende Kraft steht senkrecht auf der 


resultierenden Anstrémung {v gegen den Fliigel und hat die GréBe /K2 + K} =owT’, wobei I’ die Ge- 
samtzirkulation des Fliigels, 9 die Wasserdichte und tv die Geschwindigkeit ist, die am Ort des Fligels 
relativ zu diesem herrscht, abziiglich der von den gebundenen Wirbeln des Fliigels selbst induzierten 
Geschwindigkeit. Alle durch den EinfluB der freien Oberflache bedingten Geschwindigkeits-Anteile 
sind natiirlich in ty enthalten. Fiir diese Kraftberechnung denken wir uns die Fliigel durch in ihren 
Druck-Mittelpunkten (d. h. etwa in 1/, der Profiltiefe) angeordnete Punktwirbel gleicher Gesamt- 
zirkulation ersetzt. 
Somit sind fiir einen Einzelfliigel die Krafte gegeben durch 


K, = [up + uf) (0) + uM (OL, K, = —e [of (0) + PMO) L. (29) 
Dabei ist nach (7) 
2hg 
u2 
Ses y de 
If mm IF ® 8 
ay ee gat ta ° C+ fe oe | 
0 
2hg ? ; (30) 
if ihe Leah 2hg , 2hg ey 
ee iiecsn es [0,577 + In eee aye le re 
2hg 
g pine 
vf") (0) + v(t) (0) = — ae 0 


Die Kraftkomponente K, ist positiv, hat also den Charakter eines Widerstandes; dieser Widerstand, 
der ausschlieBlich durch den Einflu8B der freien Wasseroberflache bedingt ist, wird als Wellenwider- 


stand bezeichnet. 
Bei zwei im Abstand L > 0 hintereinander fahrenden Fligeln erhalt man die Krafte am vor- 


deren Fliigel J aus 


K0 = 0 [uy + ul? (0) + wl 0) + wy, (— DIT (31) 


KM = —Q [wv (0) + v8 (0) + »,,(—L)] Th, 
1 Siehe FuRBnote 3 von Seite 298. 
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und am hinteren Fligel 2 aus 


KQ) = Q [wy + uli? (0) + uh? (0) + uw, (L)] Ta» ! 
KP) = — 9 [of.) (0) + v2) (0) +, (DI Lh. 


(32) 


Dabei sind ul) (0) + ull (0) und v0) f voD(0) durch (30) gegeben, wahrend man aus (7) 


entnimmt 


— =e (hy + he) 
pad oth)t 4 4 MeL EMe vie ty ge fa 
ACO = eopea cee (seu Ae tiage ta) EL au ae 
ae. Gib (gL Gil ag 1 hy 
x Ee sin rp 4 C, ( Fe e°8 ua S. (Tg) sim 2 + x (=p } : 
eee hy + hy 
(ae ee Ll Gh). 4) 1 Ss g \ ) e 
ee | Day thy?  talq= hee ET ; 
ge 4g IL /g L gl Feel EN a. al. : 2 
A\24 5 sin a FeCl ges| Gs | S(Fe) sin re 5 a a +- hy)?}, “a 
g 33 
16 — a 1 Ee 1 : ews ee oye uk Soe 
OE) = oa ak WPL 2G be oe ee 
8 
=> (hy + he) 
32 gL gL gL ID) al: Sal 
x | cos a + §,(£2) cos ms uw sin eer e 
g 
il ag re) 
Dy (EL) saree i ee ee x 


2a (hy +h? +L? ' 2Qa(h,—h,)?+ 22 | 2 we 


0 


Seow Lig wed GIN hes hue 
a ear | s,( “| eos a c,(&) sin Teed e"o —] 


Uo 0 0 0 


6. Die Form der Wasseroberflache. Wir haben in Gleichung (1) eine Relation, die es uns ermég- 
licht, die Form Y(x) der freien Wasseroberflache zu berechnen. Der in (1) einzusetzende Wert u, 
ist durch (7) gegeben; und zwar mu8 u,(x,0) genommen werden, denn der ganzen bisherigen 
Theorie liegt die Annahme zugrunde, da die Wasseroberflache annahernd mit der Ebene y = 0 
zusammenfalle, und das Geschwindigkeitsfeld (7) erfillt gerade fiir y = 0 die kinematische Be- 
dingung (2) der freien Oberflache. Man hat somit 


Y(x) =— = u,(x, 0) . (1) 


Nun ist zu bedenken, da die Zirkulation y(€) unter der Voraussetzung bestimmt wurde, da die 
freie Oberfliche jedenfalls naherungsweise mit der Ebene y = 0 zusammenfalle, oder anders aus- 
gedriickt, daB die Oberflachenbedingung (2) im Abstand h von der Tragflache erfiillt wurde. Weicht 
nun die aus (1) ermittelte Form Y der Wasseroberflache in der Umgebung des Fliigels erheblich 
von der vorausgesetzten Form Y ~ 0 ab, so wird es notwendig, die Zirkulation y(&) zu verbessern 
oder sogar neu zu bestimmen. Hierfiir ware es natiirlich sehr wiinschenswert, eine Theorie zur 
Verfiigung zu haben, bei der die Bedingung (2) der freien Oberflache langs einer beliebigen Kurve 
Y(x) erfiillt werden kénnte. Aber leider stéBt eine solche Theorie auf sehr gro8e mathematische 
Schwierigkeiten, und es ist bisher nicht gelungen, hier etwas brauchbares zu entwickeln. Wir 
miissen das genannte Problem also mit Hilfe der in den vorangehenden Abschnitten dargestellten 


Theorie lésen. Dabei wollen wir die Falle Einzelfliigel und zwei hintereinander fahrende Fliigel 
getrennt behandeln. 


a) Einzelfliigel: Es hat sich gezeigt (vgl. die spateren Beispiele), daB das aus (1) berechnete 
und tiber die Umgebung des Fliigels gemittelte Y(x) nicht wesentlich von dem vorausgesetzten 
Wert (& Null) abweicht. Somit werden in diesem Fall keine zusatzlichen Uberlegungen notwendig.2 

Wir schlieBen hier noch einige Bemerkungen dariiber an, wie die praktische Berechnung von 
u,(x, 0) am besten durchgefiihrt wird. Wenn die halbe Fligellange a nicht wesentlich gréBer ist 
als die Eintauchtiefe h, so kann man den Fligel ersetzen durch einen in seinem Druckmittelpunkt 
angeordneten Punktwirbel gleicher Gesamtzirkulation J’. Bei vielen Profilen liegt der Druck- 


1 Sollte der iiber die Umgebung des Fliigels gemittelte Wert von Y(x) doch einmal wesentlich von Null 
abweichen, so kann die unter b) angegebene Methode auch hier verwendet werden. 
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mittelpunkt bekanntlich in der Nahe des sog. 1/,-Punktes, d.h. bei € ~ —a/2. Damit ergibt 
sich folgende bequeme Naherungsformel: 


hg 
Peers 0 
F g ap : es 0 wari 
u,(x, 0) See gic a gt sin X + 5 sin |X| + C,(X) cos X+ S,(X) sin X + | pee I’. (34) 


0 
q oe ag/x . é 
Dabei ist X = a e + 0,5}. Das Integral in (34) ist entweder mit Hilfe der Formel (9) oder mit 
0 : 
(25) auszuwerten, je nachdem ob §? ~ X? oder 3% <X? ist. 
Wenn andererseits die Fliigellange gro® ist gegen die Eintauchtiefe, oder wenn besondere 
_ Genauigkeitsanspriiche vorliegen, so hat man zunichst nach (7) 


+ S:(a* a") sin fe —2) 4 oe la + (96) 0,8) 48. — (5) 
0 pe + e (x > &)? ae 
0 : —a 


Soll u,(x, 0) etwa im Bereich —3a <x <7a _ berechnet werden, so wird die stetige Funktion 
O(x,&) mit Hilfe der Substitution 


*x=2a—5acost; & =—acost O22 2a) (36) 


durch ein Fourier-Polynom approximiert: 


6 6 
GOT) = 2 ee cost cosy 7, (37) 


lL 


o 


p=0 v= 


Die c,, werden mit Hilfe einer harmonischen Analyse berechnet: 


o 


 OF;, T,) COS 6 t; COS YT, 


ee 

ss 
rol 8 
~ 


k 


O(0, 0) + (— 1)’ 00, x) + (— 1)" O(a, 0) + (— 1)’ (— 1)" O@, 2)] 


ale LiMo 
a 


cos wt; [O(t;, 0) + (—1)’ O(;, 2)] (38) 


see da 
Sle 
Mea IMs 


cos » T%, [O(0, t;,) + (— 1)* O(a, t,)] 


| 8 
a 
\ 
ih 


Ae I 
Hd. 5: DM ez Il, as 6 DP bp jae =k). 


Hat entweder yu oder v den Wert null oder sechs, so ist die Halfte des sich aus Formel (38) ergeben- 
den Wertes zu nehmen; haben sowohl y als auch vy den Wert null oder sechs, so ist ein Viertel des 
sich aus (38) ergebenden Wertes zu nehmen. 

Ferner wird die Zirkulationsverteilung in der Form w(t) gemaB Gleichung (21) angesetzt. 
Damit ergibt sich die Geschwindigkeit w,(x, 0) zu 


6 1 6 

u,(t) = Te Dalene @y + pe Cup or) cosut, (x =2a—S5acost). (39) 
b) Zwei hintereinander fahrende Fligel: In diesem Fall fahrt der hintere Fliigel in dem 
Wellensystem des vorderen Fliigels; es hat sich gezeigt, dafS die durch den vorderen F'liigel in der 
Umgebung des hinteren Fliigels hervorgerufene Deformation der freien Oberfliche bei der Be- 
rechnung der Zirkulationsverteilung beriicksichtigt werden mu8. Umgekehrt beeinfluft der hintere 
Fliigel den vorderen fast gar nicht. Berechnet man nach der Theorie aus Ziff. 4 die Zirkulationen 
y, und y, der beiden Tragfliigel, so wird fiir diese Rechnung die Randbedingung (2) der freien Ober- 

flache im Abstand h, vom Fliigel 1 und im Abstand h, vom Fligel 2 erfiillt. 

21* 
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Wird nun aus Gleichung (1’), d. h. hier? 
¥ (x) = — Puy (x, 0) =P ty, (% 0) (1) 


die Form der freien Oberflache bestimmt, so wird sich im allgemeinen herausstellen, daB in der 
Umgebung des Fliigels 2 sein mittlerer Abstand von der freien Oberflache nicht h, sondern hf ist; 


1 2. He 2 3 
dabei kann h* wesentlich von h, verschieden sein; z. B [ths Zhe Die Zirkulation des Fliigels 2 


ist also unter einer nicht zutreffenden Voraussetzung bestimmt. Da keine Theorie vorliegt, die es 
erméglicht, die Oberflachenbedingung (2) langs einer beliebigen Kurve Y(x) zu erfiillen, bleiben 
zwei Méglichkeiten fiir die weitere Behandlung des Problemes. 

1. Man wiederholt die Berechnung der Zirkulationsverteilungen y,; und y, nach der Theorie 
von Ziff. 4; dabei wird lediglich h, durch h¥ ersetzt; denn da durch den ,,FremdeinfluB“ des Fli- 
gels 1 die Wasseroberflache den Abstand hy (und 
nicht h,) vom Fligel 2 hat, so muf fir letzteren 
auch die kinematische Oberflachenbedingung (2) 
im Abstand h* erfiillt werden. Mit den so neu 
berechneten Zirkulationsverteilungen y* und y7 
kénnen die an der freien Oberflache wirksamen 
Stérgeschwindigkeiten u*(x,0) und uj,(x,0) nach 
Formel (34) berechnet werden. (Fast immer wird 
ubrigens y, = yj sein.) 

Dann liefert 


¥*@) == = ux (x, 0) — - ux (x, 0) 


die verbesserte Form der freien Oberflache; diese 
ist wegen 


Abb. 3. a) 1. Méglichkeit; b) 2. Méglichkeit. u, (— co, 0) = 05 a {— co, = 0 


von dem alten Nullniveau (weit vor Fliigel 1) aus zu rechnen. So kann dieses Verfahren iterativ 
fortgesetzt werden. In vielen Fallen ist bereits Y*(x) ~ Y(x), so daB nur ein Iterationsschritt er- 
forderlich ist (vgl. Abb. 3). 

2. Die andere Méglichkeit ist folgende (ohne Iterationsverfahren). Man behalt die im ersten 
Rechengang erhaltenen Zirkulationsverteilungen y, und y, bei und verandert die Tiefenlage 
des Fliigels 2 gegeniiber dem Fliigel 1 entsprechend der Gestalt der freien Oberflache. Man hat 
damit jedoch nur ein Problem, das gegeniiber dem urspriinglich vorgegebenen abgedndert ist, 
gelist (Abb. 3). 


7. Die Y44-—3/4-Punkt-Methode. In vielen Fallen leistet in der Tragfliigeltheorie die sog. 14—3/,- 


Punkt-Methode gute Dienste, um naherungsweise die Gesamtzirkulation J" eines Fliigels zu be- | 


stimmen. Die Methode besteht bekanntlich darin, da8 man den Tragfliigel durch einen in seinem 
1/,-Punkt € = — a/2 angeordneten Punktwirbel ersetzt und die Randbedingung am Fliigel im 


3/4-Punkt x/a = 1/2 erfiillt. Hierdurch ergibt sich eine lineare Gleichung zur Berechnung von J’, | 
und man umgeht die Auflésung der Integralgleichung; natiirlich erhalt man dabei nicht die Zir- _ 


kulationsverteilung y(&). Auch fiir zwei hintereinander fahrende Fliigel 1a8t sich die Methode 
anwenden und liefert zwei Bestimmungsgleichungen fiir J, und J;,. 


Durch Vergleich mit Beispielen, die nach der genauen Theorie von Ziff. 3 und Ziff. 4 durch- | 


gerechnet wurden, hat sich gezeigt, da die Zirkulation J’ auch bei Unterwassertragfliigeln in recht 


guter Naherung nach der 1/,—*/,-Punkt-Methode bestimmt werden kann. Dieses ist z. B. wichtig | 
fiir das in Ziff. 6b) beschriebene Iterationsverfahren, da man fiir die erste Naherung gut diese | 
Methode verwenden und sich somit einmal die Auflésung des Integralgleichungssystems (24) er- 


sparen kann. 


Wir bezeichnen mit u(x, y), vp(x, y) das aus (7) sofort zu entnehmende Geschwindigkeitsfeld 
eines bei § = — a/2 gelegenen Punktwirbels der Zirkulation J’. Die Randbedingung im 3/,-Punkt 


lautet 
wl (B) Hee AS) 8) 


/ 


1 Bei der Anwendung von Formel (34) auf Fliigel 2 ist natiirlich X = 8 lx + a a,—L), Analoges gilt 
: a2 oy C4 . 
fiir Formel (39). 0 


Rey RT 


oe 
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_ und daraus ergibt sich nach einer elementaren Rechnung die Bestimmungsgleichung fiir J” 


/ ie Jp! A H 
Spe Creare a ec + [3 +2 8,4) +2" C(A)] cos A 


+ [3af’ —2 C(A) +2f" S,(A)] sinA—2(1—F + Af — 7 Af’) are ga 


A 
eet Ate 
—AH(L+5)+HP (2 eer ps (40) 


Dabei ist H = 2hg/uz und A = (ag/u2). /2/H und /1/A kénnen als Froudesche Tiefenzahl 
bzw. Langenzahl bezeichnet werden. 
Dealer hat man bei zwei Fliigeln die Randbedingungen in den */,-Punkten 


wean.) rea) =or( 2.8) nf) 
luo + ur (ee he) + ur, (L +32, —hs) fa(L + |=, i +3, ha) + on,(E +3, — he), 


und diese fiihren auf die folgenden Gleichungen zur Berechnung von J), und J): 


= woh!) = aaa + apt A en fae SFr + [Ba +2 SA) +26 C(A)] 005 4, | 
+ [Barf —2 CfA) + 2f S(A)] sin 4, —2(1 FE + AK 5 ALS aro tg 
[abate bantam +) + mh +33 Ail 
ae Fn, Ae a paean ai ae + [ +2 Spe) + 2f{ C:(Ays)] 008 Ayo 
els, — C(Ay2) + 2 fy S(Aye)] sin Ar, Fh 
wo fell + 2) = a5, ys a joe eta 4 2 | (41) 
+ [Bm +2 S(dy) + 2f; GlAy)] 008 dy + Bf, —2 Cn) + 2f3 5(4q)] sin Ay 
eS me eae + Ba +2 84 2) + 2f3 C(A,)] cos A, 
+ [32 f,—2 C(A,) +2 f% S,(A,)] sin A, — 2 ( Sa 5 Ais) arc tg a 
(44.5 +7 A fi)n BES aml + B+ Be +B 7A 
Dabei ist 
H, =e, Hy, ae, Hy = Hy ="5g, 4, =98, 4, = 5%; 
Ay, =£(a—L), 4, =2£(0+]), ie 


In Anbetracht des Naherungscharakters der 1/,—?/,-Punkt Methode sind in den Gleichungen (41) 
einige fiir die numerische Rechnung ganz unbedeutende Glieder weggelassen worden. 


8. Beweis der Integralformeln (5). Die Integralformeln (5) wurden in etwas anderer Form zuerst 
von Krienes angegeben. Da die Krienessche Untersuchung jedoch nicht veréffentlicht und nur 
schwer erhaltlich ist}, erscheint es gerechtfertigt, den Beweis der Integralformeln (5) hier zu 


1 An dieser Stelle méchte ich Fraulein E. A. Walinski von der Versuchsanstalt fiir Wasserbau und Schiff- 
bau, Berlin, dafiir danken, da sie mir liebenswiirdigerweise erméglichte, die beiden genannten Arbeiten von 
Krienes einzusehen. 
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bringen. Zudem ist der nachfolgende Beweis verschieden von dem von Krienes angegebenen; 
wahrend wir ganz im Reellen bleiben werden, benutzt Krienes komplexe Integrationswege. 


Zunachst ist 
oe) 


es cos u X du = e—# liars cos X (y + 1) dy, 


co 


—1 
— ap 
‘fe we sin X du ali sin X (y + 1) dy. 
il 


Dann wird 
oe (oo) eee, 0 co ae 
iE iy dyes =e At sin X yd y | [57 i SI <— sin X ydy + S(X) £5. 
pe ‘ 
Sa =i =X | 
Dabei gilt das + Zeichen fir X > 0 und das — Zeichen fiir X < 0. Fiir X = 0 hat das Integral 
den Wert Null. Ferner ist 


Be 
S(X) = | ee (— S(— X) = $20) 
0 
der Integralsinus. Analog haben wir 
lee cos Xydy = =i a cos X y dy +{ do = =f 1 cos X x dy — C{X), 
=f t —l —x —1 i 
fiir beliebiges X, dat 
“ | 
| coso— = 0 
=< 


ist; dabei ist 
Cros | cosa (C(— X) = C(X)) 
; x 


der Integralcosinus. Unter Verwendung dieser Formeln erhalten wir 
oe) 


“ew H _a 
| Fe cos wu X du = —e 


C,(X) cos X + $,(X) sin X 2 sin X 4 [=o exe tae) 


: = (42) 


eer : _A : Tt ree ; 
i eae oe AX du = +e si cos X — C,(X) sin X -- > cos o.¢ =a Payee oe +1) ir 


ae (43) 


Nun sind aus der Integralrechnung die Formeln 


= Be? ieee 3 Xe? 
Je POO8 DAY tr 1) OX comes Je st a) Oe rtare 
a —1l 


1 Die Integralformel gilt auch fur X = 0; denn man erhalt unter Benutzung der Relationen (9), (11), (44) 
leicht das Ergebnis 


1 1 
0,977 — ln H — H — 2 3... 
n a H 18 H 
Das gleiche ergibt auch ein direktes Ausrechnen des Integrals 
ioe) H 
: dz “l—e™% 
e—H x — dy— Ei(— H 
Wace ace, mise 
—l —H 


unter Verwendung der entsprechenden Reihenentwicklungen. 
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bekannt; aus ihnen folgt durch Integration tiber @ von 0 bis H 
oO A 
1 — e—Hy Be? 
a xX 1 = | = 
[ = cos X (y + 1) dy Lees 
nis] 0 
a (44) 


pees aia (71 I dy dd 
| ial pag) Peers : 
soy 0 


| 

if 

H ( 
Xe? | 


Durch Einsetzen dieser Relationen in (42), (43) ergibt sich (5). 


9. Zahlenbeispiele von Einzelfliigeln.!. Bei allen drei folgenden Beispielen handelt es sich um 
Einzelfliigel mit dem Plattenprofil f’(x) = konst. = 6 = — 0,1. 


Dariber hinaus sei fiir Beispiel 1 ag/u2 = 0,1, ging 0, Le, 
fur, Beispiel 2: ag/u7\== 0,1, kh glu? = 0,05, 
far Beispiel 3. ca s/uz = 0,05 bh g/u2 = 0,05, 


up//h g = 3,16; baw. 4,47. 


In den Tabellen 1 bis 3 sind die Fourier-Koeffizienten b,,, des stetigen Kernanteils der Integral- 
gleichung (19) enthalten. 


Das entspricht Froudeschen Tiefenzahlen von 


Tabelle 1. (Beispiel 1) buy 


ae 0 1 2 3 4 5 6 
ia 
1 +0,567 +0,284 —0,006 —0,009 0,000 0,000 0,000 
2 —0,131 +0,013 | =-+0,030 —0,001 —0,002 0,000 0,000 
3 —0,003 —0,028 +0,001 +0,004 0,000 0,000 0,000 
4 +0,003 —0,001 —0,004 0,000 +0,001 0,000 0,000 
) 0,000 +0,002 0,000 —0,001 0,000 0,000 0,000 
Tabelle 2. (Beispiel 2) b,, 
| 
> 0 al | 2 2 4 5 6 
bi a 
1 +0,659 +0,633 —0,019 —0,044 +0,001 +0,002 0,000 
D, —0,258 +0,037 +0,161 —0,006 —0,021 +0,001 +0,001 
3 —0,005 —0,131 | +0,009 +0,051 —0,002 —0,008 0,000 
4 +0,001 —0,004 —0,042 +0,003 +0,016 —0,001 —0,003 
5 0,000 +0,008 —0,002 —0,016 +0,001 --0,006 0,000 
Tabelle 3. (Beispiel 3) b,,» 
| 0 1 2 | 3 4 5 6 
yu i] 
i +0,341 | +0,238 —0,005 —0,008 0,000 0,000 0,000 
2, 04107) +0,009 +0,027 —0,001 —0,002 0,000 0,000 
3 | —0,002 | —0,025 -+0,001 +0,004 0,000 0,000 0,000 
4 +0,003 | —0,001 —0,004 0,000 +0,001 0,000 0,000 
5 0,000 +0,001 0,000 —0,001 0,000 0,000 0,000 


In Tabelle 4 sind die unter Beriicksichtigung der AbfluBbedingung berechneten Koeffizienten wg, 


der Zirkulationsfunktion w(t) gemaB (21) in Einheiten von uy enthalten. Rechts neben der Tabelle 


stehen die Gesamtzirkulationen nach (23); ferner haben wir diese Gesamtzirkulation in Beziehung 
gesetzt zu der Zirkulation /’,, = 0,2 a uy des gleichen Fliigels im unendlich ausgedehnten Medium, 


1 Fir die Durchfiihrung der in dieser Arbeit enthaltenen numerischen Rechnungen danke ich Frau Marzke 
und Frau Mutzek. 
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d.h. fiir h->0o. Daraus erkennt man, wie stark die Zirkulation des Fliigels und damit auch die 
Auftriebskraft herabgesetzt wird durch den Einflu8 der freien Wasseroberflache. 


Tabelle 4 
: Ip 
Beispiel 1 | Beispiel 2 Beispiel 3 Beispiel 1: J’ = 0,343 aug, rT, = 0,546. 
if 
2 SGeL ceo es Benes Beispiel 2: I’ = 0,289 a ug, 7,~ = 0,460. 
o, 40,017 +.0,025 0,017 e 
o +-0,003 0,009 0,003 oe ra itn 
i 0,000 +0,001 0,000 Beispiel 3: I’ = 0,403 a uy, To = 0,642. 
Ws 0,000 0,000 0,000 
Aus den Werten der Tabelle 4 erhalt man mit (15), (16) die Zirkulationsverteilungen y(&) namlich 
ns a—é a as ie 
Belepioltee oe) =: = (0,096 0,018 = + 0,010 =), 
pe g ee ce 
Beispiel 2: y(g) =]/2=$ (0.073 — 0,011 | + 0,022 — 0,006 5), 
ispi fa=E (0,115 — 0,018 + 0,010 = 
Beispiel 3: y(&) = \ epre: > > y > a) 


Fiir die resultierenden Krafte erhalt man aus (29): 

K, = 0,332 9 aug , 
K, = 0,0097 9 a uz, 
ie = 0,272 0.015 | 
K, = 0,0076 0 au2, 
K, = 0,387 9 aug, 
K, = 0,00749 auz. 


Beispiel 1: 
Beispiel 2: 
Beispiel 3: 


In Tabelle 5 ist die Form a Y(x) der Wasseroberflache berechnet nach Gleichung (1), (34) enthalten 
(vgl. auch Abb. 4, 5, 6). 


Tabelle 5 
pal | = =) at 
ee | | 
Y Beisp. 1 +0,085 | +0,132 | +0,163 | +0,132 | +0,039 | —0,050 |} —0,200 | —0,325 | —0,508 
Beisp. 2 +0,077 | +0,129  -+0,187 | +0,160 | +0,047 | —0,039 | —0,176 | —0,287 | —0,449 
Beisp. 3 -+0,166 +0,236 | +0,286 | +0,267 | +0,179 | +0,100 | —0,027 | 0,132 —0,307 
; I Ba 
-Y4 y sh ee 


7 ee 
V(x) 


~05 
SS ‘ 
Abb. 4. Beispiel 1. Form der Wasseroberfliche Y(x) tiber dem Fliigel. 
43/2 
O5}- 
| I ! ! | La 
au D4 =f 7 a eee 
SS V(a) 
Abb. 5. Beispiel 2. Form der Wasseroberfliche Y(x) iiber dem Fligel. 
4y/a 
Qs} 
eT a Teel EMR ree ee 
¥¥ no, 7 7 a Y a Se 
-05"- V(z) 
sreencone oe SU 


Abb. 6. Beispiel 3. Form der Wasseroberflache Y(x) iiber dem Fligel. 


es 
i: 
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é 1 
Eine Berechnung von — Y(x) nach der genauen Methode mit Formel (39) wurde ebenfalls durch- 


gefiihrt; dabei zeigte es sich, daB Abweichungen von den Werten aus Tabelle 5 in der dritten Stelle 
hinter dem Komma auftraten. Da diese Stelle in Anbetracht der Naherungstheorie der Ober- 
flachenform sowieso keine echte Bedeutung mehr hat, kann das Ergebnis aus Tabelle 5 als befrie- 
digend angesehen werden. 


Eine Berechnung der Zirkulation [’ nach der Niherungsformel (40) ergibt fiir 
Permpiela |: 2b y— 0.350 Ag, 
Beispiel: 22> 1==:0,271:a a, , 
Beispiel 3: J7=0,395au,. 


Diese Werte stimmen schon relativ gut mit den aus der genauen Theorie durch Lésung der Integral- 
gleichung erhaltenen Werten iiberein. Die Abweichung ist natiirlich fiir Beispiel 2 am gréBten, da 
dort das Verhaltnis a/h am griéBten ist. 


Zusammenfassend kann man sagen, da®B durch den EinfluB der freien Wasseroberflache auf den 
Fliigel in erster Linie ein betrachtliches Absinken der Fligelzirkulation I" (gegeniiber dem Fall 
h—> oo) bewirkt wird. Dagegen ist die Stérgeschwindigkeit (30) am Fliigel und damit auch der 
Wellenwiderstand K, relativ unbedeutend. 

AbschlieBend sei noch auf MeBergebnisse von Schuster und Schwanecke! hingewiesen. Diese 
fanden fiir einen Fliigel endlicher Spannweite (Seitenverhaltnis 5,84) mit Kreissegmentprofil und 
dem Wert ag/u? = 0,0425 bei einem geometrischen Anstellwinkel von 5,7° (6 =—0,1) fir 
h/a = 0,5 den Wert I’'/I’,. = 0,58 und fiir h/a = 1 den Wert I'/I’,, = 0,69. I’ und I’,, (A> 00) 
sind dabei die iiber Spannweite gemittelten Werte. 

Die hier vorgelegte Theorie ergibt fiir ein Plattenprofil unendlicher Spannweite bei gleichen 
iibrigen Parameterwerten /'/J",, = 0,56 bzw. l/l’. = 0,66. 


10. Zahlenbeispiele von zwei hintereinander fahrenden Fliigeln. Auch fiir die beiden folgenden 
Beispiele 4 und 5 legen wir Plattenprofile mit f;(x;) = konst. = 0; (j = 1,2) zugrunde. AuBerdem 
sollen Fligel 1 und 2 jeweils in der gleichen Tiefe. unter der Uneestorten Wasseroberflache fahren 
und die lithe Fligellange haben: a, = a, = a (Abb. 2). 7atlenmaGie sei fiir 


ee Cg ee ADA 

Beispiel 4: Pigs th hah pa 6, = 6, =— 0,1, = 20a, 
Pa Lee age hye hye 

Beispiel 5: FU Ras eg ie EN ele aa L=2W e. 


Der gewahlte L-Wert entspricht einem bei Tragfliigelbooten iiblichen Abstand. 

Die Berechnung der Zirkulationen erfolgt mit Hilfe des in Ziff. 6b) erlauterten Iterationsver- 
fahrens. Fiir die Ausgangsrechnung haben wir die folgenden b,,,-Koeffizienten der stetigen Kern- 
anteile des Integralgleichungssystems (27): 

Beispiel 4: Fir b(t.) = b(?.?) die Werte aus Tabelle 1; 

b(2) — — 0,013, bf?) = — 0,002, 6G?) = + 0,001, 
b21) — — 0,501, 62) = — 0,091, bG!) = + 0,045, 6G) = — 0,003, 6G) = + 0,001; 


alle ttbrigen 6(!,?), b:) sind vernachlassigbar klein. 


Beispiel 5: Fir b(;1) — b?.?) die Werte aus Tabelle 3: 
bt?) — — 0,030, bt) = — 0,003 , bG) = + 0,001 , 
b@-) — + 0,297, 6G) = — 0,029, bG) = + 0,014; 
alle iibrigen 6(,”), 6(?:!) sind vernachlassigbar klein. 
In Tabelle 6 sind die aus (28) berechneten Koeffizienten w{) und @{)) der Zirkulationsfunk- 


tionen @,(z,) und w,(t,) in Einheiten von uy sowie die Gesamtzirkulationen J}, /, in Einheiten von 
a uy enthalten. 


1S. Schuster u. H. Schwanecke, Schiffstechnik 4 (1957) S. 117. 
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In Tabelle 7 ist die Form = Y(x) der Wasseroberflaiche berechnet nach Gleichung (1’), (34) 


enthalten. Um das Ergebnis tibersichtlicher zu haben, ist der Zahlenwert fiir = Y(x) dabei in zwei 


Summanden aufgespalten, die den Stérgeschwindigkeiten von Fliigel 1 baw. Fliigel 2 an der freien | 
Oberflache entsprechen. 


Tabelle 6 
Beispiel 4 
Te 0,347 
ow» +0,196 
ON? +0,153 
OW +0,017 
wo? +0,003 
wp 0,000 
NS) 0,000 | 
IRs, +0,458 | 
ae +0,258 
w\? +-0,196 
wo -+0,017 
OY | +0,003 
OP 0,000 
a? 0,000 


spricht. Aus den Werten der Tabelle 7 entnimmt man}, dai es angemessen sein diirfte, h¥ = =i 
zu setzen (fiir beide Beispiele). 


alle tibrigen bt!.), b 


Man erkennt, daB in der Umgebung des Fliigels 2 der mittlere senkrechte Abstand dieses Fligels 
von der Wasseroberfliche erheblich kleiner ist als fir die Zirkulationsberechnung vorausgesetzt 
wurde. Gema unseren Uberlegungen aus Ziff. 6b) wiederholen wir also die Berechnung der Zir- 
kulationsverteilungen mit einem neuen h¥-Wert, der den tatsdchlichen Verhaltnissen besser ent- 


¥/@ 


= =e _ we 
TE 5 

Via) a 
Abb. 8. Beispiel 5. Form der Wasseroberflache Y(x) iiber den Fligeln. } 
Fiir diese Iterationsrechnung haben wir folgende Koeffizienten 5,,,: | 
Beispiel 4: Fiir bth) die Werte aus Tabelle 1, fir b(?.?) die Werte aus Tabelle 2: | 
bG?) — — 0,014, b(t?) = — 0,002, bt?) = + 0,001, | | 
bt) —= — 0,526, b@)) = — 0,096, be = + 0,048 , bg) = — 0,003 , bG) = + 0,001; | 


Ll 
1 Denn es ist — 


12 


(2,1) 
by 


Beispiel 5 


+0,411 
40,232 
+0,178 
40,017 
+0,003 
0,000 
0,000 


+0,333 
40,188 
+0,142 
+0,011 
+0,002 
0,000 
0,000 


Abb. 7. Beispiel 4. Form der Wasseroberflaiche Y(x) iiber den Fligeln. 


Tabelle 


Y/a Beispiel 4 


0,100 = + 0,085 + 0,015 


—0;028"= — 0,051 = 0,023 
lee), 176) = == (0,203.42.0.029 
2-),297 = 0329 450,032 
—0,466 = — 0,513 + 0,047 
—0,523 = — 0,603 + 0,080 


0,150 = + 0,133 0,017 
10,184 = + 0,165 + 0,019 
SO:153 = = 0533.22 0.020 

0,039 + 0,021 
= 0.59020. 0195 


—0,353 = = 0,570 40,217 
—,369.==— 0,5442- 0,175 
25,462 = —0,51340,051 


— 0,477 — 0,067 
— 0,390 — 0,267 
— 0,172 — 0,569 
+ 0,041 — 0,751 


Y/a Beispiel 5 


il tt 


169 + 0,029 | 

241 +-0,032 
291 + 0,034 
272 + 0,036 

182 + 0,038 | 

101 + 0,040 
— 0,028 + 0,045 
— 0,134 + 0,052 
— 0,312 + 0,070 
— 0,459 + 0,105 
0,578 + 0,195 


ba all 


— 0,603 + 0,237 
— 0,627 +: 0,221 
— 0,648 + 0,148 | 
— 0,668 + 0,082 | 
— 0,701 — 0,022 mI 
— 0,745 —0,185 an 
~0;757 = 0816 


Pid weed 


26 


| = Yad) mw — 0,5. 


sind vernachlassigbar klein. 


/ 
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Beispiel 5: Fir b(.1) die Werte aus Tabelle 3, fiir b(?:?) die Werte aus Tabelle 8; 


bd2) = — 0,031, bG2) = — 0,003, 602) — + 0,001, 
be; 2) — 4.0.04, be, ash oa0 be) = +.0,015 ; 
alle tibrigen b(!.), 6:1) sind vernachlassigbar klein. 
Tabelle 8. b(?: 
Se 
Ny 0 1 2 3 4 5 6 
ue 
| | 
1 | -+0,386 +0,548 | —0,016 —0,041 +0,001 -+-0,002 0,000 
2 —0,219 +0,031 | +0,151 —0,006 —0,020 --0,001 +0,001 
3) —0,004 —0,123. | +0,009 -+-0,049 —0,002 —0,008 0,000 
4 +0,001 —0,004 | —0,040 +0,003 -+-0,016 —0,001 —0,003 
5 0,000 + 0,007 —0,002 —0,016 +0,001 -+0,006 0,000 


In den Tabellen 6* und 7* sind (entsprechend Tabelle 6 und 7) die Ergebnisse der Iterations- 
-rechnung enthalten. Die Form = Y(x) der Wasseroberflaiche (vgl. Abb. 7 u. 8) nach Tabelle 7* 


unterscheidet sich nur noch wenig von derjenigen aus Tabelle 7; somit eriibrigt sich ein weiterer 
Tterationsschritt, und die Ergebnisse aus Tabelle 6* und 7* kénnen als endgiiltig angesehen werden. 


Tabelle 6* 

Beispiel 4 | Beispiel 5 
in | -+0,347 | -++0,410 
eo 40,196 | -£0,232 
a 40,153 | +0,178 
oP +0,017 | +0,017 
es +0,003 | +0,003 
ow 0,000 0,000 
oo 0,000 0,000 
r, | +0,390 | +0,285 
wo | +0,220 | +0,161 
o® | +0,174 | +0,126 
OD icia0 020m 0,019 
oo +0,012 | +0,008 
ye? +0,001 | +0,001 
oy ® 0,000 0,000 


Tabelle 7* 


Y/a Beispiel 4 


Ho we we dl 


Hoi i wt ya 


+ 0,085 + 0,014 | 


=> 0,133 + 0,015 
+ 0,165 + 0,017 
= 05133 =- 05013 


— 0,051 + 0,021 
— 0,202 + 0,024 
— 0,328 + 0,029 
— 0,513 + 0,043 
— 0,602 + 0,073 
— 0,588 + 0,174 
— 0,569 + 0,253 
— 0,543 + 0,215 
— 0,513 + 0,063 


=- 0,039 + 0,019 | 


==) ATT —=0,053. | 
= 0,390 —0,237-1 


— 0,172 —6,510 
+ 0,041 — 0,674 


Damit ergeben sich folgende Zirkulationsverteilungen: 


Beispiel 4: 


AA) == maa (0, 097 — 0,018 = = + 0,010 


+S) 


a, ee E é 
Vo(Fo) = (eos 0,102 — 0,008 (fs a 20) + 0,032 é as 


ee ola 


Beispiel 5: 


ils) = = ASE (0,17 — 0,018 = + 0,010 ay 


+0,194 
+0,268 
+0,320 
+0,303 
+0,215 
+0,136 
+0,012 
—0,089 
—0,250 
—0,365 
—0,394 
0,353 
—0,390 
—0,506 
—0,590 
oa 
—0,905 
1,031 


Y/a Beispiel 5 


+ 0,169 + 0,025 


Pa WW eh a aS eT 


0,240 + 0,028 
0,291 + 0,029 
4+ 0,272 + 0,031 
+ 0,182 + 0,033 
+ 0,101 + 0,035 
— 0,028 + 0,040 
— 0,134 + 0,045 
—. 0,312 + 0,062 
— 0,458 + 0,093 
— 0,577 + 0,183 
— 0,602 -+ 0,249 
— 0,625 + 0,235 
— 0,647 + 0,141 
— 0,667 + 0,077 
— 0,700 — 0,018 
— 0,744 — 0,161 
20,155 0.276 


(—a Se a) 


\2 A 3 
20| — 0,006 — 20) 


2 
vals) = Vo 21 4 /0,077 — 0,005 (8 — 20) + 0,019 (£20) — 0,006  £—20)) 


Sea IY 


a 


[19 <fec 21). 
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Fiir die resultierenden Fligelkrafte erhalten wir aus (31), (32) (natiirlich mit hy): 


Beispiel 4: K® = — (0,965 — 0,002) 0,347 9 ug a = 0,334.9 ua, 
K) = + (0,028 — 0,001) 0,347 9 u2 a = 0,0095 9 ub a, 
K?) = (0,920 + 0,053) 0,390 @ ug a = 0,380 @ Us Os 
K®) = + (0,035 — 0,024) 0,390 0 u2 a = 0,0045 9 uga. 
Beispiel 5: K) = (0,958 — 0,002) 0,410 9 ua = 0,392 9 upa, 
K = + (0,019 — 0,002) 0,410 9 ua = 0,0069 9 ura, 
K?) = (0,944. + 0,030) 0,285 @ ug a = 0,277 @ Lege 
K®) = + (0,013 + 0,023) 0,285 9 u2 a = 0,0103 9 ug a. 


Der zweite Zahlenwert in den Klammern stellt jeweils den Anteil der vom anderen Fliigel indu- 
zierten Geschwindigkeit dar. 


Eine Berechnung der Zirkulationen J), und J’, nach den Naherungsgleichungen (41) der 1/,—*/,- 
Methode ergibt 


Beispiel 4: J) = 0,340 au, , T= 0,30314 Ue 
Beispiel 5: J =0,40lau,, I, = 0,266 aug. 


Die Abweichung von den nach der genaueren Theorie ermittelten Werten aus Tabelle 6* ist fiir J, 
gréBer als fiir J, denn Fligel 2 hat den gré®eren a/h-Wert. a 

Hine Betrachtung der Fliigelzirkulationen [,, [, aus Tabelle 6* zeigt: Fir Beispiel 4 und 5 
stimmt J’, fast genau mit dem /-Wert des entsprechenden Einzelfliigels (Beispiel 1 und 3) uberein; 
ya dagegen ist J, bei Beispiel 4 (trotz 
si Ol Ge Selene eke Ne eae aja des geringeren Oberflachenabstan- 
ae. a Via) a 5: des h* von Fliigel 2) etwas gréBer 
s ae als I, wahrend bei Beispiel 5 J, 

Abb. 9. Beispiel 6. Form der Wasseroberflache Y(x). a aa ertohehinter F, zuriickblesaes | 
ula Diese Unterschiede sind dadurch 
Sree eee ee nt pee a bedingt, daf} die Periode des hinter 
meh : ‘2 g Fliigel 1 laufenden Wellensystems bei 
Abb. 10. Beispiel 7. Form der Wasseroberflache Y(x). Beispiel 4, halb so grok ist wie bei 


yfa Beispiel 5, wahrend der Fliigelab- 
a l l 7/2 stand L in beiden Fallen gleich groB 
0 V(a2) caper oe ist. Durch das Wellensystem hinter 

ad rare Fligel 1 wird ja der am Fligel 2 
Abb. 11. Beispiel 8. Form der Wasseroberflaiche Y(x). wirksame Anstellwinkel stark beein- 


y/% fluBt. Entsprechend ist auch der 
pees es oa eae Wellenwiderstand K() bei Beispiel 5 


3 vi | ee Ga wesentlich gréBer als bei Beispiel 4. 
— (2) oe Sicherlich laBt sich durch geeig- 
Abb. 12. Beispiel 9. Form der Wasseroberflache ¥(2). nete Wahl der geometrischen An- 
stellwinkel, d. h. der Werte von 6, 
he a und 6,, erreichen, daB I’, = I, wird, 
- a + so daB auch die Auftriebskrafte beider _ 
— -- Fligel ungefahr gleich groB werden. 


Praktisch wichtig ist die Frage, 
inwieweit bei einem bereits durch- 
by /a gerechneten konkreten Fall (etwa_ 
ele . Zt unsere Beispiele 4 und 5) die Werte 
Bass, Va) s von 6, und 6, abgeadndert werden | 

y miissen, um /, = J, =J" 2u errei- | 

chen. Dieses Problem ]a&t sich mit 


Hee Hilfe der Gleichungen (41) lésen; 
2/2 denn gibt man fir [T=J,=T,_ 
10 Va) “0m | 
a 


a einen bestimmten gewiinschten Wert — 
Abb. 15. Beispiel 12, Form der Wasseroberflache Y(s). vor, so stellt (41) zwei lineare Glei- | 


Abb. 13. Beispiel 10. Form der Wasseroberfliche Y(x). 


Abb, 14. Beispiel 11. Form der Wasseroberflache Y(x). 


es 
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: chungen zur Berechnung von 6, und 6, dar. Dabei ist jedoch zu beachten, daB der vorgegebene /-Wert 
auch eine neue Form der Wasseroberflache bedingt, so daB unter Umstanden auch der Bicherire h,-bzw 
hy-Wert entsprechend modifiziert werden mu8. Wahlt man jedoch (was praktisch sinawolll SE 
scheint) das neue J” gleich dem bisherigen J’, oder als arithmetisches Mittel der bisherigen [\-und 
iy 2 Werte, so diirfte die Formyeranderung der Wasseroberflache nicht sehr erheblich sein ie daB 
man sie jedenfalls fiir Uberschlagsrechnungen vernachlassigen kann. : 


“a Wir wenden das eben geschilderte Verfahren auf die Beispiele 4 und 5 an: 
Beispiel 4: Soll l= 7) = 0340.4 Uy werden, so ergibt sich 6, = — 0,100, 6, = — 0,092; 
eine h}-Korrektur ist dabei nicht notig. 
Soll a= Ty a + I, = 0,352 a uy werden, so ergibt sich 6, = — 0,104, 
0. = — 0,095; eine h3-Korrektur ist dabei nicht notig. 


| Beispiel 5: Soll J’= J, =0,40lau,) werden, so ergibt sich ohne eine h¥-Korrektur: 
0; = — 0,099, 6, = — 0,143; dagegen erhalt man bei einer (angemessenen) Ver- 
gréBerung des h}-Wertes um 6%: 6, = — 0,099, 6, = — 0,141. 


Soll 7 = +1, +4), = 0,334 a uy werden, so ergibt sich ohne eine h*-Korrek- 
tur: 6; = — 0,082, 6, = — 0,119; dagegen erhalt man bei einer (angemessenen) 
VergréBerung des h*¥-Wertes um 20%: 6, = — 0,082, 6, = — 0,113. 


Zum SchluB seien in Tabelle 9 noch kurz einige Ergebnisse von Beispielen mitgeteilt, die mit 
Hilfe der Gleichung (41) der 1/,—%/,-Punkt-Methode berechnet wurden. Die Form Y(x) der freien 
Wasseroberflache ist in den Abb. 9 bis 15 dargestellt. 

Bei dieser Berechnung wurde nicht das Iterationsverfahren nach Ziff. 6b) 1), sondern die Me- 
thode aus Ziff. 6b) 2) angewendet. Entsprechend ist auch in der letzten Spalte von Tabelle 9 der 
Wert von 4h, angegeben. (vgl. Abb. 3): 


26a 
l ; 
Ah, = 5 | Y (x) dx 
l4a 
Tabelle 9 
Beispiel 6,=6; ae = “1%. Lia a8 af oe a a 
0 0 0 0 

6 —0,1 0,2 20 0,2 Or | --0,29 0,37 + 0,17 
Tt — 0,1 0,2 20 0,1 0,1 | CL 0,28 + 0,15 
8 —0,l 0,05 20 0,1 0,1 | 0,47 0,38 — 0,58 
9 —0,l 0,05 20 0,1 0,05 0,47 0,32 | —0,56 
10 — 0,1 0,03 20 0,05 0,05 0,49 0,37 — 0,28 
1h —0,l 0,03 20 0,05 0,025 0,49 0,31 | —0,28 
2 — 0,1 0,03 20 0,025 0,025 | 0,42 0,32 | —0,22 


Auch aus den in Tabelle 9 enthaltenen Werten erkennt man folgendes: Fiir die Abminderung 
der Zirkulation des Fliigels 1 (gegeniiber dem Zirkulationswert fiir h > 00) ist der Wert von a g/uj 
(Froudesche Langenzahl) noch wesentlicher als der Wert von h, g/uj (Froudesche Tiefenzahl); und 
zwar ist die Abminderung um so gréBer, je gréBer ag/u2 ist. Fiir die Zirkulation des Fliigels 2 
spielt neben den eben genannten Einfliissen noch der Anstellwinkel eine Rolle, der am Fligel 2 
von dem hinter Fliigel J nachlaufenden Wellensystem induziert wird. Dieser Anstellwinkel kann 


je nach der Periode des Wellensystems {d. h. dem Wert von im << positiv (zirkulationssteigernd, 
0 
Beispiele 4, 6, 7) oder negativ (zirkulationsvermindernd, Beispiele 5, 8, 9, 10, 11, 12) sein. 
Aus dem Institut fiir Angewandte Mathematik und Mechanik 
der Deutschen Akademie der Wissenschaften zu Berlin. 


(Eingegangen am 22, Januar 1959.) 
Anschrift des Verfassers: Dr.-Ing. habil. W. H. Isay, Berlin-Dahlem, Schweinfurthstr. 90. 
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Der homogene, innerlich statisch unbestimmte Parallelfachwerktrager 


Von J. Fadle 


1. Einleitung. In diesem Aufsatz soll der Einflu8 von Temperaturanderungen auf die Stabkrafte | 
eines innerlich statisch unbestimmten Fachwerktragers untersucht werden. In dem betrachteten | 
Fall eines homogenen Parallelfachwerktragers kann eine geschlossene Lésung fiir die Stabkrafte | 
angegeben werden. Die Ergebnisse werden an einem Zahlenbeispiel in Diagrammen fiir die Stab- | 


krafte in Abhangigkeit von den Temperaturanderungen At dargestellt. 


2. Berechnung der statisch Unbestimmten. Nach den Regeln der Statik werden aus der statisch 
unbestimmten Grundfigur des Fachwerktragers die n = s + 3—2k tberzahligen Stabe entfernt 
(s die Stabzahl, k die Knotenzahl), wodurch das statisch bestimmte Hauptsystem mit s —n=2k—3 
Stiben erhalten wird. An diesem werden die eingepragten Knotenlasten P,, P,,..., P,_. (wenn 


2 Knoten als Auflagerknoten verbraucht und die ibrigen k — 2 Knoten als belastet angenommen | 
werden) als Lastsystem O angebracht, welche in den s—n Staben die Stabkrafte T; erzeugen. | 


Die Stabkrafte der n iiberzahligen Stabe werden als statisch Unbestimmte X; (i aT Os a ee 


eingefihrt. 

Hierauf werden in den beiden Endknotenpunkten und in der Richtung des i-ten iiberzahligen 
Stabes je eine Einheitszuglast 1 als virtuelle Belastung angebracht, die als das Lastsystem 1 
(i = 1,2,...n) bezeichnet wird und die Stabkrafte t;; erzeugt. Durch Uberlagerung der Stab- 
krafte T; des Lastsystems O mit den mit X; multiplizierten Stabkraften t;; der Lastsysteme 1 
(i = 1,2,...n) kommen die wirklichen Stabkrafte S; zu 


Mit diesen Stabkraften S; kann die Formanderungsarbeit des F'achwerktragers berechnet werden. 
Nimmt man den Einflu8 von Temperaturanderungen einzelner Stabe hinzu, so kommt fiir die 
gesamte Formanderungsarbeit der Ausdruck 


th ep 
i St eS alee (2.2) 
jas e€ 


Dabei bedeuten r; = 1;/E F; die Stabkonstante eines Stabes j, «,, den linearen Warmeausdehnungs- 


koeffizienten und 5’ die Summierung iiber die erwarmten Stabe S,. 
e 


Die Berechnung der statisch Unbestimmten X, erfolgt entweder nach dem Prinzip der vir- 
tuellen Arbeiten oder nach dem Prinzip vom Minimum der Formanderungsarbeit. Nach letzterem 
kommt z. B. aus (2.2) 


0A s oS; 4 . 
me OSS; So eral G=oaeae (2.3) 
Mit 
OX; ue Boe OU 


gemaB (2.1) folgen aus (2.3) nach entsprechender Zusammenfassung die sogenannten Maxwellschen 
Gleichungen 


s \ 


i? 


Sie lauten ausfihrlich angeschrieben 
Besa tii Dare ii; Ahad iD) Gn leeie 
We J =u 
=—| STH + DS eaul, AY (i=1,2,...n) @4 
ef e 


und bilden n lineare Gleichungen fiir die n statisch Unbestimmten X, bis X,, die nach bekannte 
Methoden, z. B. dem Gaufischen Algorithmus aufgelést werden kénnen. 


= (7, TE Xu) tT et eta 0 (teen) 1 


a . . . . 
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j 


3. Berechnung der statisch Unbestimmten sowie der Knotenpunktsverschiebungen mit Hilfe der 
Matrizenrechnung. In neuerer Zeit hat die Matrizenrechnung als das angepaBte Hilfsmittel auch 
in der Baustatik ein weites Anwendungsgebiet gefunden und ist in einer Reihe von Arbeiten! zur 


_Lésung von elastomechanischen Problemen angewandt worden. Deshalb soll im folgenden die 
_ Berechnung der statisch Unbestimmten allgemein in Matrizenschreibweise formuliert werden, wobei 
sich tiberdies eine Aussage iiber die Knotenpunktsverschiebungen des Fachwerks ergibt, die iiber 
das bisher in der Literatur Bekannte hinausgeht. 


Zur Berechnung der statisch Unbestimmten nach dem Prinzip vom Minimum der Formdande- 


-rungsarbeit mu man erst den Ausdruck fiir die Formanderungsarbeit nach (2.2) auf Matrizen- 


schreibweise umformen. Diese lautet 


A=53 RS +048, (3.1) 
wobei Si die Diagonalmatrix 
| re horse 0 
pe te Siege) 
pete ar, 


der Stabkonstanten r;, $’ den (gespiegelten) Zeilenvektor der Stabkrafte und qj, den Zeilenvektor 


der thermischen Lingenainderung gemaB 


SUNS ay seeO) 5 Oy Op At, Ole Aten sii et.) 


bedeuten. Spaltet man in (3.1) die Diagonalmatrix {hi auf in Sh, fiir Stabe des statisch bestimmten 
Hauptsystems und {}i, fiir die tiberzahligen Stabe, so wird aus (3.1) 


A=53RS+ Se Mee +03 (3.2) 
mit? 
$=Up,+Bp,+6r, sv =pW tye sre, 


worin Ql’ = (a; ;), B’ = (b;;),W’ = (¢,;) die gespiegelten Matrizen bedeuten. Die partiellen Differen- 
tialquotienten 


lauten in Matrizenschreibweise symbolisch @A/éx, und man erhalt bei Beachtung der Differen- 
tiationsregel fiir eine quadratische Form, weil ‘i, und i, Diagonalmatrizen sind, 


as! , 08 
op 0 = Gy Sh 8 + UE + On Gy 


Wegen 6’ /cr’ =C’, 0/cy = © folgt nach geeigneter Zusammenfassung 
(CR, C + KR.) x + W RW) p, + GW KR, B) py + ane = 0. (3.3) 


Fiihrt man in (3.3) fiir die als Federungsmatrizen bezeichneten eingeklammerten Produktmatrizen 
die Bezeichnungen 


O'R, C + KR, = Sas > CO R, 2 = Bar» O'R, B= Fae 


ein, so kann (3.3) formal nach x aufgelést werden und liefert mit der reziproken Matrix 33 das 
Ergebnis 


t = — B33 (Bar Px + Ose Diypa- arn ©) « (3.4) 


1 Th. Péschl, Ing.-Arch. 19 (1951), S. 69; J. H. Argyris, Ing.-Arch. 25 (1957) S. 174; H. Baldauf, Ing.-Arch. 
26 (1958) S. 338; R. Zurmihl, Matrizen, 2. Aufl. S. 144 ff, Berlin/Gottingen/Heidelberg 1958. 

2 Vel. R. Zurmiihl a. a. O., Hinzugenommen sind hier Horizontalkomponenten Pjx der Knotenlasten mit 
der EinfluBmatrix 2( = (a;i). 
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Aus (3.2) folgen schlieBlich auch die Vektoren u und 0 der horizontalen bzw. vertikalen Kom- 
ponenten der tk ee als die partiellen Ableitungen 


AON i, 3 + at, = (' RU) py + W RB) p, + QR, C)z + af W 
ape Op, (3.5) 
D = ee -- a Hy 8 ou Min = (B’ Fi W) pe + (BW Ri B) Py HB’ I ©) z+ an B. 
Py a 


Bezeichnet man wieder zweckmaBig die eingeklammerten Produktmatrizen als Federungs- 
matrizen mit den Abkiirzungen 


Su = WR, A), Siz = WH, B) = (BW HW = Far» S22 = (BH, B), 
Sis = WRC) = (C KR, A = Far Sas = (VW HC) = C HR, B) = Fae» 
so erhalt man nach Einfihrung von ¢ aus (3.4) in (3.5) das Ergebnis 


U = (S11 — Bis B33 Gar) Pe + (12 — Sis Gs Waa) Py + LW — Fis Gas oe Mh ! (3.6) 


D = (Gor — Bez 033 O31) Px + (Ge2— M23 033 032) Dy + (B — Bes Bag ©’) Oe, - 


Die Gleichungen (3.4) und (3.6) stellen eine Verallgemeinerung der von Zurmiihl' angegebenen 


Ergebnisse dar, wenn namlich neben vertikalen auch noch horizontale Knotenlasten und auBerdem _ 


noch Temperaturanderungen beriicksichtigt werden. 
Trotzdem die Gleichungen (3.6) einen iibersichtlichen Aufbau zeigen, sind sie fiir eine praktische 


Anwendung wenig geeignet, da (wegen 313 7033 032 = 023 033 031) Insgesamt fiinf Matrizenprodukte 
aus je drei Matrizen gebildet werden miissen, was einen sehr groBen Rechenaufwand erfordert. 
Man kann sie aber durch einen Kunstgriff vereinfachen, indem man aus der zweckmabig zusammen- 


gefaBten Ausgangsgleichung 
0A by 
ar 0 =’ (18 + ay) + HF 
den Klammerausdruck 
Fi, 3 + ay, =—O TRY 


berechnet und diesen dann in (3.5) einfiihrt. Damit kommen die wesentlich einfacheren Glei- 
chungen 
0A 


Y= Obe = QI’ (St, ) ae Oy) aad (Qt oo Je) t 
z (3.7) | 
0A , 1(Cr— 
OS ee (Rt, $ + ay) =— (V’C 1H) x 
bf 


Die mit Hilfe der Matrizenrechnung gefundenen Ergebnisse (3.4) und (3.7) eignen sich besonders | 


fiir eine Auswertung mit elektronischen Rechenautomaten, welche die Elemente von Produkt- 
matrizen und reziproken Matrizen nach festliegenden Programmierschemen errechnen und das 
Ergebnis auBerdem mit einer ausreichenden Stellenzahl ausdrucken. 


4. Der homogene Parallelfachwerktriger. Die Berechnung der statisch Unbestimmten ist im | 
allgemeinen mit sehr viel Rechenaufwand verbunden und erfordert sorgfaltige Rechenkontrollen. | 


Da in diesem Aufsatz nur die wesentlichen Ergebnisse, wie z. B. der Einflu8 von Temperatur- 
anderungen auf die Stabkrafte, herausgestellt werden sollen, beschranken wir uns auf ein ver- 
einfachtes System, das eine geschlossene Lésung zulaBt, und zwar auf einen homogenen Parallel- 


fachwerktrager. Dieser besteht aus gleich langen und parallelen Ober- und Untergurtstaben und — 
je zwei Diagonalstaben in einem Feld. Es ist also in jedem Feld ein iiberzahliger Stab vorhanden, | 
so daB der Grad n der (innerlichen) statischen Unbestimmtheit mit der Anzahl n der Felder iiber- | 
einstimmt. Auferdem sollen alle Stabe gleiche Zugsteifigkeiten EF haben; dadurch vereinfachen | 


sich die Maxwellschen Ess (2.4) auf 


1 R, cnats ara Ol:Se146: 


yo 


7 


Nee 
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wenn 
jee (2, Dal oREeS too, I, At, (4.2) 


als ,,Belastungsglied“ eingefiihrt wird. 
Die Indizes i bzw. k bezeichneten bisher die Nummer des Lastsystems, und da es in einem Feld 
nur ein Paar virtuelle Lasten 1 gibt, bezeichnen i bzw. k auch gleichzeitig die Feldnummern. Die 
s 
Auswertung der Summen >» tj, t;, 1; in (4.1) wird erleichtert, wenn eine Aufteilung in Stabkrafte 
j=l 


des Obergurtes 0,; und Untergurtes u,;, der Diagonalen d,; und Vertikalen v,; links baw. vp 44, ; 


Feld i 


Krafiplan 


Lageplan 


Abb. 1. Belastung des Feldes i durch Priiflast 1. 


rechts vom Feld k vorgenommen wird, die durch das Lastsystem 1 des Feldes i hervorgerufen 
werden. Dabei laufen i und k von 1 bis n. 


Zutolge Abb. 1 ergeben sich fiir die Stabkrafte nach obiger Aufteilung die Werte 


ae fii kL poe . 
Ope een : fiir Obergurt- und Untergurtstibe, 
0 fire eS 
17) fu Po 
a= i: is ‘ fiir Diagonalstabe, 
Oe. tur hoses 


—sina fir hk =i 
sina tur | fur Vertikalstabe. 


Ohi = Vk Li = 
: 0 fir k-+i 
Aus dieser Zusammenstellung berechnen sich jetzt die Summenfaktoren sehr einfach zu 


gine fie kei — 1 und: h=7 1) 


tt SA tg Dy EL = cos & | 
j=1 bed 0 fir k<i—1 und k>i+l1 

s te p ic ei 
jets = (oj; + uk) A + (vi + v7 41,1) A tga + (1? + di) = 


A 


cosa 


sy 


(1 + sin? « + cos? q) . 


Man ersieht daraus, daB zu den Summenfaktoren nur Felder mit einem gemeinsamen Vertikalstab 
Beitrage liefern. Dadurch vereinfachen sich weiter die Maxwellschen Gleichungen (4.1) auf drei- 
gliedrige Differenzengleichungen mit einer ,,Streifenmatrix‘‘. Multipliziert man noch (4.1) mit 
pos 1 50 kommt schlieBlich far (4.1) 


sinka A’ 


cos a B; = ; 
Meet 1) ye a (oneal VAS Arg 0) ee (4.3) 
Die charakteristische Funktion g(a) in (4.3) ist gegeben durch 
23/2 


p(a) = tee at core 1+ (z) + L +: G | (4.4) 


sin® ~ 


wenn die Winkelfunktionen ctg a =A/h, sina = 1/yi+ (A/h)? durch die Feldabmessungen A und h 


ausgedriickt werden. 
22 
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Die Gleichung (4.3) lautet jetzt in Matrizenschreibweise entsprechend (3.3) 


Tag 2 =— (Bee Py tan) =—b. (4.5) 
Hier ist 
2 (a) 1 OF ek. B, 
Sel AP ee a ie ete |eaale (4.6) 
0 0 0. +++ 29(a) B, 


und die Querstriche iiber den Matrizen bzw. Vektoren deuten an, daB diese wegen der Multipli- | 


, > Gn CECE 1 voc » ; 
kation der Ausgangsgleichung (4.1) mit ——-~ eine andere Dimension haben. 
emia sin? a A 


Aus (4.5) folgt der Vektor x der statisch Unbestimmten zu 


xX 
ee a Se ie (4.7) | 
a 


Es handelt sich also nur noch darum, die reziproke Matrix 


mit ihren Elementen g;,, zu berechnen. Dazu gehen wir aus von der Definitionsgleichung der Ein- 
heitsmatrix 


0 fir k +i 


Sas Sal = (f:;) (8 jx) =€ = (¢;,,) mit e, = - é (4.9) 
Lifarik = | 
und der Beziehung zwischen den Elementen einer Produktmatrix gema} 
ein = XSi Bin =F Q, (i, k —e il, 25407) (4.10) 
i 


als skalares Produkt der Zeilenvektoren f' von Kaa mit den Spaltenvektoren q,, von Rene 
Wegen 
fig Si Upp atise y Sit = 2 P@) » fintv=1, ASE ae =. | 


aus (4.6) kommen dann fiir die Elemente g;,, der reziproken Matrix gemaf (4.10) die Gleichungen 


j=s41 

ei, = 2) ij Bin = Bi—1.k +2 PO) Bin + 8it1k=0, (4.11) 
j=i-— 
ee 

Bye Dy Jig Sys = Bi 2 GO) Bis + B46 = 1. (4.12) | 
j=i- 


Weiter ist %33 eine symmetrische Matrix und daher auch $3! eine symmetrische Matrix mit Sik 
= g,;, in der nur Elemente mit 1 > 1 und k < n vorkommen. Demnach sind alle 
£o.—= 9 und £in4+1—9 Gk =12)..2 nH), 
welche Beziehung als Randbedingung fiir die Elemente g;, aufgefaBt werden kann, und die den 
Produktansatz 
(4%, 51 tie 
Sik — (wegen g;, = 8%) (4.13) 


ie neuen Unbekannten x; vom Anfangswert x, = 0 erlaubt. Setzt man (4.13) in (4.11) ein, so 
ommt 


A Xn +1—k [%i—1 + 2 (a) x + 4:41) =0, 


und daraus folgt die Rekursionsformel zur Berechnung der x; 


O33 = (8jx) (hel avon ott) (4.8) | 


Xi41 =—2 (4) %;,—2;_1. (4.14) | 
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Setzt man (4.13) in (4.12) ein, wobei auf Glieder mit i < k und i > k zu achten ist, so folgt 
As Apa My yet 2 rhe (OC) Ap We eda My Xa — i] = 1 


oder nach zweckmafiger Zusammenfassung 


A ie a ea + 2 ya) x;) + x; | = A (A ij 1 — Xin) Sl. (4.15) 
ee 
SEER 


Weil A eine Konstante ist, mu (4.15) fiir beliebige Werte von i erfiillt sein, und es folgt in der Tat 
sowohl fiir 1 = 0 als auch fiir i = n 
Be Le odera: tes ee (4.16) 


%Xn+1- 
Mit (4.16) und (4.13) ergeben sich jetzt die Elemente g;, der reziproken Matrix zu 


Xi Xn+1—k 


ene Bink tI 


Nun sind nach der Rekursionsformel (4.14) wegen x) = 0 alle x; proportional x,, so daB sich x, im 
Zahler und Nenner von g;, wieder heraushebt. Daher kann x, = 1 gesetzt werden, und es folgt 
das Ergebnis 


gin = (4.17) 


Mit x) = 9, x, = +1 und g(a) kénnen aus (4.14) durch Rekursion der Reihe nach die Werte 
Xp, X3,+++%X_4 1 berechnet werden. Damit sind nach (4.17) auch die Elemente g;, der reziproken 
Matrix %33 bekannt, und die statisch Unbestimmten X, folgen aus 


E=— seb oder (X,) = — (gis) (By) 
gem4B (4.7) in geschlossener Form zu 
ee Xi 


Be ee Beet SB Peak (Gal; 2en TO) (4.18) 


*n+1 *n+1k=i+1 


Da in (4.18) iiber k summiert wird, muBte darauf geachtet werden, daB die zu k <i baw. k > 1 
gehérigen Elemente aus (4.17) einzusetzen sind. 


5. Naherungslésung fiir einen mehrfeldrigen Parallelfachwerktrager. Im Falle eines mehrfeldrigen 
Parallelfachwerktragers, bei dem der Grad der statischen Unbestimmtheit entsprechend der Anzahl 
der Felder wachst, laBt sich eine ausgezeichnete Naherungslésung in geschlossener Form herleiten. 
Hierzu gehen wir aus von der Rekursionsformel (4.14), welche mit dem Clebschschen Lésungsansatz’ 


CP OG a: (5.1) 
fiir die neuen Unbekannten r die reziproke charakteristische Gleichung 


r’t2g(a)r+1=0 


mit den Wurzeln 
rm =r =— 90) —/P@)—1, => =— 90) +P) —1 (5.2) 
liefert. Entsprechend den beiden Wurzeln von (5.2) lautet der Ansatz fiir die x; 
C—O pee OM OY inal 
wobei die Konstanten C, und C, aus den Anfangswerten fiir 
ee i Se Or tC, 


1 Beziiglich der Analogie zum homogenen Durchlauftrager vgl. J. Fadle, Ing.-Arch. 17 (1949) S. 317 (Habili- 
tationsschrift). 


22" 
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ZU 
1 
—_—— ive One 
L5G) Tee val 
folgen, und es kommt 
Ree CNT ae) ie (5.4) 


Aus (5.4) folgt fiir i = n, und insbesondere fiir groBe Werte von n, fiir die wegen |r|> 1 das Glied r—” 
vernachlassigt werden kann, die Naherung x, © C,r". Damit ergeben sich schlieBlich fiir die _ 
Elemente g;, der reziproken Matrix aus (4.17) die Naherungswerte 


Xn +1 
oder nach entsprechender Zusammenfassung 
Bee eGov (PY ————— : 
Sk: gy (prt lee) = 1 Taig) pra) —1 (6.5) 
2 Vp2(«) —1 


In (5.5) erkennt man, da® die g;, wegen des Faktors (— 1)*—‘+! ein alternierendes Vorzeichen 
haben. Die Elemente g;; der Hauptdiagonalen kommen mit k =7 zu 


ea, ee ee | 
Pe por ood ed Cali te OI, (5.6) | 
2 Vp2(«) — 1 | 
und sind alle negativ. Aus (5.6) folgt das erste Element fiir 1 = | zu 
gu +r? =—[p(a) —/e() —1]. (5.7) 


Die nachstfolgenden Glieder konvergieren, wie man aus (5.6) leicht erkennt, fiir gréBere Werte | 
von 1 rasch gegen den Grenzwert 
1 
ee ee Deo 
Sii 2 Vp2(«) — = i . ( ) 
Fir i> k ist, wie aus (4.18) hervorgeht, k mit i zu vertauschen, und das bedeutet, daf die 
reziproke Matrix symmetrisch zu beiden Hauptdiagonalen ist. 


6. Ein Zahlenbeispiel. Die Rechenergebnisse sollen an dem Beispiel eines 9-feldrigen Fachwerk- | 
tragers nach Abb. 2 mit 9 ttberzahligen Staben veranschaulicht werden. 


Abb. 3. Statisch bestimmtes Hauptsystem des Briickentrigers mit Lastsystem O:p = 13,5 to, 
P = 31,5 to,undA =4m,L =9A= 36m,h=3m. 


Abb. 3 zeigt das statisch bestimmte Hauptsystem mit dem Lastsystem 0 der eingepragten | 
Lasten, die nur am Untergurt angreifen und aus 8 gleichen Lasten p,, p,,... pyg = p und aus | 
2 Einzellasten P, = P, = P bestehen sollen. 


In Abb. 2 bedeuten i ein Feld des Lastsystems i, k ein Lastfeld und r ein Stabfeld. | 
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. 


Die Knotenlasten werden den Knoten links vom Lastfeld k zugeordnet und in P2;_, fiir den 
Obergurt und P,, fiir den Untergurt aufgeteilt. Zur Kennzeichnung der Lage eines Stabes j im 


Stabfeld r bezeichnen die Indizes 
i i ee Diagonalstabe, 
| j =4r—1 Obergurtstabe, 
} j =4r—2 Untergurtstaibe und 
| ju4r— 3 Vertikalstabe (links vom Feld r). 


Die Elemente der EinfluBmatrix % fiir vertikale Knotenlasten folgen 
fn fiir den Obergurt zu 


| _ @+1—nk-}) 


Tee fir hr 
are ghey = Oy, 499, — 
Pee ry ee a 
ntg « 
- fiir den Untergurt zu 
ae) fir k<r 
b4r—2, 2k—1 = b4r_2, 0% = 
a pas wert 2 ee ye hee 
fiir die Diagonalen zu | (6.1) 
eee Ses 
| b4r,2k—1 = 547,24 = 
; shies dco fir k>r, 
m sin « 
fiir die Vertikalen zu 
Ca ee ete 1 SS iy —— tiie 1p ie 
b4,—3,2k—1 = ‘: b4,—3, 2k = oe 
aes ofa ky, pe fe er 


Die vom Lastsystem O hervorgerufenen Stabkrafte T; in den Obergurt-, Untergurt-, Vertikal- 
und Diagonalstaben eines Stabfeldes r kénnen jetzt mit Hilfe der Elemente 6; ., der HinfluB- 
matrix $ aus (6.1) sehr leicht berechnet werden. 

Fa8t man den Beitrag der 8 Lasten p,,... pig zu den Stabkraften mit Hilfe der Summen- 
formeln 

ky 


i 
Ze) a eka 1) (k, + hk, — 2), 
F@+1—) =1%—h+)RO+)—hy +h 


zusammen, so kommen mit ctg a = //h fiir die Stabkrafte die Ergebnisse 


A 
= ¢—)@+1—72 | 
4 I= 1) Paecart 41 —W Pine] =— Uae | (62) 
2 ee | 
V.= ae Pau cart? — Pakzar=—D,—isina. | 


Um den Einflu8 von Temperaturanderungen auf die Stabkrafte aufzuzeigen, wird angenommen, 
da® im ersten Fall nur der Obergurtstab O,, im zweiten Fall samtliche Obergurtstabe O, bis O, 
um At [° C] erwarmt werden. 
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Mit den Abmessungen 4 = 4 [m], h = 3 [m] von Abb. 2 folgen 
rae, == cob See =(0,6 und cosa = Z 


yee ae 


Damit ergeben sich die Stabkrafte T; nach (6.2), die in Tabelle 1 zusammengestellt sind. 


Tabelle 1. Die Stabkrafte T; des Belastungszustandes O 


1 0 +118,7 | 0 — 148,3 

2 118.7 2919.3 41), Eie80,0 — 125,8 

3 — 219.3 + 302,0 == 159 — 103,3 

4 — 302,0 + 324,7 + 62,0 — 28,3 

5 — 324,7 + 329.3 | +417,0 — 35,8 

6 — 329.3 + 274.0 | + 3,9 + 69,2 

7 e710 00 i) Sie 4 1.5 O17 

8 — 200.7 +- 109,3 — 55,0 + 114,2 

9 — 109,3 0 —68,5 | + 136,7 

10 = = eras 82,0ee4 = 

Die umstindlich zu berechnenden Summen 5 T;t;;1; kénnen fiir den homogenen Parallel- 
fachwerktrager durch einen Kunstgriff auf eine sehr einfache Form gebracht werden. Hierzu be- 
riicksichtigen wir, daB die Stabkrafte t,; der Belastungszustande t i=1,2...n nur.im jeweiligen 


Feld i auftreten, und daher auch die Produkte T;t;,/; nur im Feld 1 von Null verschieden sind. 
Sie ergeben die Summe 


D> Pj til) = — (0, + Uj) A cosa —(V; + Vi41) hsina+D; (6.3) 


cos a” 


Nun folgt aus (6.2) 


== — D;—) sin a, Vi41=— D, sina 


und aus einem Ritterschen Schnitt im Feld i 


Fiihrt man diese Beziehungen in (6.3) ein und setzt h =A tga, so kommt mit (4.4) 


‘cos ea cos a 


sin® « 


1 + sin? « + cos? « 
sin? % 


on 21D, eee 


cos a 


D;-,1+D, 


cos % 


cosa 1 


= B, die vereinfachte Form 
aA 


Damit folgt fiir die Belastungsglieder B, = 


Be = Diag + D; (a) + —— 4 2 EF 0,;1, 04, At. (6.4) 


sin® « 


Der Beitrag eines erwarmten Obergurtstabes zu den Belastungsgliedern berechnet sich fiir das 
gewahlte Beispiel mit 


O%, = 125-10 [1/° C], E = 2,1 10%; [ite/em?]j, F = 263 [cm?], 


Ory =— 0,8 [to] (e =2 ile. aes n) 
Zu 
cosa | 
‘sin? w an EE 06! e Oth At = — 20,455 At [to?] é 


5 
| 
ie 
E 
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SchlieBlich kommen mit 
4 1 in? 3 
(x) Bs + aes cos® % a9 (6.5) 

aus (6.4) die Belastungsglieder fiir die beiden betrachteten Falle: 

1. Stab O; allein erwarmt, zu 2. alle Obergurtstabe erwarmt, zu 

B, = — 11864, B, = — 1186,4 — 20,455 Ar, 

B, = — 1154,7, B, = — 1154,7 — 20,455 At, 

Bo ——— 952.9), B, = — 952,2 — 20,455 At, 

B,=— 329,7, B, =— 329.7 — 20,455 Ar, (6.6) 

B= — ~ 74,7 — 20,455 At, B,=— 74,7 — 20,455 At, 

Be 5478, B, = + 547,8 — 20,455 At. 

B= 802.83 B,= + 802,8 — 20,455 At, 

B, = + 1005,3, B, = + 1005,3 — 20,455 At, 

B, = + 1207.8, By = + 1207,8 — 20,455 At. 


Als letzter Schritt sind die Elemente g;,, der reziproken Matrix entweder aus der strengen Lé- 
sung (4.18) oder aus den Naherungslésungen (5.5) bis (5.8) zu berechnen. Einige Proben zeigen, 
daB fiir den vorliegenden Fachwerktrager mit h < Adie Naherungslésung mit der strengen Lésung 
bis zur achten Dezimalstelle tibereinstimmt. Wir ziehen daher die Naherungslésung vor und 


erhalten mit 


p(x) — Vo2(a) — 1 = 8 —// 63 = 0,06274 , 


21g —1 126 


i — 163 _ 9.96299 


die folgende geschlossene Lésung (6.7) fiir die statisch Unbestimmten X; in Abhangigkeit von den 


Belastungsgliedern B,,: 


X, = — 0,06274 B, + 0,00393 B, — 0,00024 B, + 0,00001 B, — 0,00... B; + 0,00... By — 
a dy 10.400. Bs — 0500's, 

X, = + 0,00393 B, — 0,06299 B, + 0,00395 B; — 0,00024 B, + 0,00001 B; — 0,00... Bs + 
= 0,005" by — 0:00... Fa. -13 0,002 ....B5 

X, = — 0,00024 B, + 0,00395 B, — 0,06299 B, + 0,00395 B, — 0,00024 B; + 0,00001 Bs — 
0.002) ia 0.00. 3. Bs = 0001... Bo 

X, = + 0,00001 B, — 0,00024 B, + 0,00395 B, — 0,06299 B, + 0,00395 B; — 0,00024 Bs + 
+ 0,00001 B; — 0,00... Bz + 0,00... By 

X, = — 0,00... B, + 0,00001 B, — 0,00024 B, + 0,00395 B, — 0,06299 B; + 0,00395 By — 
— 0,00024 B, -+ 0,00001 B; — 0,00... By 

X, = + 0,00... B, — 0,00... B, + 0,00001 B, — 0,00024 B, + 0,00395 B; — 0,06299 Bg + 
+ 0,00395 B; — 0,00024 B, + 0,00001 By 

X, =— 0,00... B, + 0,00... B, — 0,00... B; + 0,00001 B, — 0,00024 B; + 0,00395 Be 
— 0,06299 B; + 0,00395 Bs — 0,00024 By 

X, = + 0,00... B, — 0,00... B, + 0,00... B, — 0,00... B, + 0,00001 B; — 0,00024 Bg + 
++ 0,00395 B; — 0,06299 Bs + 0,00393 By 

X, = — 0,00... B, + 0,00... B, — 0,00... B, + 0,00... B,— 0,00... Bs + 0,00001 Bs — 
— 0,00024 B; + 0,00393 Bs — 0,06274 By 


(6.7) 
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In (6.7) erkennt man deutlich das de Saint-Venantsche Abklingungsgesetz, nach welchem der 
EinfluB entfernter Felder immer schwacher wird. 
Setzt man nun die B, von (6.6) in (6.7) ein, so folgen die statisch Unbestimmten fiir 


1. Stab O; allein erwarmt zu 2. alle Obergurtstabe erwarmt zu 
xe = + 70,14 + 0,00002 At, X, = + 70,14 + 1,208 At, | 
X, = + 64,38 — 0,00032 At , X, = + 64,38 + 1,132 At, 

X, = + 54,45 + 0,00507 At , X, = + 54,45 + 1,137 At, 

X, = + 16,87 — 0,08085 At, X, = + 16,87 + 1,136.At,; (6.8) 
X;= + 5,60 + 1,28858 At, X,=+ 5,60 + 1,136 At, | 
X, = — 31,78 — 0,08085 At , X, = — 31,78 + 1,136 At, 

X, = — 44,70 + 0,00507 At, X, = — 44,70 + 1,137 At, 

X= — 55,54 0,00032 Ar X, = — 55,54 + 1,132 Ar’, 
X, = — 72,00 + 0,00002 At , X, = — 72,00 + 1,208 At. 


Mit den statisch Unbestimmten (6.8) kommen nach (2.1) die wirklichen Stabkrafte S; zu 


Auf eine formelmaBige Zusammenstellung aller S;(4t) wie in (6.6) und (6.8) wird der Kiirze halber 
verzichtet. Dafiir wird aber der Verlauf der Stabkrafte S;(4t) in Abhangigkeit von der Temperatur- 


U(At) 


~ 440 
t 


| 
0(At) | -40 
Erwarmung siimtlicher Obergurtstabe um At | 
——— Erwarmung des Stabes 0, allein um At 


Erwarmung sdmtlicher Obergurtstibe um At 
—— — frwarmung des Stabes 0, allein um At 


Abb. 4, Temperaturabhingigkeit der Obergurtstibe. Abb. 5. Temperaturabhangigkeit der Untergurtstibe. 
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ah XO (At) V(at) 
200 


40 


| 
: 
At 
700°C 
| 
| 
| 
AG 
Ss 
-100 
| 
-120 ‘ 
-140 \ : 
Erwarmung samtlicher Obergurtstibe um At \ | 
——— Frwirmung des Stabes G, allein um At NC | 
Abb. 6. Temperaturabhingigkeit der fallenden Diagonalstibe. Sit INC: | 
3 | 
¢ Os, (At) fo GN 
200 te ee 
«4 
180 ~130 aA 
160 ~149 YyMg 


Erwarmung siimtlicher Obergurtstite um At 
— — — Frwarmung des Stabes 0, allein um At 


Abb, 8. Temperaturabhingigkeit der Vertikalstabe. 


140 


anderung At in den Abb. 4 bis 8 fir Obergurt- und 
Untergurtstiabe, fallende und steigende Diagonal- 
stabe, sowie fiir Vertikalstabe tibersichtlich gra- 
phisch dargestellt. Man erkennt deutlich aus 
den Diagrammen, daf sich die Temperatur- 
anderung eines einzelnen Stabes weniger, hingegen 
diejenige eines ganzen Stabzuges  erheblich 
starker auswirkt und da sogar bei vielen Staben 
ein Umschlagen von einer Zugkraft auf eine 
Druckkraft und umgekehrt eintreten kann. 


Zz ~160 Abb. 7. Temperaturabhingigkeit der steigenden Diagonalstabe. 


Erwirmung simtlicher Oberguristibe um At 
—— — Frwarmung des States 0, allein um At 


(Eingegangen am 27. Januar 1959.) 


Anschrift des Verfassers: Dozent Dr.-Ing. Johann Fadle, Technische Hochschule Karlsruhe, Institut ftir 
Mechanik 
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Die tordierte, einfach gekriimmte Welle mit konstanter Kriimmung* 
Von W. Roth 


1. Einleitung und Aufgabenstellung. Die Torsion einer geraden Welle kann durch ein Torsions- 
moment bewirkt werden, welches proportional der Verdrehung der Welle an der Angriffstelle des 
Momentes ist. Bei einer urspriinglich schon gekriimmten Welle, wobei vorausgesetzt werden soll, 
der Kriimmungshalbmesser R sei konstant, wird eine funktionelle Abhangigkeit zwischen dem 
Torsionsmoment und der Verdrehung bestehen, welche nicht linear ist, und die nicht mehr an Hand 
eines einfachen mechanischen Analogons angegeben werden kann. Es seien jetzt schon die beiden 
Grundaufgaben, deren Lésungen in der vorliegenden Arbeit angegeben werden, prazisiert. 


P (& 


i a 
oO 
| 7\ \e 
f A 
Ya 
My 
Abb. la und 1b, Gekriimmte Torsionsfeder mit Einzellagern. Abb, 2a und 2b. Gekriimmte Torsionsfeder mit Fiihrungsrohr. 


Abb. la zeigt einen gekriimmten Torsionsstab, welcher in seinen Endlagern A und B und seinen 
Zwischenlagern reibungsfrei drehbar so gelagert ist, daB seine Tangenten in A und B erhalten 
bleiben. Unter der Wirkung eines bei A angreifenden Torsionsmomentes M, werden Verdrehungen 
Qq und q, feststellbar sein. Wahrend das Moment M, bei dem gekriimmten Stab eine periodische 
Funktion des Winkels gy, bzw. y, sein wird, ist im Gegensatz dazu bei einer geraden Welle kein 
auBeres Moment erforderlich, um eine Drehung herbeizufiihren. 

Abb. 1b zeigt denselben Torsionsstab mit einer EKinspannung in B, so da dort die Verdrehung 
Null bleiben mu8. Die Abhangigkeit des Momentes M, vom Drehwinkel ¢, ist nun selbstverstand- 
lich eine andere, als bei drehbarer Lagerung in B. 

In beiden Fallen sei vorausgesetzt, daB die Welle hinreichend oft zwischen A und B gelagert ist, 
so daf} sowohl ein Ausweichen aus ihrer Ebene heraus als auch eine Kriimmungsanderung in ihrer 
Ebene unméglich ist. Gesucht sind die Charakteristiken der Probleme. 


2. Aufstellung der Gleichungen. Die oben gegebene Aufgabenstellung gestattet beziiglich der 
mathematischen Behandlung eine kleine Mutation. Es wurde vorausgesetzt, da durch hinreichend 
viele Lager ein Ausweichen der Welle unméglich geworden ist. Idealisiert heiBt das, die Welle 
soll in einem genau passenden Rohr liegen, in welchem eine reibungsfreie Drehung méglich ist. 
Dieses 4uBere Rohr ist als starr zu betrachten, so da es imstande ist, die von der Welle herriihren- 
den Krafte aufzunehmen, ohne daf eine Deformation des Fiithrungsrohres eintritt. Abb. 2a und 2b 
zeigen schematisch diese Analogon-Probleme. Beide Probleme unterscheiden sich nur durch die 
Randbedingungen voncinander. Es ist also das Ziel, zunachst allgemeine Gleichungen aufzustellen, 
unabhangig von den verschiedenen Aufgabestellungen. 

Im unbelasteten Zustand wird die Lage eines Elementes der Welle durch die Koordinate « 
beschrieben. Nach Aufbringung des Torsionsmomentes M, hat sich das Element um die Strecke u 
in Richtung der Wellentangente verschoben. Gleichzeitig wird eine Verdrehung g des Elementes 
um die Tangente als Drehachse zu beobachten sein. Sowohl u als auch g miissen als Funktionen 
von a aufgefaBt werden. Wahrend die Drehung ¢ keiner einschrankenden Bedingung unterworfen 
wird, soll von der Verschiebung u gefordert werden, daB® sie klein bleibt. Dieser Sachverhalt ist 
dann verstandlich, wenn man bedenkt, daB die Verschiebung u verursacht wird durch die Langs- 


* Dissertation Karlsruhe 1958, Hauptreferent: Prof. Dr.-Ing. R. Sonntag, Korreferent: Prof. O. Kraemer 
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_ kraft S in der Welle, zwischen S und u aber der Hookesche Zusammenhang wie bei einem geraden, 


on\ 


auf Zug oder Druck beanspruchten Stab besteht. Fir Metallstibe ist es aber erwiesen, daB die 
Verschiebungen u so klein bleiben, daB sie nicht in die Gleichgewichtsbedingungen des Elementes 
eingehen. 

In Abb. 3 sind nun samtliche an dem zwar um den Winkel y gedrehten, aber nicht verschobenen 
Element angreifenden Krafte und Momente eingezeichnet. Die Komponentenrichtungen sind 
orthogonal aufeinanderstehend gewahlt und zwar so, daB die erste Komponentenrichtung mit dem 
Tangentenvektor t, die zweite mit dem Hauptnormalenvektor 11 und die dritte mit dem Binormalen- 
vektor b des die Mittellinie der Welle begleitenden Frenetschen Dreibeins zusammenfiallt. Das Drei- 
bein ist nach einem Rechtssystem orientiert, so daf8 die Binormale 6 aus der Wellenebene heraus 
auf den Beschauer hin zeigt (Abb. 4). Die an den Schnittflichen des Wellenelementes angreifenden 
inneren Krafte und Momente sind derart angenommen, daf ihre positiven Richtungen an der 
Stelle « + da mit den Einheitsvektoren t + dt, n + dn und 6 + db zusammenfallen. Ihrem 
Charakter entsprechend zeigen sie dann an der Stélle « entgegen den Richtungen t, n und b. Die 
Momente drehen in einem Drehsinn, welcher sich aus der Definition des Rechtssystems ergibt. In 
Pfeilrichtung gesehen dreht ein positives Moment also rechts herum. Es bedeuten p, die Riick- 
wirkung des Fiihrungsrohres auf die Welle bezogen auf die Langeneinheit der Wellenmittellinie 


Bp, as 


Abb. 3. Krafte und Momente am Wellenelement. , Abb. 4. Begleitendes Dreibein an der Wellenmittellinie. 


in Richtung der Normalen n, p, die Riickwirkung des Fiihrungsrohres auf die Welle bezogen auf 
die Langeneinheit der Wellenmittellinie in Richtung der Binormalen 6, S die Langskraft in der 
Welle, Q, die Querkraft an der Welle in der Wellenebene, Q, die Querkraft an der Welle senkrecht 
auf der Wellenebene, M das innere Torsionsmoment um die Tangente t als Achse drehend, M, das 
innere Biegungsmoment um die Normale 1 als Achse drehend und M, das innere Biegungsmoment 
um die Binormale } als Achse drehend. Diese Krafte und Momente sind in Abb. 3 an der Steile « 
bzw. a + da so, wie sie am Wellenelement wirken mit ihren Differentialen eingezeichnet. Man 


liest aus Abb. 3 die Gleichgewichtsbedingung gegen Verschiebung 


d (tS) +d(nQ,) +d(6Q,) +n pyds +6 py ds =0 (1) 
sowie die Gleichgewichtsbedingung gegen Drehung 
d(t M) + d(u M,) + d(6 M,) + (t x n) Q, ds + (t x b) Q,ds = 0 (2) 


des Elementes ab. Bevor die Komponenten dieser Gleichungen angegeben werden kénnen, miissen 
iiber die Frenetschen Formeln 

dt du qd 
die Ableitungen der Einheitsvektoren t, 11 und 6 ermittelt werden. Da die Welle durch das Fiihrungs- 
rohr gehalten wird, kann die Wellenmittellinie ihre Gestalt nicht andern. Im unbelasteten Zustand 
ist die Wellenmittellinie aber ein ebener Kreisbogen mit dem Kriimmungshalbmesser R. Daraus 


folgt sofort, daB die erste Kriimmung 


on 
und die zweite Kriimmung 
LO 
ist. Die Frenetschen Formeln gehen mit der Beziehung ds = R da damit iiber in 
dt dn eS Lae 


ee ete 
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Beachtet man noch daB 

(tx n)=) (t x b) =—n 
ist, so gehen die Gleichungen (1) und (2) nach Division mit da iiber in die Komponentengleichungen 
fiir die Richtungen t, n und b 


dS dQ, oA dQ, ae A Us | 
Fe uO ge es Di Oe eee ee 
dM LL dM are s | 
——M,=0, S1+M—Q@R=0, 544+0,R=0. ( } | 


Diese sechs Gleichungen garantieren hinreichend und notwendig das Gleichgewicht des Elementes. 
Eine Auflésung indessen nach den einzelnen Kraften und Momenten ist noch nicht modglich, da die 
Zahl der Unbekannten nicht mit der Zahl der zur Verfiigung stehenden Gleichungen iibereinstimmt. _ 
Die Schwierigkeit ist behoben, wenn die Schnittkrafte und Schnittmomente in Beziehung zu den © 
von ihnen verursachten Verschiebungen u und p gebracht werden. Zwischen S und u besteht die 
Hookesche Gleichung 

E F du 
2 Rada 8) 
mit dem Elastizitatsmodul E und dem Wellenquerschnitt F. Ebenso kann die Abhangigkeit zwi- 
schen dem Torsionsmoment M und dem Torsionswinkel dp, um den der Wellenquerschnitt an der 
Stelle « + da relativ zum Wellenquerschnitt an der Stelle « gedreht wird, durch die elementare 
Formel fiir einen Torsionsstab von der Lange ds 


mit dem polaren Flachentragheitsmoment I, des Wellenquerschnitts und dem Schubmodul G 
beschrieben werden. Es gilt also 

dp = MR 

dg melee a 


Blickt man an der Stelle a in Richtung der Tangente t auf den Querschnitt, so soll eine Drehung p 
entgegen dem Uhrzeigersinn positiv sein. Ein positives Torsionsmoment M verursacht dann aber 
ein negatives Differential dy. Daraus resultiert das Minuszeichen in Formel (4). 

Nicht mehr so einfacher Art sind die Beziehungen zwischen dem Drehwinkel p und den inneren 
Biegungsmomenten M, und M,. Hier muf das um den endlichen Winkel ¢ ,,in sich“ gedrehte 
Wellenelement betrachtet werden. Ein Flachenelement des Wellenquerschnittes habe im un- 


Abb. 5. Lageanderung eines Flachenelementes des Wellenquerschnittes. Abb. 6. Faserelement der Welle vor und nach der Deformation. 


belasteten Zustand die Polarkoordinaten r und y, im belasteten Zustand die Koordinaten r und 
y + @ (Abb. 5). Die Drehung des ganzen Elementes um den Winkel ¢ hat zur Folge, daB eine 
Dehnung des zum Flachenelement dF’ gehérigen Faserelementes eintritt. Die Lange des Faser- | 
elementes vor der Drehung (Abb. 6) ist 


ds, = (R—rsiny) da 
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und nach der Drehung 
ds, = [R—rsin (y + 9)] da. 


Genauer ist 


ds, = da |/[R—r sin (y roe + (zy § 


_ doch kann von dem Einflu8 der relativen Drehung dp der beiden Querschnitte zucinander abgesehen 
werden, da angenommen oder vorausgesetzt werden soll, daB 
[R—rsin(y + @)P> 7 a) 
gilt, eine Voraussetzung die sicher dann erfiillt ist, wenn die Ableitung dy/dx der Bedingung 
ie roe 1| (5) 


a 


geniigt. Damit ist also eine Dehnung 
j ds,—ds, ___ r[siny —sin (v + 9)] 


ds, R—r-siny 


feststellbar. Die zu dieser Dehnung gehérige Spannung ist aber nach dem Hookeschen Gesetz 


= r [sin y — sin (vy + ¢—)] 
TAS aaa eer aes Me 


Da nur verhaltnismaBig diinne Wellen betrachtet werden sollen, gilt fir die Quotienten 
r/R < a/R <1, so daB die Entwicklung der Formel (6) nach Potenzen von r/R nach Abbruch der 
binomischen Reihe mit dem linearen Glied die einfache Formel 


o = Ez [siny —sin (y + @)] (6’) 


liefert, welche maBgebend ist fiir die Ermittlung der in der Welle sich einstellenden maximalen 
Zu- oder Druckspannung. Fir in Elemente mit den Polarkoordinaten r =a und. y = 7/2 baw. 
y =32/2 ist nach einer Drehung um g = @ eine Spannung 


2a 
Omag == 215 wR 


E 


feststellbar, welche also bei einer wirklich ausgefiihrten Welle kleiner als die FlieBspannung des 
Materials bleiben mu. Die Formel (6) ist im wesentlichen auch schon in der von R. Grammel! 
dargelegten Theorie tiber ,,Das Umstiilpen und Umkippen von elastischen Ringen“ zu finden. 
Dort wird gezeigt, daB sich ein kreisférmiger Querschnitt des Ringes nicht genau um seinen Schwer- 
punkt oder Mittelpunkt dreht, sondern um einen ,,neutralen Punkt welcher spannungsfrei bleibt 


und der in der Entfernung 
R a 
s=5 (1 yj Ss 


vom Schwerpunkt des Kreisquerschnittes zum Kriimmungsmittelpunkt des Ringes hin gelegen ist. 
Bei Beriicksichtigung nur linearer Glieder in a/R, wie es hier durchgefiihrt wurde, ist s = 0. Die 
Drehung erfolgt also dann doch hinreichend genau um den Kreismittelpunkt oder Schwerpunkt 
des Querschnittes. 

Die Spannung o wird erzeugt durch die inneren Biegungsmomente M, und M,. Waren diese 
Momente namlich nicht vorhanden, so wiirde sich das Wellenelement wie ein starrer Kérper um 
seine Achse drehen, d. h. nach einer Drehung um den Winkel ¢p stiinden die Querschnitte an den 
Stellen « und « + do nicht mehr unter dem Offnungswinkel da zum Kriimmungsmittelpunkt der 
Wellenmittellinie zueinander geneigt. Die Hypothese aber, da die Querschnitte im belasteten 
Zustand ihre im unbelasteten Zustand eingenommenen Ebenen nicht verlassen, ist ebenso einfach 
und plausibel wie die Annahme, welche in der elementaren Festigkeitslehre getroffen wird, daB die 
Querschnitte eines auf Biegung beanspruchten Tragers eben bleiben. Zwischen den Biegungs- 
momenten M, und M, und der Spannung o bestehen demnach die Beziehungen 


F BF 
M, = ff orcos(y +9) dF, M,=—fforsin(y +9) dF. 
1 R.Grammel, Z. angew. Math. Mech. 3 (1923) S. 429. 
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Fiir das Flachenelement dF ist dF =r dr dy zu setzen. Die Integration erstreckt sich tiber den 
gesamten Wellenquerschnitt. Mit (6’) erhalt man nach Auswertung der Integrale 


M, =— 52 sing, (7) 
ey) ro 8 
M, = a (1 — cos q) , (8) 


wobei I, = x a‘/2 gesetzt ist. Diese Gleichungen wurden auf anderem Wege auch von R. Sonntag? 
hergeleitet 

Damit sind die Zusammenhange zwischen den Biegungsmomenten und der Drehung des Wellen- 
elementes gefunden. Die Momente M, und M, kénnen als Komponenten eines resultierenden 
Biegungsmomentes M, aufgefaBt werden, welches dem Betrag 


nach 


ara ie 1 Gy ce 
M, = (Mt Mee i + (9) 


ist (Abb. 7) und dessen Drehachse mit der Hauptnormalen nt den 
Winkel 6 einschlieBt, der aus der Gleichung 


Mz (cos a ees 
Sai eat = tg5 


sin p 


folgt. Man erhalt 


> 


P 
Abb. 7. Lage des resultierenden Momentenvektors O) — 
im Wellenquerschnitt. 


und diese Gleichung sagt aus, daB sich die im Sinne der elementaren Biegungslehre definierbare 
neutrale Faser, welcher die Punkte mit der Biegungsspannung o = 0 zugeordnet sind, in gleichem 
Drehsinn mit der Drehung p des Querschnittes um q/2 dreht. 

Die Gleichungen (4), (7) und (8) in Verbindung mit den Gleichungen (1’) und (2’) reichen zur 
Bestimmung der unbekannten Krafte und Momente aus. Das Problem wird am einfachsten auf die 
Ermittlung der Funktion y = f(a) zuriickgefiihrt. Die maBgebende Gleichung ergibt sich, wenn die 
Ausdriicke (4) und (7) fiir die Momente M und M, in die erste der Gleichungen (2’) eingesetzt 
werden. Man erhalt fiir ~ die Differentialgleichung zweiter Ordnung 


dy d Dis ae 
aCe a ant es oe 
Ein erstes Integral dieser Gleichung ist 


dp\2 _ E 
(a) = C— FZ cosp. 


Fiir beide Grundprobleme kann dy/dx < 0 vorausgesetzt werden. Damit liegt die fiir die weitere 
Untersuchung wichtige Gleichung 


dp _ PAE Re , 
# Vc Feosp (10°) 


vor. 


3. Lésung der ersten Grundaufgabe. Die drehbare Lagerung der Welle bei B fordert M=0 | 


an der Stelle « =o. Uber (4) folgt daraus, daB auch dy/dx = 0 sein muB fiir « = &. Gleichung 
(10’) liefert damit die Integrationskonstante 


und (10’) selbst geht iiber in 


d Cy ee 
E=—|/F Ve08 He — e089 (11) 


1 R. Sonntag, Forsch. Ing.-Wes, 23 (1957) S. 216. 
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_ An der Stelle « = 0 ist das Torsionsmoment M — M a der Winkel m =: y,, somit folgt aus Glei- 
chung (11) iiber (4) 


M.R 
Ip VEG 


= cos y, — cos @, , (11’) 


und die Integration von (11) fihrt auf 
” 


4 G dy 
Oss Ve | ——., 12 
; ‘ E $ /cos p, — cos p 7) 


Die geometrische Bedingung 


Pb 


peo (ea an) 


Jcos a cos Pp 


in Verbindung mit Gleichung (11’) reicht zur Festlegung der Charakteristiken M, = F'(p,) und 
M, = F,(%,). Es ist jedoch zweckmafig, eine Umformung vorzunehmen. Die Substitution 


3 cos y/2 
sin y se, 
fihrt (12) tiber in 
ce y 
E dy 
af—= | - —— 13 
26 V1 — k? sin? y (13) 
‘1 
und (12’) in die Gleichung 
Pe 7/2 
E dy / 
a —= = 8 13 
°Y 26 | )1 — k? sin? y 2) 
af 
wobei, 
4 2 
4, = arc sin SO Pa)? -) Ves are gin eee 
und 
k = cos y,/2 
ist. 


Nun sind die Elliptischen Integrale erster Gattung nach Legendre definiert durch 


o 
dx 
0 


und speziell ist 


n/2 d 
() | 1 — F? sin? x 
. 0 
Damit geht (13) tber in 
pe . cos p/2 oA . cos @,/2 14 
a (2 = F (are sin, k F (are cine k (14) 


und (13’) in 


: Ps 
Xo Vee = K(k) — F (are sin eg 


k |. (14’) 
Beide Gleichungen erlauben eine explizite Auflésung. Man erhalt aus (14) 


Te uy, 
cos + = ksinam |. Le == F (are sin ree ‘ k)| (15) 


k 
cos“ = ksinam 


5 K(k) — % ee | , 
Die letzte Gleichung liefert den Winkel g, tiber k als Funktion des Winkels g,. Die GréBe 
cE 
Bo = % 26 


und aus (14’) 
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iibernimmt dabei die Rolle eines Parameters. Aus Gleichung (11’) erhalt man dann das dimen- 


sionslose Moment 


Oe SEES 
EO Tey hC. 
IU 
Bo = con fo =a 
Q, 
1,0000 0,00 0,00 0,000 0,00 0,000 
0,9962 10,00 2 4,90 0,278 106,88 1,129 
0,9848 20,00 46,70 0,504 178,64 1,392 
0,9063 50,00 109,68 0,990 257,66 0,926 
0,8192 70,00 137,96 1,042 266,30 0,638 
0,6428 100,00 164,14 0,888 254,76 0,298 
0,4226 130,00 176,24 0,595 229,40 0,089 
0,1737 160,00 179,78 0,245 200,00 0,018 
0,0000 180,00 180,00 0,000 180,00 0,000 
—0,1737 200,00 180,22 — 0,245 160,00 — 0,018 
— 0,4226 230,00 183,76 — 0,595 130,60 — 0,089 
— 0,6428 260,00 195,86 — 0,888 105,24 — 0,298 
— 0,8192 290,00 222,04 = 1,042 93,70 — 0,638 
— 0,9063 310,00 250,32 — 0,990 102,34 — 0,926 
— 0,9848 340,00 313,30 — 0,504 181,36 Sey 
— 0,9962 350,00 335,10 Soe 253,12 SMe 
— 1,0000 360,00 360,00 0,000 360,00 0,000 


in Abhangigkeit vom Winkel g, oder q. 
Es wurden die Tabellen fiir elliptische Funktionen verwendet, wie sie z. B. in der Hiitte! zu 


finden sind. Die Zwischenwerte wurden durch lineare Interpolation ermittelt. Abb. 8 und 9 zeigen 
fiir einige Parameterwerte /, den Winkel yy, iiber q, aufgetragen, in Abb. 10 und 11 ist Q, in Ab- 
hangigkeit von y, dargestellt. Einige charakteristische Eigenschaften dieser Kurven lassen sich 
durch elementare Formeln beschreiben. Zunichst sind fiir hinreichend kleine Momente Q, die 


Winkel gy, y, und gy, ebenfalls klein, so daB 


Pp 2 
cos p = 1—* cos g, = 1—% 
gesetzt werden darf. Gleichung (12’) hat damit die Form 
Pa 
E dp 
WE - fot 
2G ap © 
eg ae 

und die Integration liefert 

Po = Po Go] Bo - (16) 


Also ist 

dp, 

Pe — Co} By | 
an der Stelle g, =g, =9. Fiir verschwindende Momente Q, sind neben den trivialen Fallen 
Po = = 9 und oY, =% =2%, die Winkel y, = y =a einander zugeordnet. Die Welle be- 
findet sich in dieser Lage in einem instabilen Gleichgewichtszustand. Jeder Querschnitt ist um den | 
Winkel g =a gedreht. Man erkennt dies, wenn man mit dem Lésungsansatz 

(on es 
welcher die Randbedingungen p, = y, =a erfiillt, in die allgemeine Differentialgleichung (10) 
fiir den Winkel — eingeht. Diese Gleichung ist offenbar befriedigt; es ist die AuBenfaser der Welle 
zur Innenfaser geworden und umgekehrt. Die geringste Stérung in dieser Stellung jedoch wird | 
dazu fiihren, da die auf Druck beanspruchten Innenfasern und die auf Zug beanspruchten AuBen- 
fasern die Welle umschnappen lassen. In der Umgebung dieser Gleichgewichtslage kann nun | 


Y=a+e, M%=t+N, g=a+A, dy = di 
gesetzt werden, wobei A, ¢ und 7 klein sind, also auch 
2 
coop =— 145 cosy, =—1 +% 
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%a 


0° 50° 780° 270° 30° 


Abb. 8. Abhangigkeit der Verdrehungen gg und gy der Wellenquerschnitte. 


4 fa 
560° 


270° 


- 780° 


90° 


YP, 


i? 90° 780° 270° 360° 


Abb. 9, Abhangigkeit der Verdrehungen gq und yp der Wellenendquerschnitte. 
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ist. Gleichung,(12’) geht so tiber in 


12 = dh 
Ko eh cmaes hese a 
n 


und die Integration fiihrt auf die einfache Beziehung 


€=7 cos Py, 
welche aussagt, dal} 
; dy 
lim > = “%2 — cos 
n—>0 U] dp, Po 


(17) 


im Punkt y, =, = ist. Die Steigung der Kurve q, ist also -- 1 fiir By = 0 bzw. a und ver- 
schwindet fiir By = 2/2. Abb. 9 laBt dieses Verhalten deutlich erkennen. Endlich ist nach einer 


S2q 


45 


O5 


~10 


Bs i Fa 
Safe 90° 780° 270° 36G° 
Abb. 10. Charakteristikenfeld Q, der Wellen mit drehbaren Endquerschnitten. 
~49 5 ee 


0 90° 780° 270° 360° 
Abb, 11. Charakteristikenfeld Qg der Wellen mit drehbaren Endquerschnitten. - 


ad 
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ganzen Umdrehung der Welle wieder derselbe Zustand wie im Nullpunkt y, = gm, = 0 erreicht 
Demnach ist i ; 


dp, & 
dp, = Co) Bo 
fiir Pa = % = 27, eine Gleichung, die natiirlich auch wieder durch Grenziibergang hatte her- 
geleitet werden kénnen. AuBerst merkwiirdig ist der Verlauf der eigentlichen Charakteristiken 
02, der Welle. Fir kleine y, und q, erhalt man iiber (16) aus (11’) 


eine Formel, welche wegen der Periodizitat des Momentes M, auch nach einer Umdrehung der 
Welle, also fir p, =~, = 22 Giiltigkeit hat. In der instabilen Gleichgewichtslage ist nach (11’) 


: ; dp 
dQ, ae age ae - 0 
dp, 22. io) 
die Ableitung der Kurve 2,. Nach Anwendung der Bernoullischen Regel erhalt man 
dy, \? dy; 
dQ cos Va Sac O'S: Pp (32) eh Oe “gt 
dp, 9 124 
dp, 


und daraus unter Verwendung von (17) 

eo aaa! 

ip ep 
Diese Gleichung sagt aus, dafB die Steigung der Kurve 22, verschwindet fiir die Parameter Ga =O 
bzw. z und daf sie unendlich ist fiir 6, = 2/2. Schaubild (11) bestatigt diese Aussage. Wahrend 
der Verlauf der Kurven Q, fiir 0 < 6, < 2/2 ohne weiteres verstandlich ist, sind die Zusammen- 
hange bei gréBeren Parameterwerten /) komplizierter. Fiir 7/2 < fy) <a wachst mit zunehmendem 
Winkel gy, zunachst auch das Moment 2, bis zu einem Maximum; wahrend ¢, nun weiter zunimmt, 


beginnt 2, zu fallen. SchlieBlich ist an der Stelle ~, = @q maz die Ableitung 
dQ, 
Mp, 
Wie die Gleichung (11’) nach Differentiation zeigt, ist die Stelle ~, m4. gekennzeichnet durch die 
Aussage 


= —o. 


ae)? 
dp, 


Die Welle darf also nicht itiber den Winkel 9, ,,., hinaus gedreht werden. Mit weiter abnehmendem 
Moment Q, verringert sich nun auch wieder @,, bis bei verschwindendem Moment die Welle ihre 
instabile Gleichgewichtslage y, =Q, =z erreicht hat. Der weitere Verlauf der Charakteristik 
verhalt sich zu dem soeben beschriebenen Kurvenast zentralsymmetrisch, so da er bei einer 
Drehung um 180° um den Punkt g, =z als Drehpunkt mit diesem zur Deckung gebracht werden _ 
kann. ZusammengefaBt ist das Wesentliche also, daB bei einem monoton zunehmendem Winkel gy 
von Null auf x eine Drehung 9, eingeleitet werden muf, welche gréBer ist als 7; dann erst mit 
abnehmendem Q, eine Riickdrehung ¢, einsetzt, so daB also die Winkel p, und y, mit entgegen- 
gesetztem Drehsinn in die Gleichgewichtslage y, = ~ = einlaufen. Anschaulich gesehen werden 
durch die Uberdrehung gy, in der Welle bei A groBe innere Torsionsmomente erzeugt, welche nétig 
sind, um das mit seiner Drehung q, zuriickbleibende Wellenende bei B zu einer fortgesetzten 
Drehung zu veranlassen. Damit aber keine dynamischen Effekte entstehen kénnen, mu schlieBlich 
Q, plétzlich sehr stark abnehmen und die Rickdrehung einsetzen. Fiir praktische Zwecke diirften 
Parameterwerte [, > 7 nicht mehr interessieren; aus diesem Grunde werden die Charakteristiken 
fiir diese Parameter auch nicht weiter untersucht. 
23" 
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Ein ahnliches Verhalten wie die hier untersuchte gekriimmte Welle weist auch die von R. Gram- 
mel! behandelte Kreisbogenfeder auf. Infolge der ganzlich anderen Problemstellung, ist auch die 
von Grammel angegebene Lésungsmethode von der hier aufgestellten véllig verschieden. Die 
Kreisbogenfeder ist nur in ihren Endpunkten A und B gezwungen, ihre urspriinglichen Tangenten 
Geisubehalieet wahrend sie sich zwischen diesen Endpunkten zu einer raumlichen Kurve verbiegen 
kann. Der mathematische Aufwand zur Lésung dieses Problems ist demnach wesentlich kompli- 
zierter. Es ist indes naheliegend zu vermuten, daB bei nicht zu groBen Offnungswinkeln a, die 
verschiedenartigen Lagerungen der Welle nur wenig HinfluB auf die Charakteristiken haben kénnen. 
Dies beweisen auch eindeutig die numerischen Ergebnisse. Nicht zuletzt ist damit aber auch wieder 
gezeigt, da es erlaubt ist ein Problem der Mechanik so abzuandern, daB ein mathematisch gefii- 
gigeres Problem entsteht, dessen Lésung mit weniger Aufwand erzwungen werden kann. 

Theoretisch aufschluBreich ist die in den Abb. 10 und 11 eingezeichnete Grenzkurve fiir 6) = 00. 
Sie laBt sich durch eine einfache Uberlegung herleiten. Der Offnungswinkel «) geht mit dem Para- 
meter fi, gegen Unendlich. Es liegt also eine Schraubenfeder mit der Steigung Null und unendlich 
vielen Windungen vor. Bei einer Drehung der Welle um den Winkel y, bei A wird bei unendlich 
langer Wellenlange der Drehwinkel gy, bei B Null bleiben. Aus dieser Uberlegung resultiert sofort 


tiber Gleichung (11’) die Formel 
M,R 
a 1 — : 18 
= es won, (18) 


Schon jetzt kann gesagt werden, daB die durch diese Gleichung bestimmte Charakteristik iiber- 
einstimmt mit der Momentenkurve fiir die bei B unverdrehbar eingespannte unendlich lange 
Welle, da bei einem Drehwinkel y, = 0 und auch einer Ableitung dy/da = 0 in B, wie es bei der 
unendlich langen Welle der Fall sein wird, kein Einspann-Torsionsmoment auftreten kann. 

Explizite einfache Formeln erhalt man auch fiir kleine Parameterwerte {,. Hier liefert die An- 
schauung sofort einen Anhaltspunkt. Im Grenzfall «1, = 0 liegt namlich ein gerades Wellenstiick 
vor mit endlicher Lange und unendlich grof{em Kriimmungshalbmesser R. Die Aussagen gy, = 
und {, = 0 sind trivial. Fir kleine a, kann demnach 


| te=%t%, va=mte, dpade 
gesetzt werden, wobei ¢ und ¢é) kleine Winkel sind. Damit ist auch 
COS Y = COS Y, —E SING, » 


und Gleichung (12) geht tiber in 


er (Pa—"b) 
i de 
a. = ——== | wae (19) 


/2 sin y, (een Ve 
Gleichung (12’) bekommt die Form 


(Pa—%b) 
de 
x — Se , 
Waa “aaee F Ve. a 
Die Integration von (19’) liefert die Beziehung 
2 
Pa = 1 ae sin @, - (20) 


Obwohl diese Gleichung in y, transzendent ist, kann eine Auflésung nach @ durchgefiihrt werden. 
Setzt man 


Po = Va 208 


und beachtet man, dafs konsequent zu den oben angegebenen Naherungen 
Sin YP, = Pa — Ep COS Pg 


C= = sin @, » (21) 


ist, so erhalt man aus (20) 


nachdem noch in f§ linearisiert wurde. Die Charakteristiken kénnen nun auch naherungsweise 
angegeben werden. Unter Verwendung von (21) folgt tiber (11’ ) 
M,R Bo . 
___ — -—_ = gin 9 
ie VE G y2 Pa (22) 
1 R. Grammel, Ing.-Arch. 9 (1938) S. 126. 
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3 
r. 
2 


und unter Verwendung von (20) 


Ip VEG 

Die Gleichung (22) fiir die fast gerade Welle wurde von R. Sonntag schon angegeben und auf vollig 

anderem Wege hergeleitet!. Es ist nun méglich, einen Vergleich der Resultate, wie sie die strenge 

_ Theorie liefert, mit den durch obige Gleichungen (20) und (22) gegebenen Werten, anzustellen. 

- Abb. 12 und 13 zeigen fiir den Parameter £, = 0,50, dem mit E = 2. 106 kg/cm? und 

~ G =8.- 105 kg/cm? der Offnungswinkel a) = 25,60° zugeordnet ist, gestrichelt eingezeichnet die 
'Ja- Pe 


4 8 al 


(5 


ye 


2 


~-y° 


Z -6° 


i ae eee kD AN SP ORE ENE GS 


0° 90° 780° 270° 360° 


Abb. 12, Exakter Verlauf und Naherung fiir den Differenzwinkel pg — yp bei kleinem Parameter fh). 
{e 
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Loa 
~ ees ene 
aa 50° 780° 270° 560° 


Abb. 13, Exakter Verlauf und Naherung fiir das Moment Q, bei kleinem Parameter {). 
1 A.a. O. Gleichung (11). 
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nach (20) und (22) errechneten und ausgezogen die durch die exakte Rechnung erhaltenen Kurven. 
Es ist ersichtlich, daB fiir Parameterwerte 8, < 0,50 die Naherungen den durch praktische Belange 
gestellten Genauigkeitsanspriichen in auBerordentlich hohem Ma&e gerecht werden. 

Sind nun die elastischen Eigenschaften der gekriimmten Welle analysiert, so sollen noch einige 
Uberlegungen iiber die auftretenden Krafte und Momente engestellt werden. Damit die Welle 
die im tordierten Zustand ebene, kreisbogenformige Gestalt beibehalt, miissen die vom Fiihrungs- 
rohr herriihrenden Krafte p, und p, auf sie einwirken. Zur Ermittlung dieser Belastungen als auch 


der inneren Krafte und Momente der Welle, miiBte zunachst Gleichung (15) ausgewertet werden. | 


Aus ihr erhalt man den Drehwinkel y in Abhangigkeit des Argumentes 
a 

Bae Ver 
dessen Einfiihrung an Stelle von « sich empfiehlt. Ist so m als Funktion von a oder 6 bekannt, so 
sind auch iiber die Gleichungen (4), (7) und (8) die Momente M, M, und M, gegeben. Wahrend die 
erste der Gleichungen (2’) bereits zur Aufstellung der Differentialgleichung (10) verwendet wurde, 
erhalt man aus der zweiten Gleichung (2’) die Querkraft Q, und aus der dritten die Querkraft Q,. 
Die erste Gleichung (1’) ist eine Differentialgleichung erster Ordnung fiir die Langskraft S, deren 
Lésung aber sofort angegeben werden kann, wenn aus dieser Gleichung und der dritten Gleichung (2’) 
die Querkraft Q, eleminiert wird. Man erhalt 


dS ele, 
de> Rae 
und damit ist 
M, 
Sate Pag tat 
Uber Gleichung (3) bekommt man fiir die Verschiebung 
1 ¢ R 
u=—pp | Mada t pepe te. (23) — 


SchlieBlich liefert die zweite Gleichung (1’) die Streckenlast p, und die dritte die Streckenlast p,. — 


Macht man noch Gebrauch von den Beziehungen (10) und (11), dann erhalt man nach langerer 
Rechnung fiir die Krafte und Momente in dimensionsloser Form 


sha cos ; 

SSS = , > OSM, 

Ip JEG Yo 7 
VE We sin p 

Ip VEG Ce 

MR =z (1 — cos ¢) 

VEG a ec 2 . 
ee eae 


I, JEG =e sin p cos Py — COS QD, 
P 


Q, R® ele a cos cos >» — COS @ 7 
Tp VEG T 26 YP Pb YP ? 

SS Cee) 
I, JEG Cte Ip JEG ” 

pais ies ae | L— cos.g= 7) E i bode 3) ; Crake 

a 5 COS Y, —-> COS } Se 
I, VEG G 2 2G \ 2 2 P Ip VEG 
ee ee 3 Pe Eas 

Ip EG EG nels 4 26 (3 Cos ~ cos g)| : 


Wie diese Gleichungen zeigen, enthalten die Formeln fiir S und p, jeweils die Integrations- 
konstante c,. Fiir p,; wird man keine Randbedingungen vorgeben, da der Verlauf der Pressungen 
zwischen Fiihrungsrohr und Welle hingenommen wird. Die Konstante c, kénnte also durch eine 
Bedingung fiir S bestimmt werden. Naheliegend waren die Forderungen S = 0 an den Stellen 


E. 


rc 
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a = 0 und a =4a,, da aber nur eine Konstante zur Verfiigung steht, kénnen beide Bedingungen 
nicht zugleich erfiillt werden. Wie die vorliegende Rechnung nun zeigt, sind die in erster Linie 
interessierenden Charakteristiken unabhangig von einem speziellen Wert fiir c,; man wird also 
auch keinen besonderen Wert darauf legen, daB fiir S irgend welche Randbedingungen erfiillt sind. 
Vielmehr ist eine Problemstellung, welcher auch technische Bedeutung zukommt, angeboten. Durch 


0° 90° 180° 270° 360° 
Abb. 14. Verlauf der Verdrehung g der Wellenquerschnitte bei einer Welle mit unendlicher Windungszahl. 


B 
0 90° 760° 210° 0° 


Abb. 15. Verlauf des inneren Torsionsmomentes 2 bei einer Welle mit unendlicher Windungszahl. 


entsprechende Lagerung der Welle bei den Endquerschnitten ist es erreichbar, daB die Verschiebung 
u = 0 ist an den Stellen « = 0 und « =a,. Diese Bedingungen fiir u liefern tiber Gleichung (23) 
die nun festgelegten Konstanten c, und c,. Was die Querkrafte Q, und Q,, sowie die Biegungs- 
momente M, und M, betrifft, so nehmen auch sie in den Endquerschnitten, genau wie nun die 
Langskraft S, nicht mehr beeinfluBbare Endwerte an. 

Einen Anhaltspunkt iiber den Charakter von g als Funktion des Argumentes a oder f erhalt 
man durch Heranziehung der explizite vorliegenden Lésung fiir die Welle mit unendlich hoher 
Windungszahl, also fiir den Fall 6, = 00. Gleichung (12) geht namlich mit m, = 0 iiber in 


oF rw pane 
p= 0 x a pa 


und nach Ausfihrung der Integration liegt die Funktion 
y = 4Aarctg (* tg | 


vor. Das Torsionsmoment M ist damit iiber (4) 
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Fiir den speziellen Wert y, = 180° wurden beide Formeln ausgewertet. Abb. 14 zeigt den Verlauf 
von y, Abb. 15 gibt Aufschlu8 iiber die Abhangigkeit des dimensionslosen Torsionsmomentes 


Eine Betrachtung der Kurven iiberzeugt davon, daB die Lésungen fiir y und M der Welle mit 
unendlicher Windungszahl auch als sehr gute Naherungslésungen fiir Wellen mit endlichem Offnungs- 
winkel a, z. B. %) =, vorausgesetzt da dieser nur nicht zu klein ist, angesehen werden kénnen, | 
da y auBerordentlich schnell mit zunehmendem Argument a oder # abklingt. Es besteht die Még- 
lichkeit, simtliche Gleichungen (24) diesem Gedanken entsprechend auszuwerten. 

Im Gegensatz zu der Welle mit unendlicher Windungszahl, steht die fast gerade Welle mit 
kleinem Offnungswinkel a). Mit der bereits schon verwendeten Substitution 


P=O-+¢> dp = de, 
wobei ¢ ein kleiner Winkel ist, nimmt Gleichung (12) die Form 
%ag— 5 


«| = ees | de | 
2G y2 sin Py , : Ve | 


—?, 
an. Die Integration fiihrt unter Verwendung der Beziehungen (20) und (21) und nach Vernach- 
lassigung héherer Potenzen als /? auf die Gleichung fiir den Drehwinkel 


1 : 
P =P. + 5 (8? —2 8 Bo) sing, , 


welche aussagt, daB y parabolisch von ~, nach gy abklingt. Die Krafte und Momente (24) lassen 
sich entsprechend als Funktionen des Anfangswertes y, oder tiber (22) auch M, und des Argumen- 
tes § explizite angeben. Man erhalt mit Verwendung der Reihenentwicklungen fiir die sin und cos- 
Funktionen und bei Beriicksichtigaung von Potenzen bis 


ee = =a (Bo —B) > 
fee =5 G [Leos =o (F*— 2B Ba), 
nyEe = yee he O 


Eo ee : 
+ Ge VF sin pall + c05 9, 2H +} —4 BB.) ——2 


Ps R? sing, 1/E E a 
: = Vg (1 + xe coe rs) |! +3 #28 fo) 


Ip JEG 
EBay es 3 
— G6) G se 9e(gh + 3A). 
Endlich seien fiir die Welle mit beliebigem Offnungswinkel «,, jedoch unter der Voraussetzung, 


daB ~, hinreichend klein ist, die zu den obigen Gleichungen analogen Formeln zusammengestellt. 
Fir kleine Winkel geht (12) tiber in 


“f 
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und die Integration liefert mit (16) 


o =O, — aie . 
Damit haben die Gleichungen (24) die Form 

MR Pa Sin (by —f) 
EVEG (© pres. Oh 
_M,R =e E Gof (Bo —B) 

Fiji Gughen: mal G Bop’ 

MR E_ Goj? (By — B) 

[ice oe eT, 

ROR AEN 

Ip VEG HESS Se W0)F fy os? 

Q, R® ae a (1 zi 4 Sin (Bo pel 2 

Ip VEG ~ V2 26 Cop Bo’ 
S$ __ of /E Gf aR 
Ip VEG A Cof? Bo Ip VEG ’ 


p, R° = ge E a bg a E se = R 
IpVEG — & of? Bo ye (a sis 76) COG p) rae 7, EC’ 


py ee +e, oj (Bo — 8) 
eee ee 


_ wobei wiederum die Reihenentwicklungen fiir die sin- und cos-Funktionen verwendet und Potenzen 


bis yi beriicksichtigt wurden. 

Das in diesem Abschnitt behandelte Problem kann nun als gelést betrachtet werden. Die 
Giltigkeit der gesamten hier aufgestellten Theorie ist natiirlich an die vorausgesetzten Bedingungen 
oder Einschrankungen gebunden. Neben der Forderung a/R <1, mit deren Hilfe die vereinfachte 
Formel (6’) fiir die Biegungsspannung zustande kam, mu die Bedingung (5) beachtet werden. 
Nun folgt aus (11), daB 

/ ae 
a 
O05| aa Vie 2G. 
gilt und damit (5) sicher dann erfiillt ist, wenn das Verhaltnis a/R der Ungleichung 


geniigt. Bei allen technisch in Betracht kommenden Fallen diirfte diese Forderung aber tatsachlich 
erfiillt sein. 


4. Lésung der zweiten Grundaufgabe. Die Einspannung der Welle bei B hat zur Folge, daB an 
der Stelle « =, die Verdrehung gy, = 0 und ein Einspann-Torsionsmoment M, vorhanden ist. 
Uber die Gleichungen (10’) und (4) erhalt man so fiir die Integrationskonstante 


und (10’) geht iiber in 


=z Fee ees 1 — cos (25) 
Mit pg = q, an der Neer des Belastungsmomentes M, bei A ist 
M,R MR. 36 
SS 1 COS P, - ( ) 
Ip VE REG : 
Die Integration der Gleichung (25) fiihrt auf 
9 


‘Cs di 
aa [ape — o 
E jee pa 
i TEC ob cos p 
& 


i 


342 W. Roth: Die tordierte, einfach gekriimmte Welle mit konstanter Kriimmung —_Ingenieur-Archiv 


--V5 i : er) 
M? R? 
poe + 1—cosp 


erfiillt sein. Diese Gleichung in ces mit (26) reicht zur Festlegung der Charakteristik 


Speziell muB die Forderung 


M, = F,(y,) sowie zur Bestimmung des Binepaat: oe eee M, aus. Die Substitution | 
Pies T= ay. 
fihrt (27) tber in 
sels 
«V2 a i ei See (28) 
2G yi— . sinty ” 
79a 


und (27’) wird zu 


Es wurde 


WM Ae 
ee eae REG 
gesetzt. Beachtet man die Definition der Legendreschen elliptischen Integrale erster Gattung, 
dann ist in der gebrauchlichen Schreibweise (28) 


a | oak rie S®, ae (29) 


a5 [eee =k |K(s) — F (ee =k). (29’) 


Beide Gleichungen erlauben eine Auflésung. Man erhalt aus (29) 


FES Neal aa) (30) 


K(k) — +a Cs 


Diese Gleichung liefert den Winkel y, iiber k als Funktion des dimensionslosen Einspann-Torsions- 
momentes 


und (28’) 


yg =u—2am 


und aus (29’) 


Y, =u%—2am 


wobei 
By = % ral 


wieder die Rolle eines Parameters tibernimmt. Uber Gleichung (26) erhalt man schlieBlich das 


dimensionslose Belastungs-Torsionsmoment 
M,R 


Ip VEG 


a 


als Funktion von gy, oder Q2,. Unter Verwendung der Tabellen fiir elliptische Funktionen wurden 
die Zahlenwerte ermittelt. Abb. 16 und 17 zeigen fiir einige Parameter f, den Verlauf des Momen- 
tes Q, als Funktion des Winkels p,. Diese Ghacalrorenken lassen sich fiir hinreichend kleine q, in 
linearisierter Form einfach darstellen. Setzt man 


2 


cosp = 1—& cos Y, = a) 


2 
@ 


as |e =k Pop (28’) | 


a 
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Fa 


U 90° 780° G10 560° 
Abb. 16. Charakteristikenfeld Q, der Wellen mit einseitiger Einspannung. 


0 50° 180° 270° 560° 
Abb. 17. Charakteristikenfeld 2, der Wellen mit einseitiger Einspannung. 
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dann wird Gleichung (27’) zu 


Pa 
Vso = dp 
= Oe Gar s 


wenn noch iiber (26) das Moment Q, durch Q, ausgedriickt wird. Die Integration fihrt auf die 


einfache Beziehung / 
M,R Pa cS 
a = =e ot ’ 31 | 
IE VEG ive 9 Po (GY 
welche mit Hilfe von (26) auch die Form | 
M, IR ae Pa ] (31’) | 


‘IpVEG 2 Sit Bo 
hat. Die Gleichung (31) sagt aus, daB® die Torsions-Federkonstante der Welle im Nullpunkt oder | 
auch praktisch fiir kleine Winkel q, 


ist. 

Die Charakteristik fiir die Welle mit unendlicher Windungszahl erhalt man aus (27’), wenn 
dort M, = 0 gesetzt wird. Unter der Voraussetzung namlich, daB das Belastungs-Torsionsmoment 
M, endlich ist, wird infolge der Kriimmung der Welle bei unendlicher Wellenlange das innere | 
Torsionsmoment M bis zum Wellenende, also bis zur Einspannstelle hin, auf Null abgeklungen sein. 
Die Integration von Gleichung (27’) fiihrt so direkt auf Gleichung (18). Eine Betrachtung von | 
10, 11, 16 und 17 iiberzeugt davon, da die durch (18) beschriebene Kurve die Eigenschaft einer — 
Grenzkurve hat. Oberhalb von ihr liegen die Charakteristiken der bei B eingespannten, unterhalb 
von ihr die Charakteristiken der bei B frei drehbar gelagerten Wellen. Fiir die Wellen, deren Cha- 
rakteristiken durch die Grenzkurve beschrieben werden ist es also gleichgiiltig, ob eine Kinspannung 
bei B vorliegt oder nicht. 

Praktische Bedeutung haben die Wellen mit kleinem Offnungswinkel a). Als Ausgangsgleichung 
dient (27). Ersetzt man dort Q, mit Hiife von (26) durch Q,, und beachtet man, daB die Wurzel 
nach dem binomischen Satz entwickelt werden darf, so liegt nach Abbruch der Reihenentwicklung 
mit dem zweiten Glied die Gleichung 


?q 
: Ve ‘a el i Coe ay 
vor, welche integriert die Form 
2 BY2 — 22 (y.—9) ++ (Ga —9) €08 G, — F (sin y, — sin g) =0 (32) 
hat. Mit m = 0 geht diese Gleichung tiber in 
Q y/2 — Dip, ++ oq 008%, — + sing, = 0. (32’) 


So wie (32) der Gleichung (27) fiir kleine a, entspricht, ist (32’) der Gleichung (27’) als Ersatz zur 
Seite gestellt. Prinzipiell ware eine Auflésung von (32’) nach Q, zwar méglich, doch ist die Cardani- 
sche Formel ein zu schwerfalliges Werkzeug. Eine einfache Uberlegung verhilft zu einer Naherungs- 


lésung, welche an Genauigkeit nichts zu wiinschen tibrig la®t. Mit dem Ansatz 
1 
20 == (0. +4) 
By V2 


ist es namlich méglich, Gleichung (32’) zu erfiillen, wobei vorausgesetzt werden kann, daB 6 <q, 
ist. Es geniigt also, nur lineare Glieder in ¢ zu beriicksichtigen. Dies ist evident; denn mit ¢ = 0 
wiirde obige Gleichung die Charakteristik einer geraden Welle beschreiben. Die Abweichungen 
fiir kleine a) von der Kennlinie der geraden Welle sind aber so gering, da ¢ die Eigenschaft eines 
Korrekturgliedes hat. Man erhalt so aus (32’) 


Ga (sin Pa — Pa C08 Pq) 
as pe ; 


und damit ist 
= oes -- Bo (sin Q, ere %, cos Pa) 
Bo V2 2 Pa 


a 


(33) 


i 
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Entsprechend bekommt man iiber (26) dann auch 
oe Po Bo (Pq Su Pa) / 
ar — 2 . 33 
Bo V2 y2 Pa (3) 
Abb. 18 zeigt den durch die genaue Rechnung gegebenen Verlauf von Q,; die durch (33) gelieferte 


une deckt sich vollstandig mit dieser Kurve. Ebenfalls ist in Abb. 18 die uber die einfache 
ormel 


<— 


b 


Q, = 2, =" (34) 


_gewonnene Approximationsgerade gestrichelt eingezeichnet. Ein Vergleich der Kurven iiberzeugt 

“davon, daB die durch Gleichung (34) gelieferten Charakiedatiken den durch die Praxis gestellten 

Anforderungen auf Genauigkeit Baro: werden; nur in besonderen Fallen wird man auf die Glei- 
_chungen (33) bzw. (33’) Feackoresveir 

Die Ermittlung der pate aden Krafte und Momente kann in derselben Weise geschehen, wie 

sie schon fiir die bei B drehbar gelagerten Welle angegeben wurde. Nach Einfihrung des Argu- 

mentes 


=) ae 

Pome 26 
Berhalt man aus (30) die Verdrehung als Funktion von a oder f. SchlieBlich liegen nach Abbau 
der Gleichungen (1’) und (2’) die Recultate 


ie 
Byec | ec ti. 
i | he : 
Tp JEG = — |) a sing ; 
M,R 141/E 
i, \EG = Ze (1—cosq), 
Q, R? 
1, VEG = ee ee 
OC, M? R? 
i, am — 003 9) |e +1—cosp , (35) 
= SR? Va: ¢, R? 
pee (ft — cos —_— , 
miEC Ve ( Ph Ip VEG 
oe Ve . M? R? asks: c, R? 
I, EC == rae cos ~) +— ae + He cos p ++ cos p —-~> Cos || ee 
P2Re cE. ER cb 
Ip JEG = ai sing|(1 +3 ges g)— ae (eee — cosy) 


vor. Gleichung (23) beschreibt wieder die Verschiebung u und die Bestimmung der Konstanten c, 
und c, wird so erfolgen, daB u = 0 ist an den Stellen a = 0 und « = a. 
Fiir Wellen mit kleinem Offnungswinkel a, dient Gleichung (32) zur Ermittlung der Verdrehung 
als Funktion von 8. Obwohl diese Gleichung in ~ transzendent ist, gelingt eine Auflésung. Der 
_ Ansatz 
(Po — B) 


) 
Disa ears te 
erfiillt sicher (32) unter der Voraussetzung, dab 
(Po — B) 
<P 


ist, so daB also bereits quadratische Glieder in € vernachlassigt werden kénnen. Der Term & hat 
auch hier die Eigenschaft eines Korrekturgliedes. Bei gerader Welle ware nimlich 

e=9 8", (36) 
und der durch é gelieferte Beitrag kann nur ein geringer Prozentsatz von dem durch Formel (36) 
gelieferten Wert fiir y ausmachen. Nach langerer Rechnung bekommt man so, nachdem noch 
fiir Q, der Ausdruck (33) verwendet wurde, aus (32) 


ee eet oP 
See he eC |) 
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und damit wird 


Poe 4. Po. (poem By cae (37) 


Oe ae Reena sin | Bo 

Setzt man nun die Ausdriicke (33’) und (37) in die Gleichungen (35) ein, so erhalt man schlieB- 
lich explizite Formeln fiir die Krafte und Momente, welche in Abhangigkeit von § und des Anfangs- 
wertes gy, beschrieben werden. Nun haben diese Formeln allerdings eine bereits schon weniger 


A 2, 
90 


60 


40 


60 


5,0 


40 


30 


20 


10 


fa 
0 90° 780° 270° 560° 
Abb. 18. Exakter Verlauf und Niherung fiir 2, bei kleinem Parameter fy. 
gefalligen Aufbau, und man wird vielmehr die Ausdriicke (34) und (36) in die Gleichungen (35) 


einsetzen und nur in besonderen Fallen die Korrekturen€ und  beriicksichtigen. Diese so erhaltenen 
Gleichungen haben dann bei sinnvoller Vernachlassigung positiver Potenzen in f, die Gestalt 


DLR) Seo, 

IpVEG By 2’ 

See ae 

He = Veh ale 

IEC "Air tea): 

em Heft Boon lee] - 
pee Tale Pees, tae 

pias ~ Hae GoM) 85, 


Pp, R° ee oe Bo — —p 
[, JECl°s mg 4c i Po 


Vee 
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und sind zur naherungsweisen Beschreibung des Spannungszustandes in der Welle tatsachlich gut 
geeignet. Weiter unten wird gezeigt, daB die Gleichungen (35’) auch Giiltigkeit fiir Wellen mit 
beliebigem, endlichem Offnungswinkel «) haben; das Belastungs-Torsionsmoment M, mu8 dann 
_aber theoretisch unendlich grof oder praktisch sehr groB sein. 
Auch fiir Wellen mit beliebigem Offnungswinkel a, lassen sich die Krafte und Momente (35) 
als Funktionen von f und ¢, angeben. Es muf jedoch vorausgesetzt werden, da die Verdrehung ¢, 
-hinreichend klein ist und damit auch das Belastungs-Torsionsmoment M,. Gleichung (27) hat fiir 
kleine Winkel p die Form 


4 ee | en ae 
26 — J arte 
” 
Substituiert man fiir (2, den durch (31’) gegebenen Ausdruck, dann liegt nach der Integration die 
-Gleichung 
Sin (By) — 
P=%a a Be e (38) 


vor, welche beschreibt, wie die Verdrehung von ~, auf Null abklingt. Die Formeln (35) haben nun 
mit (38) und (31’) die Gestalt 


MR _ 9 &0j (Po —B8) 
Ip VEG eG y2 Sinf ’ 

M,R ga}/E GSin(fy— B) 
Pyne 2s 2G Sinpy  ” 
"M,R __ 3 1/E Gin? (B — B) 
Tyee 4/6 GiwA 

ZR gf E Sin2G6,—8) 


I,JEG 2 #6 Ginep, 
Q, Rt ae Pa ‘ Le eae 


PEG yon 26) er Cinpe- 

S R? g2 )/E Sint(B—f) , aF 
enc. wpe sterG, Sue Tyre ’ 

pk gt V2 Sin? (By — B) i E Goj2@—8)]_ 4 R? 
Pech eal GC Sin? Bo G Sin? Bo Ty EG. - 


PR _ .)/E(, 4 E\ Sin(bo—P) 
TEC © 2:)06 26] Sin By- ’ 
wenn wieder bei der Aufstellung dieser Gleichungen fiir die sin- und cos-Funktionen die bekannten 
Potenzreihenentwicklungen verwendet und Potenzen bis yj beriicksichtigt werden. 
Wahrend die bei B drehbar gelagerten Wellen die Eigenschaft haben, dafs das innere Torsions- 
moment 
aM Be 
Ip VEG 
ist, wie aus der ersten Gleichung (24) folgt, kann eine solche Schranke fiir die bei B eingespannte 
Welle nicht angegeben werden. Aus (26) ist namlich zunachst ersichtlich, dais fir M,—> co, M,=M, 
ist. Weiter wird bei einer Welle mit beliebigem aber endlichem Offnungswinkel a, das Belastungs- 
Torsionsmoment M, von einer gewissen Stelle y, = ~,; an monoton mit wachsendem —y, zunehmen, 


so daB also fiir groBe Momente M, bzw. M, die Gleichung (27’) sich vereinfacht auf 


<2 


was 


348 W. Roth: Die tordierte, einfach gekriimmte Welle mit konstanter Krimmung Ingenieur-Archiv 


und aus ihr erhalt man mit Hilfe von (39) 

p= pe | (40) 
Bo 

Ein Vergleich der Formeln (34) und (36) mit den Formeln (39) und (40) beweist véllige Uber- 
einstimmung. Daraus folgt aber, wie schon oben behauptet wurde, daB die Gleichungen (35’) 
zustindig sind fiir Wellen mit sehr kleinem Offnungswinkel «) und fir Wellen mit beliebigem Off- 
nungswinkel; jedoch mu im zweiten Fall das Belastungs-Torsionsmoment M, auch hinreichend 
eroB sein. Die durch Gleichung (39) beschriebenen Kurven sind die Asymptoten der Charakte- 


ristiken der bei B eingespannten Wellen. Das Verhalten der Asymptoten aber zeigt sofort, daf die — 


wesentliche Ungleichung (5) nur bedingt erfiillt ist. Unter Verwendung von (25) besteht (5) sicher 


dann, wenn 
E 1/M} R? |R 
Ve Vitec t2<o-} 


ist. Fiir sehr grofBe Momente M, kann diese Ungleichung auch in der einfacheren Form 


ARE | eee [= 
I, VEG a Be 
oder auch mit Verwendung von (39) in der Gestalt 
R 
Pa < . =f Xo 


angeschrieben werden. Theoretisch ist es also méglich, bei vorgegebenen maximalen M,,,,, und 
Pamax ein Verhaltnis R/a so auszulegen, daB die Ungleichung (5) erfillt ist; bei einer aber einmal 
ausgefiihrten Welle diirfen M, und damit auch ¢, nicht zu groB werden, da sonst die oben erwahnte 
Voraussetzung unbeachtet bleibt. 


5. Zusammenfassung. Die tordierte, einfach gekriimmte Welle mit konstanter Kriimmung 
wurde hinsichtlich ihrer elastischen Eigenschaften streng untersucht. Der Kompliziertheit des 
Problems entsprechend, konnten so die Charakteristiken zwar nicht explizit, aber doch exakt 


Sdq 


fal 


_— 
0° 180° I60° 540° 720° 900° 1080° 
Abb. 19. Aufteilung der Halbebene durch die Charakteristikenfelder 2g. 


angegeben werden. Fiir spezielle Falle gelang es, teils einfache Genaue-, teils einfache Naherungs- 
formeln herzuleiten. Der wesentliche Inhalt der Arbeit gipfelt neben der Angabe der Kennlinien 
in der Erkenntnis, daB zwischen beiden untersuchten Grundproblemen ein Zusammenhang besteht, 
der am treffendsten in dem Verhalten der Charakteristiken erkannt wird. Abb. 19 mag zusammen- 
fassend die Situation beschreiben. Fiir den Parameter Bo = sind dort die Charakteristiken der 
bei B eingespannten und drehbar gelagerten Wellen eingezeichnet. Beide Kennlinien liegen jeweils 
in Feldern, welche durch die Grenzkurve (18) voneinander getrennt werden. Oberhalb der Grenz- 
kurve befinden sich die Charakteristiken der bei B eingespannten Wellen; diese Kurven gehen mit 


Bee De» 
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 Y,-> 00 in ihre Asymptoten (39) tiber. Die Q,-Achse ist identisch mit der Charakteristik der 


Wellen, denen der Parameter f, = 0 zugeordnet ist oder den geraden Wellen. Unterhalb der 
Grenzkurve liegen die Charakteristiken der bei B drehbar gelagerten Wellen; sie haben die Periode 
22, und die @,-Achse fallt zusammen mit der Charakteristik der geraden Wellen mit dem Para- 
meter fy = 0. Die Kennlinien zeigen einen Verlauf, wie er auch nur ahnlich von den bisher be- 
kannten Federn nicht erreicht wird. 

Fiir die in der Welle auftretenden inneren Krafte und Momente, ebenso fiir die von dem Fih- 


-rungsrohr herriihrenden Pressungen, wurden Formeln angegeben. Die fiir die Pressungen zu- 


La, 
360° 
270° 
780° 
90° x 
Ys 
° ° ° ° ry 
0 90 780 210 560 


Abb. 20. Vergleich der MeBwerte mit den durch die Theorie gelieferten Werten. 


standigen Gleichungen haben weniger Bedeutung, da aus konstruktiven Griinden bei einer prak- 
tischen Anwendung der gekriimmten Welle einzelne Lager an Stelle des Fiihrungsrohres vorzu- 
ziehen sind. Die dadurch bewirkte Anderung der Charakteristiken ist aber so gering, da sie auch 
zugleich vernachlassigbar klein ist. So wurde mit Hilfe eines Modells, das die Daten R/a = 250, 
& = 80,50° und 6, = 2/2 hatte, und bei dem die Welle bei A und B drehbar gelagert war, die in 
Abb. 20 eingetragenen MeBwerte ermittelt. Im ganzen war die Welle nur dreimal an den Stellen 
& = 0°, a = 40° und « = 80,50° gelagert. Selbst mit nur zwei Lagern bei A und B war die raum- 
liche Ausbiegung der Welle iiberaus gering und es ergaben sich dieselben MeBergebnisse. Es wurde 
allerdings jeweils dafiir gesorgt, dafB die Endtangenten der Welle bei A und B im unbelasteten 
Zustand nach Aufbringung des Belastungs-Torsionsmomentes M, erhalten blieben. 

Verschiedene Anwendungsmdglichkeiten fiir die gekriimmte Welle liegen auf der Hand. Der 
Verfasser hat sich jedoch in dieser Arbeit bewuBt nur mit der Erforschung der Grundlagen be- 
schaftigt. 


(Eingegangen am 16. Februar 1959.) 


Anschrift des Verfassers: Dr.-Ing. Werner Roth, Karlsruhe, Technische Hochschule, 
Institut fiir Technische Mechanik. 
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Eine Erweiterung des Maxwellschen Vertauschungssatzes 
Von Mok-Kong Shen 


Unter Voraussetzung der Giltigkeit des Superpositionsprinzips lautet der Maxwellsche Ver- 
tauschungssatz folgendermaBen (Abb. 1): Die Komponente der Verschiebung des Punktes 2 eines 
starr gestiitzten Kérpers, verursacht durch die Kraft P, in der Richtung der Kraft P,, ist von glei- 
cher Grié&e wie die Komponente der Verschicbung des Punktes J, verursacht durch die Kraft P, 
in der Richtung der Kraft P,, wenn P, und P, gleich grof sind. 

Dieser Satz, der von Maxwell aus dem Arbeitsatz abgeleitet worden ist, zeigt nicht die Beziehung 
zwischen den beiden Verschiebungen selbst, sondern nur die, welche zwischen den zwei Kompo- 
nenten der Verschiebungen besteht. Es ist der Zweck der vorliegen- 
den Abhandlung, die Vertauschbarkeit der Verschiebungen weiter zu 
untersuchen. 

Unter Voraussetzung der Giltigkeit des Superpositionsprinzips 
ergeben sich fiir zwei beliebige Krafte mit den Vektoren P, und 
P, an den Punkten I und 2 die folgenden Vektorgleichungen: 


64 = P,-@, (1) 
Abb. 1. O12 aA P, PD > (2) 


worin 6, der Vektor der Verschiebung des Punktes 2 verursacht durch P, ist, ebenso 6,, der Vektor 
der Verschiebung des Punktes 1 verursacht durch P,, und ®, ’ die zwei Dyaden, welche Funktionen 
der Lage der Punkte J und 2 und unabhangig von den Kraften P, und P, sind. 

Aus dem Maxwellschen Vertauschungssatz (oder als Sonderfall des Bettischen Satzes) erhalten 
wir 
On P, = 012° Py, (3) 
woraus gemaB (1) und (2) folgt 


PD Pepe sere (4) 


Da P, und P, ganz beliebige Krafte sein kénnen, muf ®’ die konjugierte Dyade von @ sein. Schrei- 
ben wir @ in ihrer Normalform 
OP=aii+bjj +ckk, (5) 
so ist 
DP =avitbjjt+ekk. (6) 
Damit ergeben sich die Verschiebungen 6,, und 6,., wenn P, und P, die GréBe eins annehmen 
und in Richtungen der Achsen 1, 7, k und 1’, j’, k’ liegen, zu 
LoD Ot ss jo P= o7 > ka O= coke (7) 
LOOSE I fee eN he ko DO = chan (8) 
Die Achsen i, j, k und v’, j’, k’, die zwei Systeme von je drei zu einander senkrechten Richtungen 


bilden, werden im folgenden Hauptvertauschungsrichtungen genannt. 


Ist P, durch den Versor 
Qsit +jj +k 
mit P, verbunden: 
Pye Pa Oe (9) 
so erhalten wir gemaB (1) und (2) 
dy = P,- B= (P,-Q7)-(Q- @’. Q) 
eh, On Oe) ee (10) 
Fiir zwei ganz beliebige Krafte 
Pi=at+Bjtyk, 
Prawi + iy ty k 
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sind die Verschiebungen 
dy = Py-P®=aai +Bbj’+yck', 
dy. = P,-D =o ai+Pp’ bj+y'ck 


und also 
On Po =e6o a+ Bh b+yy'c=06,,- P,, 
_was zu Gleichung (3) zuriickfiihrt. 

Diese Ergebnisse kénnen folgendermafen zusammengefaBt werden: Es gibt, unter Voraus- 
setzung der Giiltigkeit des Superpositionsprinzips, fiir die Verschiebungen 6,, und 6, zweier Punkte 
an einem starr gestiitzten Kérper, verursacht durch die Krafte P, und P,, zwei entsprechende, aus 

je drei zueinander senkrechten Achsen bestehende Hauptvertauschungsrichtungssysteme, welche 
die folgenden Eigenschaften haben: 

1. Die durch zwei gleich gro®e, diesen Systemen ahnlich zugeordnete Krafte P, und P, ver- 
ursachten Verschiebungen 0,, und 0,, sind immer denselben Systemen ahnlich zugeordnet und von 
gleicher GréBe. 

2. Die durch zwei gleich groBe, in entsprechenden Hauptvertauschungsrichtungen liegende 
Krafte P, und P, verursachten Verschiebungen 6,, und 6,, liegen in denselben Richtungen und 
sind von gleicher GréBe. 

Es ist verstandlich, das unter Umstinden die Dyade @ sich in mehreren Normalformen aus- 
driicken lassen kann. In solchem Falle werden wir entsprechend mehr Hauptvertauschungsrich- 
tungssysteme haben. Andererseits kénnen die Werte a, b, c in (5) gleich Null sein, und wir werden 

dann Richtungen haben, in denen keine Verschiebungen 0,,, 6,. entstehen kénnen. 

Alle diese Gleichungen sind, wie der originale Maxwellsche Satz, auch giiltig, wenn wir, statt 
Kraft und Verschiebung, Moment und Drehwinkel an einem Punkte oder an beiden Punkten 
haben. Im allgemeinen entstehen daher die folgenden vier Gruppen von Hauptvertauschungs- 
richtungsystemen fiir zwei bestimmte Punkte J und 2: 


Gruppe | am Punkte 1 | am Punkte 2 


I Kraft und Verschiebung Kraft und Verschiebung 

II Kraft und Verschiebung Moment und Drehwinkel 
Ill Moment und Drehwinkel Kraft und Verschiebung 
IV Moment und Drehwinkel Moment und Drehwinkel 


Wenn wir den Punkt 2 mit dem Punkte 1 zusammenfallen lassen, so wird die Gruppe III gleich 
der Gruppe II sein, und in der Gruppe I bzw. in Gruppe IV werden die zwei Hauptrichtungsysteme 
miteinander zusammenfallen. [Gleichung (4) wird hier P,-®- P, = P,-@- P,, so daB @ eine 
symmetrische Dyade sein mu8.] Damit kommen wir zu der bekannten Erkenntnis zuriick: 

Es gibt, unter Voraussetzung der Giiltigkeit des Superpositionsprinzips, fiir die Verschiebung 61, 
des Punktes I eines starr gestiitzten Koérpers, verursacht durch die Kraft P, an demselben Punkt, 
drei zueinander senkrechte Hauptrichtungen, in denen P, und 6,, zusammenfallen. 


(Eingegangen am 16. Februar 1959.) 
Anschrift des Verfassers: Mck-Kong Shen, Karlsruhe, WiesbadenerstraBe 27, bei Willhauck. 
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Die Extrema der Ubersetzungen 
in ebenen und spharischen Kurbeltrieben 


Von W. Meyer zur Capellen 


1. Einleitung. In einer fritheren Arbeit! wurden die Extrema der Geschwindigkeiten in Kurbel- 
trieben und damit indirekt die Extrema der Ubersetzungen in diesen untersucht. Es lag nahe, 
die Frage nach den Extrema auch auf raumliche Getriebe, vor allem auf spharische Getriebe aus- 
zudehnen. Da die ebenen Kurbelgetriebe als Sonderfalle der spharischen Getriebe aufgefaBt werden 
kénnen, soll zunachst noch kurz auf die Ubersetzung in ebenen periodischen Getrieben ein- und in 
Anlehnung an diese Entwicklungen auf die spharischen Getriebe iibergegangen werden. 


2. Ebene Getriebe. a) Begriff des Ubersetzungsverhaltnisses. Bei ,,stetigen“’ Getrieben, 
z. B. bei den Radertrieben, wird das Ubersetzungsverhaltnis i definiert als das Verhaltnis zwischen 
Antriebswinkelgeschwindigkeit (oder -Drehzahl) und Abtriebswinkelgeschwindigkeit (oder -Dreh- 
zahl) i = Wantrich ? abtriess Dieses Verhaltnis ist bei stetigen Getrieben konstant, bei periodischen 
Getrieben aber nicht. Da zudem bei diesen, z. B. bei der Kurbelschwinge, die Winkelgeschwindig- 
keit im Abtrieb Null werden kann, also dann i nach unendlich gehen wiirde, pflegt man bei den — 
periodischen Getrieben den reziproken Wert als Ubersetzung zu definieren: Up = Wabtried? antrich> | 
wobei der Index p auf das periodische Getriebe hinweisen soll. Im folgenden wird nur von periodi- 
schen Getrieben die Rede sein, so daB der Index p nicht mitgeschleppt werden soll. 


Abb, 1. Aufbau des allgemeinen Getriebes. 


Im iibrigen mui man noch unterscheiden zwischen Abtrieb mit Dreh- und Abtrieb mit Schub- 
bewegung. ; 

1. Drehbewegung im Abtrieb. Durch einen Zwischenmechanismus werde der (gleichmaBige) 
Antrieb des Gliedes 1 (Abb. 1), mit der Winkelgeschwindigkeit jw, in eine periodische Drehbewegung | 
des Gliedes k. mit der Winkelgeschwindigkeit «, verwandelt. Dann gilt i = ,:@;, aber nach 
bekannten Satzen iiber die Drehung dreier Ebenen gegeneinander, hier der Ebene 1, der Ebene k | 
und der ruhenden Ebene 0 (Gestell), auch | 


3a Oia eae P 
oy, q Dives bg (1) 


worin d = Ay) M, = P,y P,, die konstante Lange des Steges, P,, der Relativpol der Drehung von 
k gegen I und p die Strecke p = P,, P,, = P,, A, und q die Strecke q = Pri Pro = Puy My be- 
deuten. Die Strecken kénnen auch negativ sein, so wie bei der Kurbelschwinge in Abb. 2, und wenn 


* W. Meyer zur Capellen, Ing.-Arch. 25 (1957), S. 140. 


Aa 


“a 


bs 
a. 
- 
: 
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man dann p durch — p ersetzen wiirde, so ist auch 


~ p+d? 2) 


worauf spiter zuriickgegriffen wird. 
Ist A ein Punkt des Antriebsgliedes, von dem aus die Bewegung weitergeleitet wird (wie z. B. 
bei der Kurbelschwinge in Abb. 2) und K der Punkt des Gliedes k, auf welchen durch den Zwischen- 


mechanismus die Bewegung tibertragen wird (z. B. Punkt B in Abb. 2), so gilt auch mit A, A = a 


Reo 


Abb. 2. Relative Polbahntangente beim Kurbeltrieb. 


und M, K =k sowie vy und vg als den Geschwindigkeiten dieser Punkte zunachst vy = a@,, 
Ux =k «a, und damit auch 


- _Ug/k _ a vK . 
aad ewe (3) 


d.h. i, dem Verhaltnis vg: v4 proportional, kann somit aus diesem ermittelt werden: Wenn vy 
bei den iiblichen Geschwindigkeitskonstruktionen gleich a cm gezeichnet wird, so zeigt sich i = vg,/k 
als Verhaltnis zweier Strecken, wobei vx, die gezeichnete Lange der Geschwindigkeit vg bedeutet. 

2. Schubbewegungim Abtrieb. Vollfiihrt das Abtriebsglied eine Schubbewegung, wie z. B. 
beim Schubkurbeltrieb nach Abb. 3 oder bei dem sechsgliedrigen Getriebe nach Abb. 4, so kénnen 
die Formeln (3) oder i = @,/m, nicht mehr benutzt werden, da dann @, verschwindet oder auch 
k nach unendlich geht. Es wiirde sich fiir i der unbrauchbare Wert Null ergeben. Man kann dann 


“davon ausgehen, da die Geschwindigkeit v4 in die Geschwindigkeit v, des Gliedes k verwandelt 


wird, und setzt i als das Verhaltnis 


eel 3 

int (4) 
an. Bei dem Schubkurbeltrieb in Abb. 3 erkennt man aus dem Dreieck der gedrehten Geschwindig- 
keiten (Satz von Euler), daB hier 


Nee Be 2) Ps 5 
j=BaF (5) 


gilt, wie im iibrigen auch Abb. 4 zeigt, welches das sechsgliedrige Getriebe in der Stellung zeigt, 
in welcher i ein Extremum hat. Hier bildet E den Abtrieb, und es gilt i = v;/v4 = vg/v4, sowie 
p =A, Ps, (und gleich v,; = vg, da v4 gleich der Kurbellange A, A gezeichnet wurde). 

Wenn ein Kurventrieb eine drehende Bewegung in eine Schubbewegung verwandelt (Abb. 5), 
so kann — im Gegensatz zu den soeben geschilderten Gelenkgetrieben — nicht von einer kon- 
stanten Antriebsgeschwindigkeit gesprochen werden. Es ware nach der Definition (4) und (5) 
mit der (verinderlichen) Fiihrungsgeschwindigkeit vp = rw, dann i = p/r; aber diese Definition 
befriedigt nicht, denn das Ubersetzungsverhaltnis soll nach (3) doch der Abtriebsgeschwindigkeit 
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proportional sein. Es empfiehlt sich daher die Geschwindigkeit v4 hier im Abstande des Grund- 
kreishalbmessers R, zu definieren als v4 = Ry @,, so daB dann 1 wiederum dimensionslos als 


i= (6) 


Ry a, 


erscheint. Wesentlich diirfte es sein, daB i auf jeden Fall eine dimensionslose GréBe wird und der 
Abtriebswinkelgeschwindigkeit oder der Geschwindigkeit des Schubgliedes proportional ist. 


ts 7 


Abb. 3. Relative Polbahntangente beim Schubkurbeltrieb. 


b) Das Extremum des Ubersetzungsverhidltnisses. 1. Formale Bedingung. Die in 
dem Wert (1) von i auftretenden Strecken p und q sind Funktionen des Antriebswinkels « aus 
Abb. 1, und wenn Striche Ableitungen nach « bezeichnen, so wird gemaB (1) 


i = 76) = — (7) 
oder gemaB (5) 
ek toe (8) 
Damit hat aber i in beiden Fallen ein Extremum, wenn 


P = p(s) =0 (9) 
wird. Die Berechnung von p und p’ fihrt auf langere Ausdriicke, so daB z. B. fiir die geschrankte 
Schubkurbel eine Gleichung sechsten Grades folgt!. 

2. Geometrische Deutung. Betrachtet man die Bewegung des Gliedes k gegeniiber Glied 1 
(Abb. 1), so kann diese auch aufgefaBt werden als das Abrollen der mit dem Glied k verbundenen 
Polkurve K,, Z auf der mit dem Glied 1 verbundenen Polkurve K,,,, d. h. als das Abrollen der relativen 
Polkurven K,,, und K,,, aufeinander. Sie beriihren sich im Relativpol P,, langs der gemeinsamen 


LU Meyer zur Capellen u. Mitarbeiter, Bewe haltni a 
: ; gungsverhaltnisse an der hrankt 
Forschungsbericht des Wirtschaftsministeriums von Nordrhein-Westfalen Nr. 146° Opleden- Kak 10S 
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Abb. 4. Sechsgliedriges Getriebe mit Schubbewegung im Abtrieb fiir i = 


Abb. 5. Kurventrieb, 


relativen) Polbahntangente t,,, = t,,. Denkt man sich das Antriebsglied I festgehalten — gekenn- 
zeichnet durch die Achse ¢, (Abb. 1 und Abb. 6) — und dreht den Steg, das Glied Null, um den 
Winkel « um A), so ist K,, die (relative) Rastpolkurve und stellt 4) P,,; = p(a) die Gleichung 
dieser Polkurve in Polarkoordinaten dar (Abb. 7), wobei p den Fahrstrahl und « das Argument 


TUnp| 


Abb. 6. Polkuryven und Subnormale p’. 
1 Dubbels Taschenbuch. Bd. I, S. 114. 11. Aufl. 


bedeutet. Nun stellt aber p’ = p’(«) = 
dp/dx bekanntlich die Polarsubnor- 
male dar!, d.h. die Strecke p’ wird 
auf der zum Fahrstrahl A, P,, in A) 
senkrecht stehenden Geraden durch die 
Normale, also hier durch die (relative) 
Polbahnnormale n,,. abgeschnitten; in 
Abb. 6 ist A, N = p’. 

Wenn aber p’ = 0 sein soll, so muB 
die Polbahnnormale n,,. in Richtung 
A, P,1, d. h. in Richtung des Steges 
fallen, oder es muB die relative 
Polbahntangente t,, senkrecht 
zum Steg stehen, wenn das Uber- 
setzungsverhaltniseinKxtremum 
haben soll. 


g, 
<> 


Abb. 7. p = p(a) in Polarkoordinaten, 


Berlin 1958. 
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Im iibrigen kann die Bedingung p’ = 0 auch kinematisch-geometrisch noch anders gedeutet 
werden: Fiir die auf dem Polstrahl gelegenen Punkt 4A, und M, (Abb. 6) gilt nach der Euler- 


Savaryschen Formel 


ita, (10) 
r its w 
da A, der Kriimmungsmittelpunkt der relativen Bahn von M, gegeniiber Glied I ist, mit 
r= —(p—d), 1 =—p, w= Dsing und sing = p/n,,, gemaB Abb. 6, sowie n,,) = /p? + p? 
fiir den Wendekreisdurchmesser D die einfache Formel 
—d 
D= p d Nel + (11) 


Da nun D’ = D’(a) proportional p’ ist, wie eine einfache Rechnung zeigt, so bedeutet die Bedingung 
p’ =0 auch D’ = 0, d. h. aber einen stationéren Wert des Wendekreisdurchmessers. 

Dies konnte auch daraus geschlossen werden, da der Punkt M, der (relativen) Kreispunkt- 
kurve angehort, aber auch der ja jetzt auf den Steg fallenden Polbahnnormalen, und da damit 
jene in diese Normale (und noch einen Kreis) entartet. Dieser Sonderfall der Kreispunktkurve 
bedeutet aber! einen stationaren Wert von D, fiir den dann unmittelbar mit n,,, =p auch 
D = p(p—d)/d folgt. 

Eine rein bewegungsgeometrische Deutung liefert schlieBlich noch das Folgende: Der Punkt P,, 
wandert auf dem jetzt als fest betrachteten Steg mit der Geschwindigkeit dp/dt = (dp/d«) (da/dt) 
= p’ @,, d.h. diese (beilaufig eine Komponente der relativen Polwechselgeschwindigkeit) ist pro- 
portional p’. LaBt sich nun die Wanderung des Punktes P,, als tatsichliche Bewegung deuten, so 
mu im betrachteten Augenblick die Geschwindigkeit dp/dt gleich Null sein (vgl. nachsten Absatz 3). 

3. Anwendung auf die Viergelenkkette. Die in Abb. 2 ausgefiihrte bekannte Konstruktion 
der (relativen) Polbahntangente bei der Kurbelschwinge zeigt, daB fiir das Extremum der Winkel 
tre, Ps; By) = Winkel P,) P;,; B = 93, gleich einem Rechten sein muB, d. h. daB die Kollineationsachse 
senkrecht zur Koppelmittellinie AB stehen mu, wie ja bereits fiir die Viergelenkkette friiher 
gefunden.? 

Man kann im iibrigen ohne Kenntnis des Begriffs ,,Kollineationsachse“ diese Bedingung finden: 
In Abb. 8 tragt die Koppel eine Schleife, diese paart sich mit einem Gleitstein, der um einen in 
Richtung des Steges verschiebbaren Schieber drehbar ist. Dieser Drehpunkt S fallt mit dem Rela- 


tivpol P;, zusammen und der oben genannte Wert p’ wird jetzt proportional seiner Geschwindig- _ 


keit v. Diese kann in bekannter Weise aus b = b,,, 
verschwinden, so muB P,,S = Py) P3, auf der Schleifenmittellinie, d. h. auf 4B senkrecht stehen’. 
Der Wert 7’, wie er in (7) angegeben ist, bedeutet auch die bezogene Winkelbeschleunigung des 


Roy 


Abb. 8. Getriebliche Deutung der Geschwindigkeit p’. 


1 W. Meyer zur Capellen, Ing.-Arch. 17 (1949), S. 308. 


2 Vel. FuBnote 1 auf S. 352 und F’, Freudenstein, Trans. ASME (1956) S. 779 — ASME Paper 55 — SA-20. | 
3 H. Brandenberger, Kinematische Getriebemodelle, Ziirich 1955, wies iibrigens auf das Extremum der 


Strecke p als Kennzeichen fiir ein Extremum des Ubersetzungsverhaltnisses hin, ohne jedoch diese Feststellung 
irgendwie auszuwerten. 7 


= by +0, ermittelt werden. Soll nun aber p’ | 


q 
| 


vn. Band 1960 W. Meyer zur Capellen: Ubersetzungen in ebenen und spharischen Kurbeltrieben BEY 


Gliedes k, d. h. es ist 


bei drehendem Abtrieb oder 


/ 


Sees by, P 
by a 


bei Schubbewegung im Abtrieb mit b,, b4 als den Beschleunigungen des Gliedes k und des Punktes A. 
Fihrt man beilaufig in (7) den Winkel ¢,,; aus Abb. 6 durch p’ = p ctg y,, ein und beachtet, 
daB nach (1) auch i = 1 + d/(p — d), also d/(p — d) = i —1 ist, so kann man auch 


i =i (1 —i) ctg gy (7*) 
schreiben, d.h. erhalt die bereits bei Freudenstein fiir die Viergelenkkette angegebene Beziehung. 


F Die Gleichung (7) kann nun zur zeichnerischen oder zeichnerisch-rechnerischen Ermittlung der 
s ee ; 9 ; ; 
ezogenen Winkelbeschleunigung ¢3/w} herangezogen werden, wenn der Relativpol noch erreichbar 


Abb. 9. Konstruktion von tt’, Abb. 10. Kopplung zweier Viergelenkketten. 


und damit die relative Polbahntangente einschlieBlich der Polarsubnormalen leicht angegeben 
werden kénnen. In Abb. 9 gilt z. B. nach eingetragener Konstruktion 


“7 p’ d dtgp MGC) 


re ee) pe 

Besonders einfach wird nach (8) die Konstruktion von i’ und damit hier des Verhaltnisses 
bp/b4 = b3/b4 bei der Schubkurbel in Abb. 3: Es wird i’ und damit bp, proportional der Strecke p’, 
bei dem benutzten Mafistab aber gleich bz. 

c) Kopplung zweier Viergelenkketten. Bei Kopplung zweier Kurbelschwingen oder 
allgemeiner zweier Viergelenkketten (Abb. 10) kann man von der Multiplikation der Ubersetzungs- 
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verhaltnisse ausgehen: Es ist 


7) Ws S gat Ae ; 
oe eye tet) lee (12) 


109 M39 10 


wenn 


(13) 


Abb. lla. Relative Polbahntangenten bei 
Kopplung zweier Viergelenkketten. 


Abb. 11b. Konstruktion der Winkel 2, und ju. 


gesetzt wird. Bildet man jetzt i’ = di/dx unter Beachtung, daB i, zunadchst eine Funktion von f 
ist, so folgt i’ = ij(a) i, + i, 12(B) B’(x) , worin B’(~) = i, ist, und man hat nach einigen Rechnungen 


7 “> (Py dy Py d, 
iis. a : 14 
: (5, Pi ar Ps a ( ) 


Werden jetzt noch die Winkel ~, und @, der beiden relativen Polbahntangenten t,,,. = ts, und 
tse = !53 eingefiihrt (Abb. lla), so kann man die Bedingung i’ = 0 ein fiir Extremum in der 
Form 


d, ctg Q, af d, Cte Py 0 
P, 1 


(15a) 
oder 
C2 Ug tS 0. (15b) 


schreiben, worin “4, und “1, die in Abb. lla konstruierten Winkel darstellen. Greift man dann die 
Strecken y, = / tg uw, und y, = / tg uw, mit beliebiger MaBstabseinheit J ab, Abb. 11b, so trage man 


die Summe y = y, + y, = f(a) in der Nahe des zu erwartenden Extremums als Kurve iiber « auf 
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und bestimme ihren Schnittpunkt mit der a-Achse (vgl. auch Absatz d) oder bringe die Kurven y, und 
—Y2 zum Schnitt, um die Stellung des Extremums zu gewinnen.! 
Wird die Bedingung aus (15a) in der Form 


PitsPr GotS Po 
Tele Secor eae ee) 


geschrieben, so lautet sie auch mit den in Abb. 12a, b konstrvierten Winkeln y, und v, abnlich 


tg%, +tgv,=—0 oder 4+2,=0. (15d) 


Abb. 12 a,b. Konstruktion der Winkel », und », bei der Kopplung zweier Viergelenkketten. 


Beilaufig laBt sich schlieBlich noch der Bedingung i’ = 0 auch die Form 


1 \ : 
ie a 1) etg P = (1, — 1) ctg p, (16) 
geben. 
Ist das Abtriebsglied 5 ein Schubglied, so gilt 
: U5 V5 U > Ww bas oe 
j= >= OS Lae a ae (17) 


| 

mit B, B = b, is, = 1, = p,/q, und i, = i, = po/b. Die gleichen Uberlegungen wie vorher fihren 
dann zu der Bedingung 

Piteg 


re + tg go, =0 (18a) 


oder auch 
(1/t; — 1) ctg gp, + ctg gy, = 0. (18b) 
d) Ermittlung der Extremalstellung. Will man bei vorgegebenem Getriebe die Stellung 
fiir das Extremum der Ubersetzung ermitteln, so kann man sich der direkten Methoden bedienen, 
d. h. man ermittelt in der Nahe des zu erwartenden Extremums die Winkel-, bzw. die Schub- 
beschleunigung des Abtriebs und tragt sie tiber dem Kurbelwinkel « auf, und der Schnitt der so 
erhaltenen Kurve mit der «-Achse liefert den gesuchten Wert a. Bei den indirekten Methoden 
kann man in dhnlicher Weise als Funktion von « gewisse Strecken € und 7, wie friiher beschrieben? 
oder die soeben entwickelte Summe y = f(a) oder die Strecke p’ = p(x) benutzen. Bei Ermittlung 
der Nullstelle konnte man naturgemaB durch Iteration Verbesserungen erreichen und zwar unter 
Benutzung der Regula falsi, da ja nur die Funktionswerte f(a) zeichnerisch ermittelt werden kén- 
nen, nicht ihre Ableitungen’. 


1 Bei der Konstruktion der Relativpole treten haufig flache Schnitte auf. Deshalb liefert dieser Weg nicht 
immer die genaue Stellung des Extremums. 

2 Vel. FuBnote 1 von S. 352. 

3 W. Meyer zur Capellen, Ing.-Arch. 27 (1959). S. 53 u. 73. 
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Auch wenn rein rechnerisch vorgegangen wird, so sind die Nullstellen einer gewissen Fnktion f(a) 
zu bestimmen (z. B. der rechnerisch ermittelten Winkelbeschleunigung). Da aber hierbei, von 
Sonderfallen abgesehen, sich Gleichungen héheren Grades in sin « oder cos « ergeben, so wird man 
nur bei einer Gleichung zweiten Grades die Lésung auch zeichnerisch finden kénnen, bei Glei- 
chungen dritten Grades und vierten Grades nur rechnerisch, im allgemeinen aber durch Iterations- 
methoden. In diesem Sinn wird man die Stellungen im allgemeinen nicht mathematisch genau 
angeben kénnen, wenn auch die praktische Genauigkeit der genannten Verfahren durchaus aus- 
reichen diirftet. 


3. Sphirische Getriebe. a) Begriff des spharischen Getriebes. Beim ebenen Gelenk- 
getriebe sind die vier Gelenkachsen zueinander parallel, beim sphirischen Gegenstiick gehen sie 
durch einen Punkt, und denkt man sich diesen als Mittelpunkt einer Kugel vom Halbmesser 
R = ,,eins‘*?, so sind die Mittellinien der Glieder jetzt durch GroBkreise und ihre ,,Langen“ durch 
Bogen auf der Einheitskugel, d. h. durch Winkel gegeben (vgl. Abb. 13%). Wir bezeichnen diese 
Langen der Glieder 1, 2, 3, 4, entsprechend dem ebenen Getriebe mit a, c, b, d, die Drehwinkel der 
Glieder entsprechend mit « und /. 

b) Das Ubersetzungsverhaltnis. Wahrend bei der ebenen Bewegung Gangpolkurve und 
Rastpolkurve, bzw. kérperlich gesehen, die entsprechenden Zylinder aufeinander abrollen, so be- 
kanntlich bei der Bewegung um einen Punkt, d. h. hier bei den sphirischen Getrieben, Gangpolkegel 
und Rastpolkegel. Z. B. dreht sich in Abb. 13 das Glied 2, die Koppel, um die Achse 20, 
welche den Schnitt der Ebenen 10—21 und 30—32 gegeniiber dem Gestell bildet. Das Ubersetzungs- 
verhaltnis i ist hier ebenfalls durch das Verhaltnis der Winkelgeschwindigkeiten «.,) = d{/dt und 
19 = da/dt gegeben. Nun gilt doch fiir die drei Winkelgeschwindigkeiten bei der Bewegung zweier 
K6rper k und i gegeneinander hier ahnlich wie in der Ebene (Abb. 14) 

Oxo = Wig + Wu; - (19a) 
aber auBerdem 


Wy9: Oo = sin a; : sin O,. (19b) 


Wendet man diese Beziehung auf das Gelenkviereck (Abb. 13) an, so gilt mit den Indices i = | 
und k = 3, a; = p als Winkel zwischen den Achsen 10 und 31 sowie 0, = q = p + d als Winkel 


(bb. 13. Spharischer Kurbeltrieb — Erweiterung der Koppel durch 
eine Schleife zur Bestimmung des Extremums des Ubersetzungsverhalt- 
nisses, Abb, 14, Addition der Winkelgeschwindigkeiten, 


' Dies muBte hier betont werden, da bei einem Vortrag des Verfassers vor Ingenieuren tiber die GroBt- 
und Kleinstwerte der Ubersetzungen in Getrieben in der Diskussion bedauert wurde, daf} man die Stellung 
nicht ,,genau‘* finden kénne und diese Ungenauigkeit im Verfahren vermutete und nicht in dem Auftreten 
Gleichungen hoheren Grades. 

2 W. Meyer zur Capellen, VDI-Ber. Bd. 29 (1958) S.91; W. Meyer zur Capellen, Werkstatt und Betrieb 
91 (1958) S. 435; W. Meyer zur Capellen, Werkstatt und Betrieb 91 (1958) S. 723; R. Beyer u. E. Schérner, 
Raumkinetische Grundlagen. Miinchen (1953). 

5 In meiner Vorlesung ,,Raumgetriebe* hat sich die Verwendung einer ,,Kugeltafel‘* als zweckmaBig 


to) 


erwiesen. Die Anregung hierzu gab H. Mohrmann durch sein Buch ,,Nicht euklidische Geometrie“, Leipzig 


1930, fur das ich seinerzeit als Amateurphotograph einige Aufnahmen machte, u. a. auch von einer Kugeltafel. 
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zwischen den Achsen 30 und 31 


pes OE re ee (21a) 
do M9 sin q sin (p + d) 
oder auch 
~ cos d + sindetg p. (21 b) 


c) Besondere Werte des Ubersetzungsverhaltnisses. Bevor wir zu den Extrema iiber- 
gehen, sei noch auf einige besondere Werte von i hingewiesen. 

1. In der Steglage fallt das Glied @ in die Richtung des Stegs, und es ist dann p = 4 fir die 
auBere, p = —a fiir die innere Steglage, d. h. 


1 = tae eh oe 1 7 Ca ters 

= = cos d-+sindctga bzw. yo) e008 d—sindctga. (22 a) 
1 2 

Diese Werte sind, ebenso wie bei dem ebenen Getriebe nur von a und d abhangig, nicht von der 

GréBe der anderen Glieder, wobei das ebene Getriebe mit a/d =A die Werte 


1 1 1 1 
— — f = i 2 
i, 1+ 7 bzw. i, ] 4 (22b) 
liefert. Addiert man die beiden Gleichungen (22a), so ist beilaufig das harmonische Mittel aus i, 
und i, gleich cos d baw. gleich eins beim ebenen Getriebe, wie im iibrigen auch der Ubergang vom 
Kugelradius R = 1 zum Kugelradius R-> oo zeigt. 
2. Beidem Kreuzgelenk oder der spharischen Doppelschleife! ist a = 6 = ¢ = 2/2 (Abb. 15a), 


und es gilt beilaufig bei der im Bild angegebenen Winkelzahlung 


cosdtgB =tga, (23) 
und in den Steglagen (Abb. 15b) folgt sofort 
i =i, = Le (24a) 
cos d 
Fiir a = 90°, der .,Steglage“ fiir die Abtriebskurbel (Abb. 15c) wird p = — (d + 6) = — (d 4+- 2/2), 
q =p +d =—a/2, d.h. aber nach (21a) 
i =cosd, te) 


als dem reziproken Wert zu i,, im itibrigen in Ubereinstimmung mit den durch Differentiation 
von ( gefundenen Werten. 
3. Die spharische Kreuzschleife oder der ,,Hookesche Schliissel fiihrt als das dem Kreuzgelenk 


benachbarte Getriebe! mit b = ¢ = d =a/2 gema® (21a) auf den einfachen Wert 
somit auf i, = tga = —i,. Dies ist auch leicht aus Abb. 15c, b zu erkennen, wenn man sich das 
kleine Glied jetzt als Steg vorstellt. 


4. Die spharische Schubkurbel (vgl. auch Absatz d) fiihrt mit 6 =d = /2 auf die gleichen 
Werte fiir ¢ in (25) und fiir 1, bzw. i. ; 
d) Das Extremum des Ubersetzungsverhaltnisses. 1. Bedingung fiir das Extremum. 
Nach (21a) folgt 
CL ee Dee (26) 
da sin? (p + d) 
und dieser Wert stellt beilaufig auch die bezogene Winkelbeschleunigung i’ = &,/m? dar. Soll nun t 
ein Extremum haben, so muB — formal wie bei den ebenen Getrieben — 
pes (27) 
Le Tita y 
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Abb. lia. Abb. 15b. Abb. 15c. 


Abb, 16. Abb, 17. Abb. 18. 


Abb. 15a. Die sphirische Doppelschleife, das 
Kreuzgelenk. 


Abb. 15b. Steglage der Antriebskurbel beim 
Kreuzgelenk, 


Abb. 15c. Steglage der Abtriebskurbel beim 
Abb. 19a. Kreuzgelenk. Abb. 19b. 


Abb. 16. Extremum des Ubersetzungsverhiilt- 
nisses beim spharischen Kurbeltrieb. 


Abb. 17. Zwei benachbarte Lagen H und H, des 
Relativpoles. 


Abb. 18. Spharische Polarsubnormale. 
Abb. 19a. Schwingende spharische Kurbelschleife. 
Abb. 19b. Umlaufende spharische Kurbelschleife. 

Abb, 20a. Spharische Schubkurbel. 


Abb. 20b, Extremum des Ubersetzungsverhiilt- 
nisses bei der spharischen Schubkurbel. 


Abb. 20a. Abb. 20b. 
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sein. Diese Beziehung kann man nun ohne Zuhilfenahme von spharischen Polarkoordinaten zu- 
nachst in gleicher Weise kinematisch deuten, wie an Hand von Abb. 6 geschildert wurde. Erweitert 
man analog in Abb. 13 die Koppel durch eine Schleife, welche sich mit einem in H angebrachten 
Gleitstein paart, der absolut auf dem durch den Steg gezogenen Grundkreis gefiihrt wird, so 
kann p’ als die Geschwindigkeit des Punktes H bzw. des Gleitsteines angesehen werden, welche 
tangential an den genannten Kreis verlauft. Ermittelt man dann p’ = v wie beim ebenen Getriebe 
aus Fithrungs- und Relativgeschwindigkeit, so mu auch jetzt entsprechend der GroBkreis durch 20 
und 31 — die ,,spharische K ollineationsachse“ — senkrecht zur spharischen Koppel, d. h. senkrecht 
zum GroBkreis durch 32 und 21 baw. 31 stehen, wie in Abb. 16 angedeutet. 

Mit Hilfe spharischer Polarkoordinaten ergibt sich das Folgende. Seien in Abb. 17 H und fae 
die benachbarten Lagen des Relativpoles und bilde der durch HH, gezogene GroSkreis mit dem 


»Polstrahl* 4AjH, den Winkel y, so folgt mit den Winkeln ZH, = dp, Tl =x zunachst aus dem 
spharischen Teildreieck A, Z, H nach dem Sinus-Satz 


Peis : Re: 
sin x : sin p = sin (da): sin = 
oder 

sin x = sin p da, (28) 


sowie beim ebenen Getriebe x = p da ist. Aus dem rechtwinkligen Teildreieck HZ,H, ergibt sich 
ferner nach den Neperschen Regeln 


TU IU = 
cos (3 — dp) =O & = *) ctg yp 


oder 


tex sin x 


t a= = =s ee 
By sin dp cos x: dp 


Da x klein ist, so hat man unter Einsetzen von (28) ahnlich wie in der Ebene 


sin p da sin p 
tgy=— dp a p’ : (29) 


to} 


Zeichnet man nun in Abb. 18 durch H den Tangentenkreis an die relative Polkurve und senkrecht 
dazu den Normalenkreis, so treffen diese den zum ,,Polstrahl“ p senkrecht gezogenen GroBkreis 


in T bzw. N, und der Bogen bzw. Winkel A\N stellt die spharische ,,Polarsubnormale“ s,, (und 


beilaufig A,T die spharische ,,Polarsubtangente“) dar. Wendet man auf das sphiarische recht- 
winklige Dreieck HA,N die Nepersche Regel an, so folgt 


EL =| Tt (ate = 
cos (5 — P) = tg (5 — )ete(5 — 3] 


oder unter Beachtung von (29) schlieBlich 


tgs, = Pp (30) 


gegeniiber s, = p’ in der Ebene. ; 

Im Falle eines Extremums von 7, d. h. fiir p’ = 0 ist somit s, = 0 oder: das Ubersetzungs- 
verhaltnis i hat ein Extremum, wenn die durch die Polachse 0—k, zum (relativen) 
Polkegel gezogene Normalebene mit der Stegebene (0 —10—k,) zusammenfiallt 
bzw. die Tangentialebene zur Stegebene senkrecht steht — in Ubereinstimmung mit 
dem oben gefundenen Ergebnis. ; 

Im iibrigen lat sich fiir i’ ein ahnlicher Ausdruck wie beim ebenen Getriebe in (7)* angeben. 
Ks ist 
; z. 


= (1 —icosd)t 


cos p 


l" Bi 
wenn x = 2/2 —y in Abb. 18 gesetzt wird. Denn nach den Neperschen Regeln gilt fiir Dreieck 
NAH doch ctg x tg s, = sin p, und setzt man mit tgs, = p’ den hieraus folgenden Werten fir tg x 
sowie i aus (21a) ein, so erkennt man die Form 7’ aus (26). + 

LaBt man den Kugelhalbmesser R—> 00 gehen, so gehen cos p und cos d nach eins, und man 
erhalt Gleichung (7*), wenn man beachtet, daB dort ~,, = 2/2 — x ist. 
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2. Beispiele. Beispiel 1: Beim Kreuzgelenk sind die in Abs. ¢ 2) erwahnten Stellungen fir 


a = 0 baw. « =z und a = 2/2 baw. « = 3/2 diejenigen fiir die Extrema, da die vorstehende | 
Bedingung, wie leicht aus Abb. 15a und b zu erkennen, erfiillt ist. | 

Beispiel 2: Bei der spharischen Kreuzschleife sind, wie aus Abs. c 3) zu erkennen, die Werte i 
in den Steglagen bereits die extremalen. 

Beispiel 3: Bei der spharischen Kurbelschleife’ werden b =¢ =a/2 (Abb. 19a fiir die schwin- | 
gende und Abb. 19b fiir die umlaufende Schleife, d.h. a< d bzw. a> d). Da auch hier in den 
Steglagen die Bedingung fiir das Extremum von i erfiillt ist, sind die Werte fiir i in (22a) bereits” 
die extremalen — wie bei der ebenen zentrischen Kurbelschleife. | 

Beispiel 4: Bei der spharischen Schubkurbel nach Abb. 20a (vgl. Abs. ¢ 4) — liegt die Stellung 
extremaler Ubersetzung zwischen a = 2/2 und —z/2, wie in Abb. 20b angedeutet. Die rech- 
nerische Behandlung fiihrt auf eine Gleichung achten Grades in sin?a, worauf an anderer Stelle 
ausfihrlich eingegangen wird?, wahrend sich fiir das ebene Getriebe (R—> oo) nur eine Gleichung 
dritten Grades in sin? a ergibt (weswegen auch der Winkel « im spharischen Getriebe in gleicher 
Art gezahlt wurde wie beim ebenen Getriebe). | 


4, Zusammenfassung. Nach Definition des Ubersetzungsverhaltnisses fiir ebene Getriebe, der 
Ermittlung der Bedingungen fiir ein Extremum und Anwendungen auf Gelenkgetriebe, insbesondere | 
auch auf die Viergelenkkette wurden die Ergebnisse auf die spharische ,,Viergelenkkette“ iiber- 
tragen und durch Beispiele erlautert. Im Ubrigen kénnen die Bedingungen fiir ein Ubersetzungs- | 
extremum in einem mehrgliedrigen ebenen Getriebe sinngema® in gleicher Weise auf spharische 
mehrgliedrige Getriebe ibertragen werden. | 


(Eingegangen am 24, Februar 1959.) 


Anschrift des Verfassers: Professor Dr.-Ing. W. Meyer zur Capellen, Aachen, Techn. Hochschule. 
Lehrstuhl fur Getriebelehre. 


' Vgl. FuBnote 2 von Seite 360. | 
» W. Meyer zur Capellen u. W. Rath, Kinematik der spharischen Schubkurbe!, Forsch i | 
Landes Nordrhcin-Westfalen (Opladen 1960). : ameter | 
3 Vel. FuBote 1 von Seite 354. | 
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Technische Strémungslehre 


Von Dr.-Ing. BRUNO ECK 


Neudruck der fiinften Auflage 


Mit 407 Abbildungen. X, 422 Seiten Gr.-8°. 1958. Ganzleinen DM 29,40 


AUS DEN BESPRECHUNGEN 


»»Die vorliegende finfte Auflage des bekannten Fachbuches ist ein nahezu unveranderter Nachdruck der 
erst 1954 erschienenen vierten Auflage. Die ,Technische Strémungslehre‘ von Eck wird seit ihrem ersten 
Erscheinen im Jahre 1940 vor allem von den in der Praxis stehenden Ingenieuren wegen ihrer verstand- 
lichen Darstellungsweise geschatzt. Der umfangreiche Stoff der Strémungslehre ist zweckmaBig gegliedert 
und besonders auf die Gebiete beschrankt, die fiir den Maschinenbau von Bedeutung sind... 

Der Verfasser verzichtet bei der Behandlung dieses umfangreichen Stoffes auf verwickelte Ableitungen und 
stellt die wichtigsten Erscheinungsformen und Gesetze der Strémungslehre knapp und unter Verwendung 
einpragsamen Bildmaterials dar. Dem an speziellen Problemen interessierten Leser steht ein umfangreiches 


Literaturverzeichnis zur Verfiigung.‘* Die Technik 


Technische Hydro- und Aeromechanik 


Von Dr.-Ing. habil. WALTHER KAUFMANN, em. o. Professor der Mechanik an der Technischen 
Hochschule Miinchen 


Zweite, verbesserte und erganzte Auflage 


Mit 266 Abbildungen. VIII, 386 Seiten Gr.-8°. 1958. Ganzleinen DM 37,50 


AUS DEN BESPRECHUNGEN 


»,Dieses nach wenigen Jahren bereits in 2. Auflage herausgegebene Buch kann als zusammenfassendes 
Lehrbuch alles dessen bezeichnet werden, was man heute unter Strémungslehre versteht, zunéchst ohne 
Riicksicht auf den Unterschied von Gasen und Fliissigkeiten. Nachdem in zwei kleineren Abschnitten kurz 
die besonderen Eigenschaften der strémenden Medien und die Hydro- bzw. Aerostatik behandelt werden, 
ist der weitaus gréBte Teil der Bewegungslehre gewidmet. Nach der Stromfadentheorie und der Darstellung 
der besonderen Energieverluste werden die Theorien fiir zwei- und dreidimensionale Strémungen, ideelle 
und mit Zahigkeit behaftete Medien behandelt. SchlieBlich sind kurze Abschnitte der Grundwasserbewe- 
gung, der Schmiermittelreibung, dem Tragfliigel und seiner Anwendung im Gitter gewidmet. Den Ab- 
schlu8 bildet ein Abschnitt iiber Gasdynamik. 

Neben dem Studierenden wird es auch dem Konstrukteur von Strémungsmaschinen, dem Wasserbau- 
Ingenieur usw. viel Interessantes und Grundsitzliches bieten, auch wenn Erfahrungswerte eher sparlich 


angegeben sind.“ Schweizer Archiv f. angew. Wissenschaft u. Technik 
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Soeben erschien: 


Ingenieur-Mathematik 


Von Dr. R. SAUER, Professor an der Technischen Hochschule Miinchen 


Erster Band: Differential- und Integralrechnung 


Mit 178 Abbildungen. VIII, 304 Seiten Gr.-8°. 1959. Ganzleinen DM 24,— 


INHALTSVERZEICHNIS 


Differential- und Integralrechnung fiir Funktionen von einer Verinderlichen: Reelle Zahlen. Funktionen 
von einer Veranderlichen; Stetigkeit. Spezielle Funktionen; Kurvendiskussion. Lineare analytische 
Geomettie der Ebene. Analytische Geometrie der Kegelschnitte. Grenzwert. Grundziige der Differential- 
rechnung. Mittelwertsatze der Differentialrechnung. Numerische und graphische Differentiation; Inter- 
polation. Grundziige der Integralrechnung. Graphische und numerische Integration. Logarithmus und 
Exponentialfunktion; Hyperbelfunktionen. Rechenschieber und logarithmische Papiere. Unendliche 
Reihen. TAYLOR-Entwicklung von Funktionen in Potenzreihen. Anwendungen der TAYLOR-Ent- 
wicklung fiir das numerische Rechnen. Komplexe Zahlen. Elementar integrierbare Funktionenklassen. 


Differentialgeometrie der ebenen Kurven. Anwendung der Differentialgeometrie auf die Getriebelehre. 


Differential- und Integralrechnung fiir Funktionen von mehreren Veranderlichen: Funktionen von meh- 
reren Veranderlichen. Graphische Darstellung von Funktionen mehrerer Veranderlicher (Nomographie). 
Vektorrechnung. Determinanten und Systeme linearer Gleichungen. Lineare analytische Geometrie des 
Raumes. Analytische Geometrie der Flachen 2.:Ordnung. Grundziige der Differentialrechnung bei 
Funktionen von mebreren Veranderlichen. TAYLOR-Entwicklung fiir Funktionen von mehreren Ver- 
imderlichen, Anwendung der TAYLOR-Entwicklung bei mehreren Veradnderlichen fir das numerische 
Rechnen. Allgemeine Abbildungen und allgemeine Koordinatensysteme. Integraldarstellung von Funk- 
tionen. Mehrfache Integrale. Differentialgeometrie der Kurven und Flachen im Raum. Anwendungen 

der Integralrechnung in der Mechanik - Anhang - Sachverzeichnis. 


ZUR INFORMATION 


Das Buch richtet sich an Studierende der Ingenieurwissenschaften und der Physik und gibt eine Einfiihrung 
in die Mathematik in einem Umfang, der etwa den beiden ersten Semestern der mathematischen Kurs- 
vorlesungen an den Technischen Hochschulen entspricht (Differential- und Integralrechnung samt Grund- 
ziigen der Vektorrechnung und der analytischen und Differentialgeometrie). Es wird durchwegs versucht, 
die Begriffe und Methoden in einer Weise zu entwickeln, die der aufs Anschauliche gerichteten Denkweise 
des Ingenieurs angemessen ist. Numerische und graphische Methoden nehmen einen verhiltnismaBig 
breiten Raum ein. Auch eine Einfiihrung in die Nomographie ist eingefiigt sowie ein Exkurs in die 

Getriebelehre. 


Ein spater folgender zweiter Band wird sich mit den Differentialgleichungen und der Funktionentheorie 


befassen. 
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